AD' A 1 3 1  425  NUCLEAR  WARFARE  WATER  CONT  AMINAT  I ONf U )  SCIENCE 

APPLICATIONS  INC  SCHAUMBURG  II  JC  PHILLIPS  ET  AL 
01  MAY  82  DNA-TR-81-127  DNAOO T - Q 1 -C - 02 34 


"3 


I 


UNCLASSIFIED 


F/G  15/6 


NL 


OTIC  FILE  copy  AO  a  l  8  l  4 


DNA-TR-81  -127 


NUCLEAR  WARFARE  WATER 
CONTAMINATION 


Science  Applications,  Inc 
Chicago  Office 

One  Woodfield  Place  Building 
1701  E.  Woodfield  Road 
Schaumburg,  Illinois  60195 


1  May  1982 


Technical  Report 


CONTRACT  No.  DN A  001-81-C-0234 


APPROVED  FOR  PUBLIC  RELEASE; 
DISTRIBUTION  UNLIMITED. 


THIS  WORK  WAS  SPONSORED  BY  THE  DEFENSE  NUCLEAR  AGENCY 
UNDER  RDT&E  RMSS  CODE  X384081469  Q22QAXNA00001  H2590D. 


Prepared  for 
Director 

DEFENSE  NUCLEAR  AGENCY 
Washington,  DC  20305 


DTIC 


AUG  i  0  1393 


i; 


0  r-r 


i. 


Destroy  this  report  when  it  is  no  longer 
needed.  Do  not  return  to  sender. 

PLEASE  NOTIFY  THE  DEFENSE  NUCLEAR  AGENCY, 
ATTN:  STT I ,  WASHINGTON,  D.C.  20305,  IF 
YOUR  ADDRESS  IS  INCORRECT,  IF  YOU  WISH  TO 
BE  DELETED  FROM  THE  DISTRIBUTION  LIST,  OR 
IF  THE  ADDRESSEE  IS  NO  LONGER  EMPLOYED  BY 
YOUR  ORGANIZATION. 


UNCLASSIFIED 


REPORT  DOCUMENTATION  PAGE 


4  TITLE  (and  Subtitle) 


NUCLEAR  WARFARE  WATER  CONTAMINATION 


RF  AD  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


RECIPIENT’S  C»TALOC  NUMBER 


J  type  of  report  a  period  covered 


Technical  Report 


6  PERFORMING  ORG  REPORT  NUMBER 


7  AUTHORS, 


J.  C.  Phillips,  J.  A.  Roberts,  R.  H.  Sievers, 
and  R.  Gminder 


■  CONTRACT  or  GRANT  NUMBERFO 


DNA  001-81-C-0234 


9  PfcRtORMINL*  ORGANIZATION  NAMF  AnD  A  D  RTSS 

Science  Applications,  Inc.,  Chicaoo  Office 


»□  PROGRAM  ELEMENT.  PROJECT  TASK 
AREA  A  WORK  UNIT  NUMBERS 


One  Woodfield  Place  Bldg.,  1701  E.  Woodfield  Road  Task  Q22QAXNA-00001 
Srhaiimhnrn .  Tllinrnc. 


Schaumburg,  Illinois,  60195 

M  CONTHOLl'NG  OFFICE  NAME  AND  AGORtSS 


12  REPORT  DATE 

1  May  1982 


Director  L  May 

Defense  Nuclear  Agency  '3  number  of  paces 

_ Washington.  D.  C. .  20305  _ 222 _ 

monitoring  AGEN-Tv  N  AM  E  A  A  DON  ?.  v'.<  i  /  jiNrrcnf  trorn  ('on  rr<tj/ mj  Oil,,  ?  >  !  I  S  b  t  C.  U  R1  T  Y  CLASS  (o  t  this  report) 

|  UNCLASSIFIED 

j  15«  ~DECL  ASSIFIC  ATION  OOWNGRAOlN^G 
SCHEDULE 


I 


<6  Distribution  statement  (oi  th i.  RePc 


Approved  for  public  release;  distribution  unlimited. 


17  D|STRi0lt,ON  ST  ATEMtNT  ( o(  the  abttmet  entered  tn  Mock  10.  it  different  trom  Report) 


10  SuPPL  EMLN  T  AR'i  NOtES 


This  work  was  sponsored  by  the  Defense  Nuclear  Agency  Under  RDT£E  RMSS  Code 
-  X 384081469  Q22QAXNAOOO01  H2590D. 


T9  KEY  WORDS  Canrrnue  on  reverse  eiJe  if  net-*  •.  :.  «ir  ■»  end  identify  hy  block  number > 


Water 

Pal  lout 

Contamination 


Nuclear  Warfare 
Water  Purification 
Water  Decontamination 


Radiation  Monitoring 


i  20  AflSTH  A.‘  T  Cmlinu*  on  revert*  tide  .1  ne>»>.\*r\  ami  toentifv  hi  block  number 

The  fallout  contamination  of  watersheds  and  water  supplies  resulting  from 
surface  burst  nuclear  weapons  is  sufficiently  high  to  require  the  use 
of  water  purification  equipment  to  produce  potable  water  that  meets  the 
current  radiological  water  quality  standards.  Problems  with  existing 
radiation  detection  and  monitoring  equipment  are  discussed.  A  water 
contamination  model  is  described.  . 


1473  EClTiCN  OF  1  NOV  AS  i$  OBSOLETE 


UNCLASSIFIED 


SECURITY  Classification  OF  Tnis  PAGE  I  »»>*•«  />**•  Fmered) 


TABLE  OF  CONTENTS 


Section  Page 

1.  INTRODUCTION .  5 

2.  THREAT  CHARACTERIZATION .  7 

2.1  Introduction .  7 

2.2  Contamination  Sources .  7 

2.3  Water  Sources .  11 

2.4  Water  Contamination  Processes .  14 

2.4.1  Fallout  Production,  Transport,  and  Deposition .  15 

2.4.2  Radioactive  Decay .  28 

2.4.3  Dissolution  of  Radionuclides .  28 

2.4.4  Transport  of  Radionuclides  by  Precipitation 

Runoff .  30 

2.4.5  Aquatic  Mixing  and  Transport .  33 

2.5  Water  Contamination  Model .  34 

3.  THREAT  RECOGNITION .  45 

3.1  Introduction .  45 

3.2  Water  Quality  Standards .  45 

3.3  Radiation  Monitoring  Equipment .  49 

3.4  Radiation  Monitoring  Procedures .  52 

4.  THREAT  COUNTERMEASURES .  55 

4.1  Introduction .  55 

4.2  Water  Purification  Equipment .  55 

4.3  Field  Operations  Policy  and  Procedures .  57 

5.  CONCLUSIONS .  59 

6.  REFERENCES .  60 

APPEND  I C I E  S 

A.  Water  Source  Information .  A-l 

B.  WSWCM  -  Watershed  Water  Contamination  Model .  B-l 


LIST  OF  ILLUSTRATIONS 


Figure 


Page 


1  Idealized  Fallout  Pattern  Parameters .  9 

2  Scenario  Area .  13 

3  1-135  Water  Contamination  -  Problem  1 .  36 

4  Sr-90,  Y-90  Water  Contamination  -  Problem  1 .  37 

5  Te-131m,  1-131  Water  Contamination  -  Problem  1 .  38 

6  Total  (40  Radionuclides)  Water  Contamination  -  Problem  1...  39 

7  1-135  Water  Contamination  -  Problem  2 .  41 

8  Sr-90,  Y-90  Water  Contamination  -  Problem  2 .  42 

9  Te-131m,  1-131  Water  Contamination  -  Problem  2 .  43 

10  Total  (40  Radionuclides)  Water  Contamination  -  Problem  2...  44 

11  TB  MED  229  Requirements .  46 

12  STANAG  2136  (MED)  -  Requirements .  48 


2 


1. 


LIST  OF  TABLES 


Tabl e  Page 

1  Characteristics  of  Idealized  Fallout  Patterns .  10 

2  Radionuclide  Ground  Concentrations .  17 

3  Percent  of  Organ  Dose  Commitment  for  Specific 

Radionuclides .  21 

4  Percent  of  External  Dose  Rate  for  Specific 

Radionuclides .  25 

5  Selected  Distribution  Coefficients .  31 

6  Responses  to  Request-For-Information .  51 


3 


SECTION  1 


INTRODUCTION 


In  the  event  of  nuclear  warfare,  the  raw  water  supplies  available 
to  U.  S.  Army  field  units  could  become  radiologically  contaminated. 
Water  is  essential  to  the  health  and  hygiene  of  personnel; 
radiologically  contaminated  water  poses  a  physiological,  and 
pyschological ,  threat  to  personnel  who  are  expected  to  perform 
effectively  in  a  radiologically  contaminated  environment.  To  address 
the  potential  threat  of  contaminated  water,  the  U.  S.  Army  has 
guidelines  and  standards  for  water  quality,  equipment  and  procedures 
for  the  radiological  monitoring  of  water,  and  equipment  and  procedures 
for  the  purification  of  contaminated  water. 

This  assessment  of  the  nuclear  warfare  water  contamination  threat 
is  principally  concerned  with  the  equipment  and  procedures  for  the 
monitoring  and  purification  of  radiologically  contaminated  water.  The 
assessment  does  not  explicitly  address  water  quality  standards; 
however,  the  standards  are  discussed  in  terms  of  the  technical 
performance  criteria  that  they  provide  for  the  water  monitoring  and 
purification  equipment. 

The  objective  of  this  water  contamination  threat  assessment  are: 
(1)  to  characterize  the  contamination  threat,  (2)  to  identify  field 
methods  for  recognizing  the  contamination  threat,  and  (3)  to  indicate 
the  effectiveness  of  possible  threat  countermeasures.  The  assessment 
is  focused  on  the  water  contamination  threat  associated  with  nuclear 
warfare  in  Europe  and  is  primarily  concerned  with  the  contamination  of 
rivers  and  streams  by  the  fission  products  produced  by  a  nuclear 
weapon.  With  respect  to  the  degree  of  detail  of  technical  analysis, 
the  assessment  is  limited  to  a  scoping  level  with  the  intent  of 
separating  problems  from  non-problems. 
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Section  2  of  this  report  addresses  the  subject  of  water 
contamination  threat  characterization.  It  discusses  water  sources, 
contamination  sources,  and  the  water  contamination  processes;  a  model 
used  to  determine  the  time-dependent  concentration  of  fission  product 
radionuclides  dissolved  in  water  is  presented.  The  subject  of  water 
contamination  threat  recognition  is  addressed  in  Section  3.  The 
ability  to  measure  radiological  water  contamination  in  order  to 
satisfy  specified  water  quality  requirements  is  discussed.  Water 
contamination  threat  countermeasures  are  discussed  in  Section  4  of 
this  report.  The  effectiveness  of  water  purification  equipment  and 
procedures  are  addressed.  Section  5  of  this  report  summarizes  the 
conclusions  that  were  reached  during  this  assessment.  Materials  cited 
as  reference  within  the  report  are  identified  in  Section  6. 
Supporting  material  is  provided  in  two  appendices:  Appendix  A  is  a 
report  on  the  water  sources  and  characteristics  within  the  geographic 
area  selected  for  study;  Appendix  B  is  a  report  that  documents  the 
water  contamination  model  developed  for  this  assessment. 
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SECTION  2 


THREAT  CHARACTERIZATION 


2.1  Introduction 


Threat  characterization  refers  to  a  determination  of  the 
magnitude  and  duration  of  the  radiological  water  contamination  th^t 
results  from  a  specific  set  of  circumstances  and  is  expressed  in  terms 
of  the  time-dependent  activity  concentration  of  radionuclides  in 
water.  For  a  given  contamination  source  and  water  source,  the 
radiological  water  contamination  is  determined  by  considering  the 
dissolution,  mixing,  and  transport  processes  that  affect  the 
contaminant  material.  A  simplified  computer  model  has  been  developed 
to  characterize  the  water  contamination  of  rivers  and  streams  by  the 
fission  products  produced  by  a  nuclear  weapon. 

2.2  Contamination  Sources 


Within  the  broad  concept  of  nuclear  warfare,  three  general  types 
of  radiological  contamination  sources  can  be  considered:  nuclear 

weapons,  radiological  warfare  agents,  and  nuclear  sabotage  or 
terrorist  actions.  Although  all  three  sources  were  initially 
considered,  this  assessment  quickly  focused  on  the  nuclear  weapon 
contamination  source  as  the  most  significant  and  relevant  to  the 
overall  objective  of  the  assessment. 

A  nuclear  weapon  employed  in  the  surface  burst  mode,  either 
intentionally  or  accidently  due  to  a  firing  system  error  or  failure, 
can  produce  radioactive  fallout  over  a  widespread  area.  The  intensity 
and  extent  of  the  fallout  area  is  principally  determined  by  the 
fission  yield  of  the  weapon  and  the  meteorological  conditions  that 
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occur  during  the  fallout  deposition  period.*  An  idealized  fallout 
pattern  (see  Figure  1),  based  on  the  parameters  given  in  "The  Effects 
of  Nuclear  Weapons,"  provides  a  convenient  way  of  portraying  a  fallout 
area.^**  Table  1  shows  some  of  the  characteristics  of  idealized 
fallout  patterns.  Note  that  a  10-KT  nuclear  weapon  can  contaminate  an 
area  of  approximately  9000  Km  with  a  reference  dose  rate  equal  to,  or 
exceeding,  1  R/Hr  at  H+L  hour.  Clearly,  a  nuclear  weapon  has  the 
potential  for  contaminating  a  large  area  that  encompasses  hundreds  of 
individual  watersheds  that  could  service  U.  S.  Army  water  supply 
points.*** 


Radiological  warfare  was  discussed  in  the  popular  press  in  the 

1950's  and  the  1960's.  One  concept  involved  using  special  materials 

within  a  nuclear  weapon  to  increase  the  amount  of  radioactive  material 

produced.  Another  concept  involved  the  localized  dispersal  or 

emplacement  of  gross  fission  products  or  selected  radionuclides. 

Although  scientifically  feasible,  such  concepts  have  considerable 

engineering  and  logistical  difficulties;  in  addition,  the  military 

utility  of  these  concepts  is  not  clear,  nor  is  it  clear  how  the 

concepts  would  fit  within  a  nuclear  warfare  strategy  that  emphasizes 

the  offense.  Furthermore,  based  on  discussions  with  representatives  of 

the  Defense  Intelligence  Agency,  it  was  learned  that  no  specific  or 

(2) 

credible  radiological  warfare  threat  is  currently  projected.  There¬ 
fore,  no  further  consideration  was  given  to  the  use  of  radiological 
warfare  agents  as  a  possible  water  contamination  source. 


r  In  general ,  tactical  or  low-yield  nuclear  weapons  generate  their 
yield  by  the  fission  process;  strategic  or  high-yield  nuclear 
weapons  use  both  the  fission  and  the  fusion  processes.  If  part  of 
the  weapon  yield  is  obtained  from  the  fusion  process,  the  fallout 
area  will  be  affected  since  the  fusion  process  does  not  produce 
radioactive  material  to  the  extent  that  the  fission  process  does. 

**  The  number  in  the  parentheses  denotes  a  reference  that  is 
identified  in  Section  6. 

***Typical ly,  watersheds  in  the  European  area  considered  in  this 
assessment  had  areas  in  the  range  of  10  Km  to  50  Km  . 
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DOWNWIND  DISTANCE  FOR  I  RAD/HR 


GROUND  ZERO 


Figure  9  93.  Illustration  of  idealized  unit  time  dose-rate  pattern  for  early  fallout  from  a 
surface  hurst  (The  contour  dimensions  arc  indicated  for  a  dose  rate  of  1 
rad'hr  ) 
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faliout  pattern 


Th<?-,e  characteristics  are  based  on  the  fallout  pattern  parameters  given  in  Glasstone,  S.  and  P.  J.  Dolan,  "The  Effects  of  Nuclear  Weapons 
U.  5.  Departments  of  Defense  and  Energy,  Washington,  D.  C.,  1977,  page  430. 


Nuclear  sabotage  or  terrorist  action  has  been  the  subject  of  much 
discussion  for  the  past  several  years.  Studies  of  postulated 
accidents  at  nuclear  facilities  tend  to  focus  on  groundwater 
contamination  rather  than  surface  water  contamination;  the  resulting 
level  of  water  contamination  and  the  time  scale  on  which  the 
contamination  occurs  indicate  that  such  events  are  long  term 
environmental  problems  rather  than  problems  of  concern  in  a  nuclear 
warfare  environment.^  Conceivably,  material  stolen  from  a  nuclear 
facility  could  be  used  to  cause  the  contamination  of  a  water  supply 
source,  for  example,  a  spent  reactor  fuel  assembly,  stolen  from  a 
power  plant  or  during  shipment,  could  be  dumped  into  a  water  reservoir 
or  lake.  The  potential  level  of  water  contamination  associated  with 
such  an  event  could  be  rather  high,  on  the  order  of  300,000  p C i / 
(gross  fission  products);*  however,  the  extent  of  the  water 
contamination  could  be  quite  localized.  Since  nuclear  weapons  have 
the  potential  for  equally  severe  levels  of  water  contamination  over 
wide  scale  areas,**  there  appeared  to  be  no  particular  benefit  to  any 
further  study  of  hypothetical  nuclear  sabotage  or  terrorist  actions. 

2.3  Water  Sources 

In  general,  the  variety  of  raw  water  sources  includes  springs  and 
wells,  rivers  and  stream,  lakes  and  ponds,  and  ocean  bays  and  harbors. 
However,  within  the  context  of  nuclear  warfare  in  Europe,  the 
imposition  of  strict  radiological  defense  measures  will  limit  the 
selection  of  raw  water  sources  to  those  sources  that  have  sufficient 

*  A  spent  fuel  assembly  contains  roughly  600,000  Ci  of  gross  fission 
products  after  a  cooling  time  of  one  year. (4)  It  is  estimated  that 

0.1%  of  the  fission  products  would  be  released  into  the  water  (see 
Reference  5).  A  small  reservoir  or  lake  having  an  average  depth  of 
10  ft.  and  a  water  surface  area  of  0.25  mi?  contains  2  X  10  9;  of 
water. 

**In  Section  2.5,  a  model  is  presented  that  correlates  the  level  of 
water  contamination  with  the  level  of  fallout  contamination.  For  a 
fallout  level  of  1  R/Hr  at  H+l  hour,  the  water  contamination  is 
initially  around  10 8  pCi/c  and  drops  to  about  10 ^  pCi /  ;  in  about  10 
days. 


capacity  to  support  the  operation  of  an  engineer  water  supply  point. 
In  addition,  only  those  raw  water  sources  that  are  locally  available 
to  the  field  forces  are  of  interest. 


For  this  water  contamination  threat  assessment,  the  focus  was  on 
a  preselected  nuclear  warfare  scenario  area  in  West  Germany  bounded  by 
Marburg,  Giessen,  and  Frankfurt  am  Main  on  the  west,  and  the  Fulda 
River  Valley  on  the  east  (See  Figure  2).  Descriptive  information  on 
the  water  source  characteristics  of  this  area  are  contained  in 
Appendix  A. 

Within  the  scenario  area,  the  primary  raw  water  sources  are 

rivers  and  streams.  Based  on  a  map  study  (scale  1:50000)  of  the  area 

and  in  accordance  with  U.  S.  Army  Engineer  School  doctrine,  thirty 

water  supply  points  were  identified  and  the  specific  watershed  area 

for  each  water  supply  point  was  delineated.  The  area  of  the 

2  2 

watersheds  ranged  from  8  Km  to  60  Km  ,  with  an  average  value  of  20 
2 

Km  .  For  a  typical  watershed,  the  stream  or  river  area  was  about  0.1% 
of  the  total  watershed  area  and  the  estimated  average  stream  velocity 

at  the  watershed  exit  was  0.6  m/s.  In  the  absence  of  precipitation 

runoff,  the  typical  volumetric  flow  rate  ranged  from  38  £  /s  (36,000 
gph)  to  100  c.  /s  (95,000  gph),  with  the  lower  flow  rate  occurring  in 
July,  August,  and  September  and  the  higher  flow  rate  occurring  in 
February,  March  and  April. 

Precipitation  occurs  frequently  (on  almost  half  the  days)  within 
the  scenario  area  with  the  largest  amounts  of  rainfall  typically 
occurring  in  June,  July,  and  August.  When  rainfall  sufficient  to 
promote  surface  runoff  occurs,  the  volume  of  water  associated  with  the 
precipitation  runoff  can  be  a  factor  of  2  to  5  times  more  than  the 
volume  of  water  present  on  the  watershed  during  dry  conditions.  This 
volume  of  surface  runoff  water  will  flow  off  the  watershed  into  the 
stream  or  river  over  a  period  of  four  days. 
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Figure  2.  Scenario  area 


2.4  Water  Contamination  Processes 


The  radiological  contamination  of  rivers  and  streams  that  would 
occur  due  to  the  radioactive  fallout  from  a  nuclear  weapon  detonation 
is  the  result  of  several  different  processes.  Fallout  material  is 
produced,  transported  by  the  winds,  and  then  deposited  on  the 
watershed  area.  The  fallout  material  that  is  deposited  on  the  water 
surface  will  experience  dissolution,  mixing,  and  transport,  along  with 
radioactive  decay.  The  fallout  material  that  is  deposited  on  the  land 
surface  will  initially  experience  only  radioactive  decay,  but 
subsequent  precipitation  could  dissolve  the  material  and  transport  it 
to  the  watershed  stream  or  river  where  it  mixes  with  the  channel  water 
and  is  transported  out  of  the  watershed  area. 

The  processes  that  must  be  considered  to  determine  the 
radiological  water  contamination  thus  include: 

•  fallout  production,  transport,  and  deposition, 

•  radioactive  decay, 

•  dissolution  of  radionuclides  in  fallout  material, 

t  transport  of  radionuclides  by  precipitation  runoff,  and 

•  aquatic  mixing  and  transport. 

These  water  contamination  processes,  with  the  exception  of  the  fallout 
production,  transport,  and  deposition  process,  have  been  incorporated 
into  a  simplified  computer  model  developed  for  this  assessment.  A 
detailed  description  of  the  water  contamination  model  is  provided  in 
Appendix  B. 

The  following  material  provides  a  discussion  of  the  water 
contamination  processes  identified  above,  with  emphasis  on  the 
approach  used  to  incorporate  the  processes  into  the  model.  It  should 
be  noted  that  this  water  contamination  model  addresses  the  water 
concentration  of  soluble  radioactive  material;  insoluble  radioactive 
material  is  not  considered  since  such  material  can  be  easily  removed 
from  the  water  by  simple  filtration  methods. 
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2.4.1  Fallout  Production,  Transport,  and  Deposition 


The  process  of  fallout  production,  transport,  and  deposition 
determines  the  initial  surface  contamination  of  the  watershed. 
Typically,  the  output  of  a  fallout  model  specifies  the  unit-time 
reference  dose  rate  (e.g.,  10  R/Hr  at  H+l  hour)  within  the 
fallout  field.  For  the  purposes  of  the  water  contamination  model,  it 
is  necessary  to  know  the  specific  radionuclide  ground  concentrations 
(e.g.,  Ci/Km  at  H+l  hour  for  Sr-90)  that  correlate  with  the  unit-time 
reference  dose  rate.  Since  the  prime  concern  is  with  the  radionuclide 
composition  of  the  fallout,  rather  than  the  specific  dimensional 
characteristics  and  orientation  of  the  fallout  pattern  dose  rate 
contours,  it  is  appropriate  to  use  idealized  fallout  patterns,  like 
those  discussed  in  Section  2.2  and  Reference  1,  to  determine  unit-time 
reference  dose  rates  and  then  employ  a  separate  procedure  to  convert 
the  dose  rates  to  radionuclide  ground  concentrations.* 

A  methodology  to  correlate  radionuclide  ground  concentrations 
with  fallout  radiation  dose  rates,  recently  developed  by  SAI  for  use 
on  DNA's  Nuclear  Test  Personnel  Review  (NTPR)  program  is  contained  in 
a  computer  code  named  FIIDOS.^  The  radionuclide  inventory  and 
associated  gamma  radiation  emission  spectrum  are  determined  for  a 
mixture  of  fallout  material  that  includes  fractionated  fission 
products,  light  elements  which  were  made  radioactive  by  neutron 
capture  or  activation,  and  heavy  elements  or  actinides  which  were  part 
of  the  nuclear  weapon  or  resulted  from  neutron  capture  by  weapon 


♦Initial ly,  consideration  was  given  to  the  use  of  the  DELFIC  fallout 
code  (Reference  6)  to  calculate  radionuclide  ground  concentrations 
directly.  However,  since  DELFIC  outputs  mass  chain  depositions 
rather  than  specific  radionuclide  depositions,  the  code  would  have 
required  modifications.  In  addition,  the  use  of  a  sophisticated 
fallout  code  like  DELFIC  could  have  obscured  the  important 
relationship  between  radionuclide  ground  concentrations  and  fallout 
dose  rate  contours  by  focusing  too  much  attention  on  the  mechanics  of 
the  fallct  process  and  the  details  of  the  fallout  pattern. 
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material.  By  using  gamma  radiation  exposure  factors  to  determine  the 
above-ground  dose  rate  that  results  from  the  calculated  gamma 
radiation  emission  spectrum  in  a  planar  source  geometry,  the 
relationship  between  the  dose  rate  and  the  surface-distributed 

radionuclide  inventory  is  established.  It  should  be  noted  that  to  use 
FIIDOS  the  user  must  provide  fallout  material  characterization  data  to 
treat  the  fractionation  of  fission  products  and  weapon  test  data  to 

specify  the  inventory  of  activated  light  elements  and  actinides. 

The  fallout  radiation  dose  rate-radionucl ide  ground  concentration 
correlation  developed  for  this  assessment  is  based  on  a  source 

consisting  of  unfractionated  U-235  fission  products.  It  is  recognized 

that  because  of  fission  product  fractionation  the  composition  of  the 
fission  product  material  within  a  fallout  field  is  not  uniform  and  the 
relative  abundance  of  a  radionuclide  in  a  fallout  field  is  not  the 
same  as  that  observed  in  laboratory  experiments  or  calculated  by  basic 
fission  product  buildup  and  depletion  computer  codes.  However,  the 
inclusion  of  fractionation  effects  would  represent  an  added  complexity 
that  is  not  necessary,  or  appropriate,  for  the  purposes  of  this 
assessment  and  could  impact  the  results  of  the  assessment  in  a 
non-conservative  manner  by  causing  some  important  elements  (e.g., 
iodine  and  strontium)  to  be  partially  depleted  from  the  radionuclide 
inventory.  The  source  term  used  for  the  correlation  does  not  include 
activated  light  elements  or  actinides;  although  such  radionuclides 
generally  account  for  less  than  10T  of  the  radioactivity  in  the 
fallout  material,  there  are  some  radionuclides  (e.g.,  isotopes  of 
plutonium)  that  could  possibly  be  of  significance  to  radiological  dose 
estimates.  However,  the  inclusion  of  such  elements  in  the  source  term 
would  require  gross  assumptions  about  nuclear  weapon  design  or  access 
to  classified  design  and  intelligence  information. 

Table  2  lists  the  radionuclide  ground  concentrations  (Ci/Km  at 
H+l  hour)  for  a  mixture  of  unfractionated  U-235  fission  products  that 
is  referenced  to  a  fallout  deposition  contour  that  has  an  above-ground 
external  radiation  exposure  rate  of  1  R/Hr  at  H+l  hour.  For 
completeness,  the  table  includes  all  the  radionuclides,  with  the 
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Table  2.  Radionuclide  ground  concentra ti ons 


(cont*  d) 


Activity 
Rad i  on uc  1  t  de  {C  i  /  Mi  -  ) 


Pd- 109m 

3.01  -  01 

I'd-  1 1  ] 

2.  36  *  01 

Pd-  ?  1  Ir.t 

8.26  -  0. 

Pd- 1  i  r 

2.2]  -  00 

Pd-  t  1  3 

1.57  -  U9 

Pd- 1  M 

3.0!  -  06 

A'j- 1 0'tr 

1.13  *  01 

A.;-  i  in 

7.01  -  03 

Ad- 1 1 i 

7.60  -  m 

Ad-  i  1  liv 

3.46  '  Cl 

Aq-  1  17 

4.4M  -  01 

Ac- !  1 5 

6. 13  -  03 

A.-Il  it- 

6.6“  -  Id 

A.:- 11.) 

.1.10  - 

Ac-  1  H 

1  .  S:  *  0 1 

Ac  - 116 

9.44  -  06 

Ad-n; 

1.67  -  1, 

Cd-  111: 

2.%  -  O'- 

Cd- 1 1  3- 

i.-io  -  o<: 

Cd-  !  1  5 

6.1:-  -  01 

Cd  -  1  1  S,» 

4.14  -  0:  'i 

Cd- !  1  7 

7.07  -  Co 

Cd- 1 1 7t 

3.07  -  00 

Cd- 1  in 

7. IF  ♦  01 

Cd-  •  i 

1.67  . 

Cii- :  id¬ 

ea-  ... 

i.  •  ■  -  i. 

In-  Ml 

-  ii 

In-  i  p 

1  IT 

In-;  if. 

7. 4  5  -  r.i. 

It.-  !  ! 

1  _ 

I  n  -  11 

7.64  -  DO 

I n-  i  i<i 

7.  38  *  01 

In  -  1  1 1 

9.40  -  01 

In- 1  Id- 

1.66  *  01 

In-  170 

1.67  -  \7 

Act i vi ty 
Radionucl  ide  [Ci/Xrr- ) 

In- 120m  7.7s  -  )9 

Acti vi ty 

Radipnucl_i<fe  [CijXn-  ) 

Tc- ! 25:n  3.19  -  07 

In-]71i:i 

5.57  -  04 

Tc- 12  7 

1.09  -  01 

I  n  - 1 7  3'" 

3.73  -  70 

Tt—  1 27m 

9.89  -  05 

Sn- 7  J  9- 

9. 77  -  06 

Tt-129 

9.07  +  01 

Sn- 171 

1. 66  -  00 

To- 129m 

7.02  -  02 

Sr.- 1.' Ir¬ 

4.17  -  Or 

To- 131 

2.61  +  03 

on- 17  i 

9.  14  -  03 

Tc—  1  31 m 

1.31  *  01 

Sn- 1.3 

1.67  *  01 

Te- 1 32 

7.11  +  Cl 

6n- 176 

1. 01  -  0! 

Te-133 

1.18  ♦  03 

Sn  -  176.- 

5.64  -  00 

Te- 1 33m 

1.24  4  03 

Sn- 176 

1.98  -  07 

Te- 1  34 

4.03  *  03 

Sn- 1 7/ 

8.47  *  <M 

1-128 

8.42  -  04 

Sn-177m 

7.00  -  07 

1-179 

2.7!  -  )C 

Sn  - 17; 

3.64  *  07 

i-no 

6.56  -  03 

Sn- 17'. 

1.30+  01 

I  - 1 30m 

7.92  -  03 

Sn- 1 7  \,n 

6.99  .  04 

I- 131 

9.54  -  00 

Sn- 130 

3.43  -  01 

i  - 1 37 

3.03  ♦  01 

Sn-  1 ;;; 

r  •  3 

1-133 

3.03  4  02 

Sn  - 1 

.  '  a  -4  -  2  3 

1-134 

4.36  4  03 

'!■-  17. 

7.33  -  06 

1-134::. 

1.20  -  Cl 

i77>! 

6.  71  -  0° 

1-135 

9.86  4  02 

St  - 1  7  4 

1 . 1 1  -  05 

1-136 

1 .  94  -  08 

SI  -  174  .-i 

1.54  -  17 

I  - 1 36- 

4.  ■”  -  18 

Sb- 176 

C.oo  .  ,  . 

Cs-  1  34 

1.49  -  06 

St-  176 

7  ,c-7  -  0  I 

tS-1  3.4- 

6.02  -  03 

Si-1.  >■- 

2.67  -  01 

Cs-1  <6 

2. 38  -  09 

SI- 17  7 

2.16  -  0: 

C;.-  '■  3::r 

7.67  -  a; 

51—173 

5.71  -  {10 

Cs-136 

6.09  .  02 

55- 1 7.—. 

4.  31  +  07 

Cs- 7  37 

:  -  O' 

Sb- 1 79 

7.46  *  o7 

C  -. 

6.53  )  03 

Sb  - 1  .v.' 

1.16  *  01 

C  ,  -  i  .  .!• 

3.19  -  03 

St,  1  .ill" 

1.07  ‘  03 

C'-.-pd 

5.97  4  02 

s(,-ni 

1.6.)  *  '1 

Cr-UO 

4.43  -  12 

Sb-I 37 

1.97  -  04 

64-1  .V- 

6.81  -  06 

Sb-13?n 

8.51  -  01 

Ba - 1 36n 

*.77  -  03 

Sb- 1 3 1 

7.07  -  01 

Bd-  1  37;,, 

7.68  -  07 

To- 1 7  l;r 

5.06  -  11 

Ra-139 

3.71  c  03 
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Table  2.  Radionuclide  ground  concentrations  (cont'd) 


Radionucl  i  do 

Acti  vi  ty 
( C  i  /  Km  - ' ) 

Radionucl ide 

Acti  vi ty 
(Ci/Km-') 

Radionucl ide 

Acti vi ty 
( C i / Km‘ ) 

Ba-140 

2.36 

4 

01 

Pm- 149 

5.20  -  04 

Tb-162m 

6.18  -  04 

Ba- 141 

2.50 

4 

03 

Pm- 150 

7.95  -  02 

Tb- 1 63 

1.36  -  03 

Ba-142 

7.53 

4 

0? 

Pm- 151 

5.83  -  02 

Tb- 1 64 

8.93  -  08 

La- 138 

2.47 

- 

17 

Pm- 152 

2.43  -  03 

Dy-165 

5.76  -  05 

La-140 

6.32 

- 

01 

Pm- 1  S2in 

1.32  -  01 

Dy-165m 

3.61  -  17 

La-141 

1.49 

+ 

03 

Pm- 153 

1.19  -  01 

Dy-166 

1.15  -  05 

La-142 

3.01 

+ 

03 

Pm- 1 54 

9 . 84  -  05 

Ho-166 

1.40  -  07 

La-143 

1.54 

+ 

03 

Pm- 1 S4m 

3.99  -  08 

Ho- 166m 

3.66  -  13 

La- 144 

6. 02 

- 

22 

Pm- 157 

5.39  -  14 

Ce- 141 

9.33 

- 

01 

Sm- 151 

3.26  -  09 

Co- 142 

2.05 

- 

12 

Sm-153 

9.61  -  03 

Ce- 143 

1.95 

+ 

0? 

Sin- 155 

2.33  -  09 

Ce-144 

9.42 

- 

01 

Sm-156 

7.02  -  01 

Ce- 145 

3.29 

- 

01 

Sm- 157 

3.61  -  01 

Ce- 1 46 

7.98 

+ 

02 

Sm- 1 58 

3.61  -  00 

Ce-147 

5. IS 

- 

11 

Sm- 1 59 

1.17  -  05 

Ce- 148 

9.82 

- 

21 

Sm- 160 

9.65  -  03 

Pr-142 

2.08 

- 

06 

Eu-152 

2.96  -  11 

Pr-142r 

5.04 

- 

06 

Eu- 1 52m 

3.36  -  07 

Pr-143 

2.98 

- 

01 

Eu-154 

2.87  -  08 

Pr-144 

9.33 

- 

01 

Eu- 1 55 

7.18  -  07 

Pr- 144r 

1.20 

- 

02 

Eu- 1 56 

2.59  -  04 

Pr-  145 

6.86 

+ 

02 

Eu- 1 57 

2.03  -  02 

Pr- 146 

2.66 

4 

03 

Eu- 1 58 

3.78  -  01 

Pr- 1 4  7 

5.01 

4 

02 

Eu- 1 59 

2.73  -  01 

Pr-143 

8.61 

- 

05 

Eu- 160 

1.39-20 

Pr-144 

4.80 

- 

04 

Eu-162 

2.18  -  04 

Nd- 147 

1.05 

f 

01 

Gd-153 

1.63  -  12 

Nd- 144 

5.22 

4 

02 

Gd-159 

2.16  -  03 

Nd- 151 

8.98 

4 

01 

Gd- 161 

1.67  -  05 

Nd- 1 5? 

4.  99 

4 

01 

Gd- 162 

1,24  -  02 

Nd- 1  S3 

1.08 

- 

12 

Gd- 163 

3.71  -  12 

Nd- 1 54 

4.29 

- 

01 

Gd- 164 

8.72  -  03 

Nd-  1 56 

1.14 

- 

16 

Gd-165 

3.73  -  12 

Pm- 1 4  7 

9.84 

- 

08 

Tb-160 

1.07  -  08 

Pm- 148 

5.34 

- 

06 

Tb-161 

4.59  -  06 

Pm- 1 48m 

6.97 

- 

07 

Tb- 162 

5.76  -  05 
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exception  of  the  isotopes  of  the  noble  gas  krypton  and  xenon,  that  are 
present  after  a  decay  time  of  one  hour. 


It  is  not  necessary  to  consider  all  the  radionuclides  listed  in 
Table  2  to  adequately  characterize  radiological  water  contamination. 
The  most  significant  radionuclides  have  been  identified  by  a  screening 
process  that  considered  the  relative  contribution  of  each  radionuclide 
to  potential  ingestion  doses.*  Based  on  unfractionated,  U-235 
time-dependent  fission  product  inventories  obtained  from  F 1 1  DOS  and  a 
50-year,  ingestion  dose  conversion  factors  obtained  from 
0RNL/NUREG/TM-190,^  the  percent  of  the  total  dose  commitment  to 
specific  body  organs  has  been  calculated  for  each  radionuclide.  The 
results  of  such  calculations  are  shown  in  Table  3  where  those 
radionuclides  that  contribute  more  than  1%  to  the  organ  doses  are 
identified.  Clearly,  a  host  of  radionuclides  could  be  significant  to 
the  ingestion  dose  from  contaminated  water  and  the  relative  importance 
of  specific  radionuclides  is  both  time  and  organ  dependent. 


For  the  radionuclide  screening  process  described  above,  the  focus 

was  on  the  internal  dose  associated  with  the  ingestion  of 

radiological ly  contaminated  water.  Similar  calculations,  using  dose 

(91 

rate  conversion  factors  for  external  exposure,  have  been  used  to 
determine  the  significance  of  radionuclides  to  the  external  dose  rate 
associated  with  immersion  in  radiological  ly  contaminated  water.  The 
external  exposure  pathway  is  not  of  concern  to  this  assessment; 
however,  the  external  dose  rate  does  provide  an  indication  of  the 
radiation  dose  rate  that  might  be  measured  by  radiation  instruments. 
Table  4  shows  the  results  of  the  external  exposure  calculations  for 
both  beta  and  gamma  radiation. 

*It  Ti  recognized  that  fission  product  fractionation  during  fallout 
formation  and  the  effects  of  solubility  during  the  contamination 
process  will  affect  the  radionuclide  composition  of  the  fission 
product  debris  in  water.  It  is  also  noted  that  the  use  of  50-year, 
dose  conversion  factors  might  not  be  appropriate  for  addressing 
wartime  doses.  Nevertheless,  this  approach  is  still  adequate  as  a 
screening  process  for  the  identification  of  important  radionuclides. 
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Table  3.  Percent  of  organ  dose  commitment  for  specific  radionuclides*. 


r  ; :  rt 


■no 


Ti:w 

1  2*2  [j*i  0*’. 

I.'  IS  1  .  2  i 

o  ■  1!  )? 

1 


W  tM  v.l  Y"  IQ 

29  4S  IS  ) 

17  .17  7 "  92  94 


♦Organ  dose  coirani tments  were  calculated  for  a  mixture  of  unfractionated 
U-235  fission  products  based  on  50-year  ingestion  dose  conversion  factors. 
Calculations  were  performed  for  specific  decay  times:  H-hours,  D-days. 
M-munths,  Y-years.  Totals  might  not  give  100  due  to  rounding. 
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Table  3.  Percent  of  organ  dose  commitment  for  specific  radionuclides  (cont) 
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Table  4.  Percent  of  external  dose  rate  for  specific  radionuclides  (cont) 


2.4.2  Radioactive  Decay 


The  process  of  radioactive  decay  is  important  to  considering  the 
time  dependency  of  the  radiological  water  contamination.  Since 
radionuclides  are  typically  part  of  a  chain  of  isotopes  connected  by 
the  process  of  radioactive  decay,  the  decay  of  one  radionuclide  is 
frequently  the  source  of  another  radionuclide. 

In  the  water  contamination  computer  model  used  for  this 
assessment,  radionuclides  were  treated  as  parent  or  daughter  members 
of  a  two-member,  radionuclide  chain.  Since  some  radionuclide  decay 
chains  contain  more  than  two  members,  this  approach  sacrifices 
accuracy  for  modeling  efficiency  and  simplicity.  However,  frequently 
only  two  members  of  the  chain  are  radionuclides  that  are  actually 
important  to  water  contamination. 

2.4.3  Dissolution  of  Radionuclides 


The  dissolution  of  radionuclides  in  fallout  from  the  solid  phase 
to  the  liquid  phase  is  probably  the  most  important  water  contamination 
process.  The  degree  to  which  the  radionucl ides  in  fallout  material 
are  soluble  or  insoluble  in  water,  and  the  rate  at  which  the 
radionuclides  dissolve  are  basic  to  the  modeling  of  radiological  water 
contamination.  Unfortunately,  the  available  information  on  the 
solubility  and  dissolution  rates  cf  radionuclides  in  fallout  is  quite 
1 imited. 

Most  of  the  statements  found  in  the  literature  regarding  fallout 
solubility  refer  to  fallout  as  basically  insoluble.  The  statements 
are  typically  based  on  measurements  of  fallout  solubility  in  terms  of 
the  gross  beta  radiation  activity  before  and  after  water  washing  of 
fallout  particles  collected  during  nuclear  weapons  tests.  Generally, 

IT  to  3T  of  the  gross  beta  activity  is  removed  from  the  fallout 

*  .  ,  .  .  ..  (10,11,12) 

material  by  such  washings. 
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There  are  very  few  references  in  the  literature  regarding  the 

specific  radionuclides  in  fallout  that  are  soluble.  A  group  of  six 

radionuclides  (i.e.,  Sr-89,  Sr-90,  Ru-106,  1-131,  Cs-137,  and  Ba-140) 

have  been  labeled  as  soluble  in  some  water  contamination 

studies. However,  no  substantial  discussion  of  the  solubility 

of  specific  radionuclides  in  fallout  has  been  found.  Similarly,  there 

are  few  references  in  the  literature  regarding  the  dissolution  rates 

of  radionuclides  in  fallout.  Solubility  rates  for  actual  fallout 

particles  are  discussed  by  Larson^1"^  and  some  experiments  with 
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artificial  fallout  particles  have  been  performed.  However,  no 

substantial  models,  theories,  or  results  have  been  found. 

In  the  absence  of  the  necessary  fallout  solubility  data,  the  use 
of  radionuclide-specific  distribution  coefficients  provides  a 
convenient  and  reasonable  approach  for  fallout  solubility  modeling. 
The  concept  of  the  distribution  coefficient,  Kd,  was  originally 
developed  from  ion-exchange  theory  to  represent  the  equilibrium 
distribution  of  a  trace  constituent  between  the  solid  exchanger  and 
the  solution.  Currently,  the  distribution  coefficient  concept  is 
being  used  to  quantify  the  chemical  interaction  of  radionuclides  with 
soils  and  minerals  for  the  assessment  of  the  contamination  of  water 
bodies  by  releases  from  nuclear  power  plants  and  for  the  analysis  of 
the  impacts  of  potential  releases  of  radioactive  material  from  nuclear 
waste  repository  facilities. 

A  distribution  coefficient  is  defined  as: 

Kd  =  amount  of  radionuclide  sorbed  on  solid  phase 
amount  of  radionuclide  left  in  solution 

Since  the  solid  phase  activity  is  usually  expressed  in  units  of  Ci/g 
and  the  liquid  phase  activity  in  units  of  Ci/m?  ,  Kd  typically  has 
units  of  ms  /g.  For  a  specific  element,  the  value  of  Kd  is  dependent 
upon  the  chemical  state  of  the  element,  the  type  of  solid  matrix  in 
which  it  exists,  the  physical  characteristics  of  the  solid  and  liquid 
phases,  and  the  nature  of  the  dissolution  process;  however,  the  actual 


relationship  of  the  value  of  Kd  to  these  parameters  is  generally  not 
known.  Values  of  Kd  are  normally  determined  by  laboratory  or  field 
experiments  and  can  exhibit  a  wide  range;  for  example,  the  Kd  value 
for  Zr  ranges  from  1000  to  10000  m  ■  / g  and  the  value  for  Sr  ranges 
from  8  to  4000  m  ?/g. 

Table  5  shows  selected  values  of  distribution  coefficients  for 
those  elements  whose  radioisotopes  are  considered  in  the  water 
contamination  model.  As  cautioned  in  Reference  16,  it  is  a  gross 
generalization  to  prescribe  single-valued,  non-specific  distribution 
coefficients;  although  it  is  recognized  that  such  parameters  are 
needed  by  computer  modelers  for  preliminary  or  scoping  studies. 
Reference  16  suggests  that  median  Kd  values,  like  those  shown  in  Table 
5,  should  be  considered  to  vary  by  a  factor  of  10  for  those  values 
greater  than  100  m  c /g. 

The  distribution  coefficient  refers  to  the  phase  distribution  of 
a  radionuclide  at  equilibrium.  For  the  fallout  material  deposited  on 
the  water  surface,  it  is  assumed  that  equilibrium  is  reached  within  a 
time  period  corresponding  to  the  time  constant  of  the  mixing  tank 
model  of  the  watershed  (see  Section  2.4.5).  For  the  fallout  material 
deposited  on  the  land  surface,  it  is  assumed  that  equilibrium  is 
reached  i nstantaneously  and  maintained  for  the  four  days  that  the 
runoff  water  remains  on  the  watershed.  These  assumptions  about  the 
time  to  achieve  equilibrium  are,  in  effect,  assumptions  about  the  rate 
at  which  the  radionuclides  in  the  fallout  dissolve.  As  pointed  out 
above,  there  is  no  strong  basis  for  supporting  such  assumptions; 
nevertheless,  such  assumptions  are  necessary  in  order  to  model  the 
water  contamination  processes. 

2.4.4  Transport  of  Radionuclides  by  Precipitation  Runoff 


The  importance  of  precipitation  runoff  with  respect  to  water 
contamination  was  not  generally  recognized  until  the  late  1960's. 
Prior  to  that  time,  the  focus  had  been  on  the  pollution  caused  by 
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Table  5.  Selected  di stribution  coefficients 


Kd 


Element 

(m :J31 

Ba 

500 

Ce 

10,000 

Cs 

1  ,000 

I 

10 

La 

500 

Mo 

25 

Nb 

10,000 

Nd 

10,000 

Pm 

10,000 

Pr 

10,000 

Rh 

5,000 

Ru 

5,000 

Sb 

100 

Sr 

1 ,000 

Tc 

1 

Te 

100 

Y 

1  , 000 

Zr 

1  ,000 

'  _  amount  of  radionuclide  sorbed  on  solid  phase  \,n/g) 
amount  of  radionuclide  left  in  solution  (Ci/m.) 

The  above  Kd  values  were  taken  from  References  16  and  18  and  represent 
mean  or  best-estimate  values.  Kd  values  greater  than  100  inr/g  should 
be  considered  to  vary  (+)  by  a  factor  of  10. 
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point  sources  as  opposed  to  nonpoint  sources  Now,  it  is  estimated 

that  nonpoint  sources  of  pollution  account  for  more  than  half  of  the 

(19) 

total  water  quality  problems  in  the  United  States. 

Radiological  water  contamination  by  nonpoint  sources  has  received 

little,  if  any,  attention  in  the  current  literature  on  the 

environmental  impacts  of  nuclear  facilities,  since  such  facilities 

correspond  to  point  pollution  sources.  The  principal  source  of 

radiological  nonpoint  pollution  has  been  the  atmospheric  testing  of 

nuclear  weapons;  however,  with  the  cessation  of  such  testing  by  the 

major  nations,  the  subject  has  received  little  attention. 

Nevertheless,  the  environmental  contamination  of  water  supplies  by 

fallout  has  been  noted  as  a  source  of  background  radiation^^  and  a 

possible  concern  in  connection  with  the  long-term  impact  of  strategic 
(21 ) 

nuclear  warfare. 

The  transport  of  radionuclides  from  the  land  surface  into  the 
local  river  or  stream  by  precipitation  runoff  is  an  important  water 
contamination  process  in  this  assessment.  As  pointed  out  in  Section 
2.3,  and  discussed  further  in  Appendix  A,  within  the  selected  nuclear 
warfare  scenario  area  sufficient  rainfall  occurs  with  such  frequency 
that  precipitation  runoff  is  a  common  occurrence.  Furthermore,  for  a 
typical  watershed  area,  the  stream  or  river  area  is  about  0.1T  of  the 
total  watershed  area,  indicating  that  the  land  surface  of  the 
watershed  area  provides  a  large  area  source  for  potential  water 
contamination. 

The  volume  of  surface  runoff  water  that  results  from  a  given  rain 
is  determined  by  the  area  of  the  watershed,  the  amount  of  water 
deposited  by  the  rain,  and  the  prior  precipitation  history  of  the 
watershed.  This  volume  of  water  flows  off  the  land  surface  into  the 
local  water  channel  is  a  non-linear  fashion  that  continues  for  about 
four  days.  The  subject  of  precipitation  runoff  is  discussed  in 
various  civil  engineering  texts,  the  specific  approach  used  for  this 
assessment  is  discussed  in  detail  in  Appendix  A. 
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While  the  potential  precipitation  runoff  water  remains  on  the 
land  surface,  it  effects  a  phase  separation  of  the  radionuclides 
present  in  the  deposited  fallout  material.  The  radionuclides  that  are 
partitioned  into  the  liquid  phase  are  subsequently  transported  off  the 
watershed  into  the  local  river  or  stream  by  the  runoff  water  flow, 
assuming  uniform  mixing  with  the  water.  The  radionuclides  that  are  in 
the  solid  phase  are  assumed  to  remain  on  the  land  surface  and  be 
subjected  to  further  dissolution  by  subsequent  rains. 

It  should  be  noted  that  this  approach  omits  two  processes  that 
would  actually  affect  the  water  contamination:  (1)  the  transport  of 
radionuclides  down  into  the  soil  depths  by  precipitation  that 
infiltrates  the  soil,  and  (2)  the  transport  of  radionuclides  as 
particulates,  or  solid  phase  material,  off  the  watershed.  The 
omission  of  radionuclide  leaching  into  the  soil  will  cause  the  model 
to  overpredict  the  water  contamination;  however,  this  effect  is 
expected  to  be  rather  minor  and  only  of  concern  at  late  times,  e.g., 
several  weeks  after  the  initial  fallout  contamination.  The  omission  of 
particulate  transport  is  not  likely  to  be  of  much  significance  and  is 
acceptable  as  long  as  it  is  understood  that  the  model  is  mainly 
concerned  with  the  radionuclides  dissolved  in  the  water. 

2.4.5  Aquatic  Mixing  and  Transport 


A  host  of  aquatic  mixing  and  transport  models  are  described  in 

/  1  c  oo  \ 

the  technical  literature.'  ’  Most  of  these  models  are  designed 
for  situations  in  which  the  water  contaminant  is  introduced  by  a 
continuous,  point  source.  Sophisticated  models  that  can  address 
various  mixing  and  transport  parameters  and  provide  time-dependent 
modeling  are  also  available.  In  general,  however,  the  available 
models  appeared  to  be  too  detailed  and  too  complex  to  warrant 
application  to  this  assessment;  in  addition,  it  was  not  obvious  that 
modifications  to  such  models  to  accomodate  a  nonpoint  contamination 
source  would  be  warranted. 
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For  the  purposes  of  this  assessment,  the  stream  or  river  within 
the  watershed  is  assumed  to  be  a  "mixing  tank"  with  instantaneous  and 
uniform  mixing  of  the  radioactive  contaminant  with  the  water.  The 
concentration  of  the  radioactive  contaminant  at  the  outlet  of  the 
mixing  tank  (which  corresponds  to  the  location  of  the  water  supply 
point)  is  then  determined  by  the  intitial  concentration  in  the  tank 
and  the  number  of  volume  changes  per  unit  time  (the  flow  rate  divided 
by  the  tank  volume).  The  initial  concentration  of  the  radioactive 
contaminant  is  obtained  from  fallout  calculations  that  give  the  areal 
density  of  deposited  radionuclides.  Information  on  the  water  volume 
and  flow  rate  for  various  watershed  is  obtained  from  the  water  supply 
point  study  given  in  Appendix  A. 

2 . 5  Water  Contami n at ion  Model 

A  simplified  computer  model  has  been  developed  to  characterize 
radiological  water  contamination.*  The  Watershed  Water  Contamination 
Model  (WSWCM)  calculates  the  time-dependent  activity  concentration  of 
fission  product  radionuclides  dissolved  in  water  that  could  result 
from  the  deposition  of  nuclear  weapons'  fallout  on  a  watershed.  WSWCM 
considers  both  the  prompt  water  contamination  that  would  result  from 
the  fallout  material  deposited  directly  in  the  water  and  the  delayed 
water  contamination  that  would  result  from  the  fallout  material 
initially  deposited  on  the  land  surface  and  subsequently  transported 
to  the  water  by  pr  ec i  p i  ta  ti  on  runoff.  All  activity  is  assumed 
initially  to  be  associated  with  solid  particulate  fallout.  The 
activity  may  leave  the  watershed  only  by  radioactive  decay  or  by  being 
dissolved  in  water  which  flow  past  the  water  supply  point. 


*A  detailed  description  of  the  model,  WSWCM,  is  provided  in  Appendix  B. 
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The  principal  characteristics  of  WSWCM  are:  (1)  the  watershed  is 
modeled  as  a  mixing  tank,  (2)  radionuclide-specific  distribution 
coefficients  are  used  to  address  fallout  solubility,  and  (3)  the  model 
treats  radioactive  decay  including  daughter  in-growth.  It  is 

important  to  note  that  WSWCM  addresses  radioactive  material  in 

solution  but  does  not  incorporate  any  modeling  of  particulate  or 
sediment  transport. 

The  input  data/information  used  by  WSWCM  falls  into  three  main 
areas:  (1)  fission  product  radionuclide  data,  (2)  information  on 

watershed  and  precipitation  characteristics,  and  (3)  data  and 
information  on  the  solubility  modeling.  The  data  necessary  to  perform 
water  contamination  calcu1 ations  for  the  selected  scenario  area  is 
included  in  WSWCM.  The  principal  product  of  WSWCM  are  plots  of  the 
activity  concentration  in  water  ( pC i /  )  as  a  function  of  the  time 
since  the  fallout  material  was  deposited  (hr),  for  specified 
parent-daughter  radionuclide  pairs. 

WSWCM  was  applied  to  a  situation  in  which  the  initial  fallout 
contamination  of  a  specific  watershed  was  1  R/Hr  at  H+i  hour  and  a 
rain  occurred  on  the  6th  day  after  the  fallout  was  deposited.  The 
results  of  the  model  are  shown  in  the  following  figures.  As  can  be 

seen  in  the  figures,  three  typical  types  of  radionuclide  behavior  are 

observed.  Figure  3  shows  an  example  of  a  radionuclide  (1-135)  that 
decays  so  rapidly  that  it  presents  an  initial,  but  not  a  delayed, 
water  contamination  effect.  Figure  4  shows  an  example  of  a 
radionuclide  parent-daughter  pair  (Sr-90,  Y-90)  that  presents  a  rather 
constant  level  of  water  contamination.  Figure  5  shows  an  example  of  a 
radionuclide  parent-daughter  pair  (Te-131m,  1-131)  that  presents  both 
an  initial  and  a  delayed  water  contamination  effect.  Finally,  Figure  6 
shows  the  composite  effect  of  all  the  radionuclides  considered  for  the 
specific  problem.  Clearly,  both  the  initial  and  the  delayed  water 
contamination  processes  are  of  possible  significance.  In  addition, 
radionuclides  having  long  half-lives  and  low  solubilities  could  cause 
a  rather  persistent  water  contamination  problem. 
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Problem  1 


Watershed  No.  1.  J 

Fallout  =  1  R/Hr  at  H+l  hoi}r. 
Single  rain  at  H+144  hour. 


C i  »!i  t ai,;i  na  oi  on  CaujoJ  hy 
R.jnof f  of  fallout  Ik'posi ted 
i‘n  1  and 


Fiqure  3.  1-135  water  contamination  -  problem  1. 


Problem  1 . 

Watershed  No.  1.  : 

Fallout  =  1  R/Hr  at  H+l  hour 
Single  rain  at  H+144  hour. 


Figure  4.  Sr-90,  Y - 90  water  contamination  -  problem  1. 


Problem  1. 


Watershed  No.  1. 

Fallout  ^  1  R/.Mr  at  FM  hour 
Single  rain  at  H  + 1 4 4  hour. 


Figure  (.  Total  '10  radionuclide'.]  water  contamination 


problem  1 


To  better  illustrate  the  above  concerns,  WSWCS  was  applied  to 
another  watershed  that  was  contaminated  at  the  level  of  1  R/Hr-  at  H+l 
hour  but  with  21  rains  occurring  over  a  period  of  about  40  days.  The 
results  of  this  simulation  are  shown  in  Figures  7  through  10  for  the 
same  radionuclides  considered  in  the  previous  problem.  The  presence  of 
persistant  water  contamination,  albeit  at  a  level  much  below  the 
initial  water  contamination  level,  is  clearly  illustrated. 

As  seen  in  the  figures,  WSWCM  provides  a  convenient  tool  for 
analyzing  water  contamination  problems.  However,  it  should  be 
recognized  that  the  model  is  really  quite  simple  and  contains  several 
rather  basic  assumptions  about  the  behavior  of  fallout  particle 
radionuclides  in  water.  WSWCM  should  be  considered  as  a  scoping 
method  that  provides  order-of-magni tude  results. 
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Problem  2.  i 

Watershed  No.  8.'  1 

Fallout  =  1  R/Hr  at  H+l  hour 
21  rains  over  40  days. 

Note  time  scale  -  hour  x  10] 
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1-135  water  contamination  -  problem  2. 


Y-10  water  contamination  -  nroblen  2 
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Problem  [2.  ; 

Watershed  No.1  8. 

Fallout;  =  1  F^/Hr  at  H+l  hour,  | 
21  ra  injs  over  40  d4ys. 

Note  t i hie  scale  -  hour  x  10.  j 
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!es)  water  contamination  -  problem  2. 
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THRLAT  RECOGNITION 

3 . 1  Introduction 

Threat  recognition  refers  to  the  detection  and  measurement  of 
radiological  water  contami nat ion.  Technical  performance  criteria  for 
radiation  detection  and  measuring  equipment  are  based  on  the 
requirement  to  show  conformance  with  specified  water  quality 
standards.  F i el d- portabl e  radiation  detection  and  monitoring  equipment 
currently  available  and  under  development  does  not  appear  to  have  the 
capability  to  adequately  monitor  water  contamination  at  low  levels  of 
radioactivity.  Procedures  for  radiation  detection  and  measurement 
could  be  used  to  indicate  compliance  with  water  quality  requirements, 
in  lieu  of  demonstrating  compliance. 

3 . 2  bate*"  Quality  Standards 

esently,  standards  and  requirements  that  relate  to 
rad  i  o  I  og  i  ca  1  1  y  contaminated  water  are  contained  in  a  1J.  S.  Army 
technical  bulletin  and  a  NATO  standard i zat i on  agreement.  As  discussed 
below,  both  of  the  documents  are  currently  being  revised  and  their 
final  forms  cannot  yet  be  ascertained. 

The  most  specific  water  quality  standards  regarding  water 
contaminated  with  radioactive  material  are  given  in  TB  MFD  229.  “’Vb 
MTD  220  addresses  both  fixed  installation  and  field  water  supplies; 
for  fi-  Id  water  supplies,  radiological  water  contamination  criteria 
are  given  for  short-term  and  long-term  usage  (see  Figure  11).  The 
short-term  criteria  are  not  specific  with  regard  to  acceptable  levels 
of  water  contami nat i on  and  no  numerical  criteria  are  given.  The 
long-term  criteria  do  provide  specific  numerical  Ciiteria  (i.e.,  1000 
pC i/  gross  beta  activity  and  10  pC i /.  stront. i um-90)  by  imposing  the 
criteria  for  fixed  installations  on  the  field  water  supplies. 

4  b 


TB  MED  229  Requirements 
Radiological  Contamination  -  Field  Supplies 


Short  Term 
(less  than  7  days) 


Long  Term 
(more  than  7  days) 


For  short  term  consumption,  no  absolute 
numerical  standard  is  recommended  or 
considered  necessary  This  is  based  on 
the  conclusion  that  if  the  external 
radiation  hazard  permits  occupancy 
of  the  water  point,  tne  water  is 
suitable  for  consumption  during 
occupancy  not  exceeding  the  one-week 
period 

(Same  as  for  fixed  installations) 

Gross  Beta  Activity  1000  pC/ 

Stronti um  90  10  pC / . 


Figure  11.  TB  MED  22')  requirements 


i 
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TB  MED  229  is  currently  being  revised  and  this  revisioning  will 
result  in  two  new  technical  bulletins.  FB  MED  576  will  address  the 
fixed  installation  portion  of  TB  MED  229;  TB  MED  577  will  address  the 

field  installations  portion  of  TB  MED  229.  For  fixed  installations,  TB 

MED  576  will  adopt  new  drinking  water  standards  that  impose  stringent 
limits  on  the  allowable  rad i oact i v i ty  (i.e.,  50  pCi/;  gross  beta 
activity  and  8  pCi/e  strontium-90). 

At  present,  it  is  not  known  what  regu ; rements  will  be  imposed  on 
field  instal lations  by  TB  MED  577.  Since  the  field  installations 
portion  of  TB  MED  229  will  remain  in  effect  until  TB  MED  577  is 
published,  the  new  water  contami nat i on  limits  for  fixed  installations, 
to  be  given  in  TB  MED  5/6,  will  not  be  automat  ical  ly  adopted  for  the 

field  i  nst  al  lat  ions.  It  is  possible  that  TB  MED  57  7  will  not  specify 

any  numerical  criteria  for  radiological  contaminated  water. 

Radiological  water  contamination  is  also  addressed  in  the  NATO 
standardization  agreement  S’ANAG  2136  (MED)  -  Minimum  Standards  -of 
Water  Potability.  Figure  12  gives  the  present  standard  and  proposed 
N A'O  and  US  revisions  to  the  standard.  Specific  numerical  criteria 
for  radi o 1 ugi ca I ly  contaminated  water  are  not  provided  by  the  current 
standard  nor  the  proposed  revisions.  It  is  important  to  note  that  the 
absence  of  specific  numerical  criteria  does  not  imply  that  the 
consumption  ot  radi ol og i ca 1 ly  contaminated  water  is  acceptable.  On 
the  contrary,  the  philosophy  behind  the  non-numer i cal  criteria 
statements  is  that  any  radiation  ex'osure,  beyond  that  attributable  to 
normal  background  radiation,  should  be  avoided,  if  at  all  possible. 

At  present,  the  requirements  of  TB  MED  229  are  in  effect  and  thus 
provide  technical  performance  criteria  for  radiation  detection  and 
measuring  equipment.  Specifically,  such  equipment  must  have  the 
capability  to  demonstrate  that  water  intended  for  consumption  does  not 
curtain  more  thar  1000  pC i / .  of  gross  beta  activity  and  more  than  10 
pU  :  /  .  of  stront i um-90. 
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STANAG  2136  (MED)  -  Requirements 
***  Present  Standard  *** 
d .  Radiologi cal  Standards 

For  short  term  consumption  (1  to  7  days)  no  absolute  maximum 
tolerance  is  recommended  or  considered  necessary.  This  is 
based  on  the  consideration  that  if  the  risk  of  external 
radiation  is  such  as  to  allow  the  source  to  be  used,  then  the 
water  will  be  suitable  for  drinking  during  occupancy  not 
exceeding  one  week 

***  NATO  Proposed  Standard  *** 
d  Radiological  Standards 

It  is  undesirable  to  drink  water  contaminated  with  radioactive 
substances.  Consideration  should  be  given  as  to  whether  a 
source  of  water  is  likely  to  be  contaminated  In  some  circum¬ 
stances,  sources,  such  as  underground  water,  may  be  safely 
assumed  to  be  uncontaminated  Filtered  water  will  be  free  of 
insoluble  particlates  but  may  still  contain  soluble  radioisotopes 
such  as  iodine.  If  there  is  any  doubt,  monitoring  of  the  water 
should  be  attempted  This  can  be  achieved  simply  either  by  the 
use  of  a  dose  rate  monitor  placed  close  to  a  large  sample,  e.g., 
bucketful  1  or  preferably  by  drying  down  a  sample  and  using  a 
contamination  meter  Any  water  sample  showing  a  reading  above 
background  should  only  be  used  if  no  better  source  is  available 
and  the  use  is  essential. 

***  US  Proposed  Standard  *** 

d  Radiological  Standards 

A  standard  in  the  normal  sense  of  definite  limits  is  not  appropriate 
in  this  case,  however,  the  following  procedures  should  be  employed: 

1.  Areas  Having  Received  Fallout  -  For  short  term  consumption  (1  to 

7  daysT  no  absolute  maximum  tolerance  is  recommended  or  considered 
necessary  This  is  based  on  the  consideration  that  if  the  risk 
of  external  radiation  is  such  as  to  allow  the  source  to  be  used, 
then  the  water  will  be  suitable  for  drinking  during  occupancy  not 
exceeding  one  week. 

2.  Areas  Not  Having  Received  Fallout  -  For  short  term  consumption 
(1  to  7  days )  any  water  sample  showing  a  reading  above  back¬ 
ground,  as  measured  with  a  dose  rate  meter  or  other  suitable 
method,  should  only  be  used  if  no  better  source  is  available 
and  the  use  is  essential  This  is  based  on  the  consideration 
that  personnel  should  not  be  subjected  to  unnecessary  radiation 
exposure. 


Figure  12.  STANAG  2136  (MED)  -  requirements 


3. 3  Radiation  Monitoring  Equipment 


Field-portable  radiation  detection  and  measuring  equipment  for 
monitoring  radiological  water  contamination  has  been  developed  for 
both  military  and  civilian  applications.  As  discussed  below,  neither 
the  military  equipment  nor  the  civilian  equipment  appears  to  be 
capable  of  monitoring  radiological  water  contamination  at  the  levels 
of  1000  p C i / v  (gross  beta  activity)  and  10  pC i /  •  (strontium-90). 

The  Radiac  Set  AN/PDR-27  is  a  low-range,  beta-gamma  instrument, 
standard  for  all  Services  and  used  for  personnel  and  equipment 
monitoring.  It  has  two  Ge i ger-Mue 1 1 er  tubes  --  a  large  tube  for 
low-range  detection  and  a  small  tube  for  high-range  detection.  A  beta 
shield  on  the  large  G-M  tube  permits  the  measurement  of  beta  activity 
by  using  the  difference  between  unshielded  (B  +'■■')  and  shielded  (t 
only)  measurements.  Ihe  AN/PDR-27  covers  a  range  of  0  to  500  mR/hr  in 
four  decade  steps;  beta  measurements  by  the  difference  technique  are 
only  possible  on  the  two  lowest  ranges  (0-0.5  mR/hr  and  0-5  mR/hr). 

The  AN/PDR-27  can  be  used  to  give  a  qualitative  indication  of 

beta-gamma  water  contamination  by  holding  the  large  G-M  tube,  without 

the  beta  shield,  about  one-half  inch  from  the  surface  of  the 
( 24 ) 

water. x  ’  A  reading  above  background  indicates  a  concentration  of  an 
unidentified  beta-gamma  emitter  in  excess  of  10^  pCi/.'.  ^ 

Another  water  monitoring  procedure  with  the  AN/PDR-27  has  been 
described  by  D.  C.  Lindsten.^)  jn  this  procedure,  the  large  G-M 
tube,  without  the  beta  shield,  is  encased  in  a  rubber  surgical  glove 
and  immersed  into  the  water.  The  minimum  level  of  detection  is 
estimated  to  be  5  x  10  pCi/v  of  mixed  fission  products  (beta-gamma 
emitters)  based  on  correlation  data  developed  using  fallout  material 
from  nuclear  weapons  testing. 

A  new  beta-gamma  survey  meter,  Radiac  Set  AN/VDR-1  is  being 
developed  as  a  replacement  for  the  AN/PDR-27.  When  used  for  surface 
monitoring  of  contaminated  water,  the  AN/VDR-1  is  expected  to  be  able 
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to  detect  a  beta  emitting  isotope  at  the  level  of  3.5  x  10&  pC i /  -  At 

present,  the  status  of  the  AN/VDR-1  is  uncertain,  and  further 

(27) 

development  or  production  might  be  discontinued. 

Rad i at  1  on  detection  and  measuring  equipment  developed  for 
civilian  applications  is  described  in  numerous  textbook  and  other 

pu  b  1  i  c  a  t  i  ons  .  )  in  order  to  supplement  such  references  arid 

determine  the  current  state-of -the-art ,  a  request-for- i  nfoririat  i on  was 
sent  to  43  firms  identified  as  providing  equipment  or  services  related 
to  the  detection  and  measurement  of  radiological  water  contami r.at. i on. 
The  identification  of  the  firms  was  based  on  their  cited  capabilities 

( 33 ) 

listed  in  the  Nuclear  News  1981  Buyers  Guide.'  '  The 

request-for-information  letter  expressed  an  interest  in  equipment  that 
could  {Measure  beta  radiation  in  a  mixed  beta-gamma  contami nated  water 
sample.  The  range  of  interest  was  given  as  from  3  x  10b  pCi/i  down  to 
3  pCi/  ,  and  if  was  noted  that  the  gamma  radiation  could  be  a  factor 
of  10  above,  or  below,  the  beta  radiation.  An  interest  in  both  gross 
beta  radiation  measurements  ana  specific  radionuclide  identification 
was  expressed.  It  was  noted  that  the  principal  interest  was  with 
equipment  that  could  be  used  for  real-time  monitoring  o+  a  process 
stream;  however,  there  was  also  interest  in  1 aboratory-type  equipment 
if  such  equipment  could  be  trans;  ortable  to  the  f  ield  and  provide  an 
analysis  within  a  couple  of  hours. 

Table  6  identifies  the  firms  that  were  contacted  and  indicates 
their  responses.  Of  the  43  firms  contacted,  24  did  not  respond  to  the 
request.  Of  the  19  responding  firms,  9  indicated  that  they  did  not 
have  the  capabilities  for  measuring  the  water  contamination  of 
interest;  10  firms  indicated  that  they  did  have  some  capability  and 
provided  relevant  literature. 

Based  on  an  assessment  of  the  supplier  literature  and  the 
references  cited  above,  it  appears  that  for  process  monitoring 
applications  the  minimum  level  of  detection  is  about  10°  pCi/c  gross 
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Table  6.  Responses  to  request- for- in  format  ion 


Firm 


No 

Response 


Negative 

Response 


Posi ti ve 
Response 


Anacon,  Inc/Aero  Vac  Products  x 

Applied  Health  Physics,  Inc  x 

Applied  Physical  Technology,  Inc  x 

Aptec  Nuclear  Inc  x 

The  Aston  Company  x 

Baird  Corp 

Berthol d-Beta  Analytical,  Inc  x 

Canberra  Industries,  Inc. 

Cedar  Grove  Operations 
Centronic  Inc 

Don  L  Col  1 i ns  &  Assoc.  x 

Dionex  Corp.  x 

Dosimeter  Corp  of  America 
Eberline  Instrument  Co.  x 

EG&G  Ortec  Inc. 

Electrometer  Corp  x 

Evans  Nuclear  Consulting  Services  x 

Foxboro  Analytical 

Gamma-Metrics  x 

General  Atomic  Co  x 

The  Harshaw  Chemical  Co. 

High  Voltage  Engr.  Corp  x 

IRT  Corp 

Kaman  Sciences  Corp. 

National  Nuclear  Corp 
Nuclear  Data  Inc 

Nuclear  Equipment  Chein  Corp  x 

Nuclear  Instrument  Co.  x 

Nuclear  Measurement  Corp.  x 

Nuclear  Research  Corp.  x 

Princeton  Gamma-Tech  ,  Inc 
Reuter-Stokes ,  Inc. 

J.  L  Shepard  &  Assoc 

Technical  Assoc  x 

Technology  for  Energy  Corp.  x 

Teledyne  Analytical  Inst 
Tennelec,  Inc 

Tera  Corp.  x 

United  States  Testing  Co  ,  Inc.  x 

Victroneen  Instrument,  Inc. 

Westinghouse  Electric  Corp.  x 

Weston  Components 

Xetex  Inc  x 


x 


x 

x 

x 


x 


x 


x 


x 

x 

x 

x 


x 

x 

x 


x 


x 


x 

x 


(43) 


(34)  (9)  (10) 
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beta  activity  and  specific  radionuclide  detection  is  not  possible.* 
It  is  possible,  however,  to  measure  gross  beta  activity  at  a  level 
of  1  pC i / •  and  determine  specific  radi onucl i des ,  like  strontium-90,  at 
levels  of  1  to  10  pCi/i  by  using  1 aboratory-type  equipment  and 
procedures.  However,  the  necessary  equipment  with  its  radiation 
shielding  and  supporting  utilities  is  not  field-portable  and  the 
procedures,  which  include  chemical  processing  of  the  sample  and 
counting  times  of  24  to  48  hours,  are  too  time  consuming  and 
complicated  to  be  used  by  troops  in  the  field  situation. 

3.4  Radiation  Monitoring  Procedures 

As  discussed  above,  radiation  detection  and  monitoring  equipment 
d'^es  not  have  the  capability  to  detect  radiological  water 
contamination  at  the  levels  required  to  meet  current  standards  (i.e., 
1 U 0 0  pC i / ■  gross  beta  activity  and  10  pC i /  •  strontium-90).  D.  C. 
Linsten  has  suggested  a  procedure,  termed  "supply-side  nuclear  water 
monitoring",  that  could  be  used  to  indicate  compliance  with  the  water 
quality  standards. ^ ^ 

lindsten's  procedure  is  based  on  monitoring  the  radiological ly 
contaminated  water  as  it  passes  through  the  stages  of  purification 
rather  than  monitoring  the  finished  product  water.  As  discussed  in 
Section  4,  the  water  purification  stages  include:  coagulation  and 

filtration,  reverse  osmosis,  and  ion  exchange.  The  coagulation  and 
filtration  stage  is  expected  to  remove  all  of  the  insoluble  fallout 
material  from  the  raw  water,  but  not  affect  the  soluble  material.  The 
soluble  material  is  expected  to  be  removed  with  efficiencies  of  99% 
and  99.9%  by  the  reverse  osmosis  and  ion  exchange  stages,  respectively. 


*It  should  be  noted  that  a  variety  of  radiation  detection  and  measuring 
equi  'lent  (and  data  analysis  capability)  is  available  for  gamma 
radiation  monitoring.  It  might  be  possible  to  use  this  technology 

for  water  monitoring  to  achieve  lower  minimum  levels  of  detection  or 
simplify  the  water  monitoring  problem. 
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Based  on  the  above  stated  removal  efficiencies,  water  entering 
the  ion  exchange  stage  with  a  gross  beta  activity  of  10^  pCi/c  will  be 
decontaminated  to  a  level  of  1000  pCi/?  .  Similarly,  water  entering 
the  reverse  osmosis  stage  with  a  gross  beta  activity  of  100  x  10^  pCi/- 
will  be  decontaminated  to  a  level  of  10^  pCi/;.  These  levels  of 
activity  for  the  water  prior  to  the  purification  stages  are 
sufficiently  high  to  be  measured  with  existing  radiation  detection  and 
measuring  enuipment. 

According  to  Lindsten,  water  that  meets  the  criteria  of  1000  pCi/? 
gross  beta  activity  will  meet  the  critera  of  10  pCi/-  strontium-90  for 
the  first  200  days  after  the  nuclear  explosion.  His  argument  is  based 
on  the  fact  that  for  gross  fission  Droducts  the  strontium-90  activity 
is  less  than  IT  of  the  total  fission  product  activity  for  decay  times 
less  than  0.7  years.  Actually,  the  comparison  should  be  made  between 
the  strontium-90  activity  and  the  total  activity  of  those  soluble 
radionucl ides  present  in  the  process  stream  between  the  purification 
stages.  ised  on  the  water  contamination  model  described  in  Section 
2,  is  ak  ears  that  the  strontium-90  activity  is  about  \l  of  the  total 
activity  of  the  radionucl  ides  in  solution,  thus  indicating  that  the 
strontium-90  criter  would  be  met  if  the  gross  activity  criteria  were 
met.* 

Lindsten's  "supply-side  nuclear  water  monitoring"  approach  is 
technically  feasible.  However,  as  noted  by  Lindsten,  to  implement  the 
procedure  the  decontamination  capabilities  of  the  water  purification 
processes  must  be  precisely  known.  In  particular,  it  appears  that  it 
would  be  necessary  to  measure  and  validate  the  removal  efficiencies  of 

*It  shoul d  5?  pointed  out  that  it  is  not  really  correct  to  compare 
these  activities  in  such  a  simple  fashion.  Actually,  one  should 
compare  the  activities  that  would  be  measured  by  the  radiation 
monitoring  equipment. 


each  process  on  field  equipment  while  actually  deployed.  This  does 
not  mean  that  calibration  activities  would  necesarily  be  required 
during  a  wartime  situation,  but  it  does  mean  that  the  necessary 
calibration  procedures  would  have  to  be  developed  and  applied 
routinely  to  monitor  the  status  of  the  water  purification  equipment  to 
assure  its  combat  readiness. 


SECTION  4 

THREAT  COUNTERMEASURES 

4. 1  Introduction 

Threat  countermeasures  include  water  purification  equipment  and 
methods,  and  field  operations  policy  and  procedures.  The  water 
purification  equipment  of  interest  is  the  Reverse  Osmosis  Water 
Purification  Unit.  The  field  operations  policy  and  procedures  include 
such  measures  as  water  point  selection,  water  treatment  scheduling, 
water  storage,  etc.  As  discussed  below,  this  assessment  has  been 
primarily  concerned  with  the  water  purification  equipment. 

4 . ?  Water  Purification  Equipment 

To  meet  the  requirements  for  a  multi-purpose  water  purification 
unit  to  provide  potable  water  in  the  field,  the  U.  S.  Army  is 
developing  the  Reverse  Osmosis  Water  Purification  Unit  (RWPU).  The 
ROWt'j  is  intended  to  purify  raw  water  contaminated  with  biological, 
chcical,  or  radiological  alenals.  The  ROWPU  will  be  highly  mobile 
and  available  in  f-OO  Cl'i!  (gallons  per  hour),  1500  GPH,  and  3000  GPH 
units. 

The  ROWPU  produces  potable  product  water  from  radiological ly 
contaminated  raw  water  by  a  series  of  three  purification  processes: 
(1)  coagulation  and  filtration,  (2)  reverse  osmosis,  and  (3)  ion 
exchange.  The  coagulation  and  filtration  process  is  ’ntended  to 
remove  all  of  the  insoluble  radioactive  material  from  the  raw  water; 
soluble  radioactive  material  will  pass  through  the  coagulation  and 
filtration  process  unaffected.  The  reverse  osmosis  p'rocess  is 

expected  to  remove  9^"  of  the  soluble  radioactive  material;  this 
process  will  also  back-up.  the  coagulation  and  filtration  process  by 
removing  any  soluble  radioactive  material  that  is  present.  The  ion 
exchange  process  is  expected  to  remove  99.9%  of  the  soluble 
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radioactive  material.  Overall,  the  ROWPU  is  expected  to  remove  any 
insoluble  radioactive  material  present  in  the  raw  water  and  to  reduce 
the  concentration  of  soluble  radioactive  material  in  the  raw  water  by 
a  factor  of  10^. 

5 

Based  on  a  decontamination  factor  of  10  and  a  product  water- 
criteria  of  1000  pc i / ?  gross  beta  activity,  the  ROWPU  can  effectively 
handle  raw  water  radiological ly  contaminated  up  to  the  level  of  10” 
pC i /  .  with  soluble  radioactive  material.  Based  on  the  water- 
contamination  model  discussed  in  Section  2,  a  watershed  contaminated 
by  fallout  at  a  level  of  1  R/Hr  at  H+l  hour  would  yield  water 

g 

contaminated  at  a  level  in  excess  of  10  p C i /  :  for  about  8  to  10 
hours.  This  peak  level  of  contamination  would  diminish  rapidly  and,  in 
the  absence  of  water  contamination  introduced  by  precipitation  runoff, 
would  be  less  than  10°  pCi/;  within  3  to  4  days.  Precipitation  runoff 
could  also  case  the  level  of  water  contamination  to  increase  to  10^  to 
107  pCi /  within  this  3  to  4  day  period;  precipitation  runoff  could 
also  cause  the  level  of  water  contamination  to  remain  at  10  4  to  10  5 
pC i / r  for  several  weeks  after  the  fallout  had  been  deposited. 

Assuming  that  the  problem  of  the  initial  (peak)  radiological 
water  contamination  can  be  avoided  by  water  storage  or  rationing,  the 
ROWPU  could  initially  handle  water  from  an  area  contaminated  by 
fallout  at  a  level  of  1  to  10  R/Hr  at  H+l  hour.  After  about  4  days, 
the  ROWPU  could  handle  water  from  areas  contaminated  by  fallout  at  a 
level  of  about  100  R/Hr  at  H+l  hour.  For  several  weeks  after  the 
fallout  deposition,  the  ROWPU  would  still  be  needed  to  purify  the  raw 
water  to  acceptable  potable  water  criteria. 

As  was  shown  in  Section  2,  surface  burst  nuclear  weapons  can 
produce  fallout  contamination  of  rather  large  areas.  However,  to 
estimate  the  extent  and  the  intensity  of  such  fallout  contamination  in 
the  event  of  nuclear  warfare  requires  major  assumptions  regarding  the 
nature  of  the  nuclear  strikes,  the  yields  of  the  weapons,  and  the 
prevailing  meteorological  conditions.  Statements  about  the  adequacy  of 
ROWPU  to  provide  potable  water  in  the  nuclear  warfare  environment  can 
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only  be  made  in  the  context  of  a  specific  nuclear  warfare  scenario. 
Since  this  assessment  does  not  address  such  scenarios,  no  specific 
comments  on  the  adequacy  of  ROWPU  are  offered.* 

However,  this  assessment  has  identified  two  design-related  areas 
that  merit  special  mention.  First,  the  ROWPU  should  include  within 
its  design  those  features  and  support  equipment  that  provide  a 
capability  for  in-field  testing  of  the  removal  efficiencies  of  the 
purification  processes;  this  capability  will  be  of  paramount 
importance  if  concepts  such  as  "supply-side  nuclear  water  monitoring" 
are  adopted.  Second,  the  ROWPU  should  have  an  availability 
characteristic,  achieved  through  low  equipment  failure  rates  and  short 
maintenance/repair  times,  adequate  to  ensure  that  the  equipment  can  be 
operated  nearly  continuously  for  several  weeks. 


Operations  Policy  and  Procedures 


Field  operations  policy  and  procedures  could  complement  the  water 
purification  equipment  to  provide  countermeasures  to  the 
potential  radiological  water  contamination  threat.  The  policy  and 
procedures  inc'ude  water  point  selection,  water  treatment  scheduling, 
water  storage,  etc. 


The  effectiveness  of  such  countermeasures  can  only  be  examined  in 
the  context  of  specific  nuclear  warfare  scenarios  and  force 


♦Earlier  in  this  assessment,  a  nuclear  warfare  scenario  including  force 
deployments  was  developed.  For  this  scenario,  radiological  ware  con¬ 
tamination  was  not  a  significant  problem  because:  only  a  few  surface 
bursts  occurred,  the  prevailing  wind  was  blowing  away  from  the  area 
where  the  U.  S.  forces  were  deployed,  and  the  wind  shear  was  so  small 
that  the  fallout  pattern  exhibited  little  width.  However,  during  the 
development  and  analysis  of  the  scenario  it  became  clear  that  a 
single,  unique,  hypothetical  scenario  does  not  provide  an  adequate 
basis  for  making  judgments  on  the  worth  or  utility  of  water  purifica¬ 
tion  systems  or  related  policy  and  procedures. 
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deployments.  Since  this  assessment  does  not  address  such  scenarios, 
no  specific  comments  on  the  utility  of  field  operations  policy  and 
procedures  are  offered. 


However,  this  assessment  has  identified  two  aspects  of  field 
operations  that  merit  special  mention.  First,  sufficient  product  water 
storage  capability  should  be  available  to  satisfy  requirements  if  the 
water  purification  units  have  to  curtail  operations  for  about  a  day 
because  of  an  inability  to  handle  the  intial  (peak)  radiological  water 
contamination.  Second,  strict  radiological  defense  measures  should  be 
maintained  in  effect  until  proper  radiological  water  contamination 
monitoring  has  determined  that  the  problem  no  longer  exists;  such 
procedures  are  important  because  the  water  could  remain  radiol ogical ly 
contaminated  even  after  area  radiation  monitoring  does  not  detect  any 
military  significant  fallout. 


c  o 

J  ( ) 


SECTION  5 


CONCLUSIONS 


In  the  event  of  nuclear  warfare  with  nuclear  weapons  employed  in 
a  surface  burst  mode,  the  fallout  contamination  of  watersheds  and 
water  supplies  would  be  sufficiently  high  to  require  the  use  of  water 
purification  equipment  to  produce  potable  water  that  meets  the  current 
water  quality  standards.  This  problem  of  radi ol ogically  containi nated 
water  could  persist  for  many  weeks. 

The  existing  radiation  detection  and  measuring  equipment  is  not 
capable  of  verifying  that  suspect  water  actually  meets  the  current 
water  quality  standards;  in  fact,  the  minimum  detectable  level  of 
radioacti vi ty  in  water  with  the  existing  equipment  and  procedures  is 
about  a  factor  of  1000  above  the  current  radiological  water  quality 
standard.  An  indication  of  acceptable  water  quality  can  be  obtained 
by  a  procedure  that  involves  measuring  the  activity  of  the  water  prior 
to  processing,  provided  the  decontamination  efficiency  of  the  water 
purification  system  is  known. 

The  water  purification  equipment  currently  under  development  can 
decontaminate,  to  the  current  radiological  water  quality  standards, 
water  from  a  watershed  contaminated  by  fallout  at  a  level  of  1  to  10 
R/Hr  at  H+l  hour.  After  about  4  days,  the  equipment  could  effectively 
decontami nate  water  from  a  watershed  contaminated  at  a  level  of  about 
100  R/Hr  at  H+l  hour.  The  water  purification  equipment  could  be  needed 
for  several  weeks  after  the  fallout  deposition. 
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INTRODUCTION 

This  appendix  provides  descriptive  information  on  water  sources 
within  the  Western  European  area  of  concern.  This  area  is  bounded  by 
Marburg  (Lahn),  Giessen,  and  Frankfurt  am  Main  on  the  west  and  the 
Fulda  River  valley  on  the  east.  The  information  is  used  to  identify 
typical  water  point  sites;  characterize  the  watersheds,  region,  and 
climate;  and  provide  an  approach  to  watershed  modeling. 

The  assistance  of  the  following  individuals  in  obtaining 
appropriate  references,  determining  applicable  factors,  and 
identifying  potential  resources  is  gratefully  acknowledged: 

Mr.  William  Abbe,  Terrain  Analysis  Center, 

Engineer  Topographic  Laboratories,  Ft.  Bel  voir,  VA 

Captain  Molzahn,  Department  of  Military  Engineering, 

U.  S.  Army  Engineer  School,  Ft.  Bel  voir,  VA 

Mr.  Henry  Zoller,  Library,  U.  S.  Geology  Survey, 

Reston,  VA 

Ms.  Laurie  Stackpole  and  Mrs.  Rouse,  Library,  National 
Oceanic  and  Atmospheric  Administration,  Rockville,  MD 

Mr.  Richard  Farnsworth  and  Dr.  Michael  Hudlow, 

Hydrology  Laboratory,  National  Oceanic  and  Atmospheric 
Administration,  Silver  Springs,  MD 

Mr.  Bruce  White  and  Mr.  Charles  White,  New  England  River 
Forecast  Center,  U.  S.  Weather  Service,  Bloomfield,  CT 
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SECTION  A-2 


WATER  SUPPLY  POINTS 

Water  supply  points  are  selected  to  support  troops  units  on  an 
area  basis.  Thus,  a  water  point  would  serve  all  water  supply  trucks 
and  trailers  coining  to  it  and  not  just  support  a  specific  unit.  The 
typical  field  water  purification  equipment  is  oriented  for  use  of 
surface  water  supplies.  The  Army  organization  includes  well  drilling 
capability  but  this  would  not  normally  be  used  in  an  area  such  as 
Germany.  The  Division  is  equipped  with  five  1500  gallon  per  hour 
water  point  units.  Additional  1500  gph  units  are  operated  by  Corps 

Combat  Battalions,  and  3000  gph  units  are  operated  by  the  Corps  Water 
Supply  Company. 

Water  supply  points  are  sited  on  the  bases  of: 

•  Adequate  Source,  minimum  of  2000-3000  gph  (to  allow  for 

waste  due  to  backflushing)  (equivalent  to  approximately 
2-1/2  liters  per  second)  per  1500  gph  unit.  These  values 
would  be  doubled  for  the  3000  gph  treatment  unit. 

•  Adequate  Road  Net,  to  permit  all  weather  vehicle  access  to 

the  fill  station  in  the  immediate  vicinity  of  the  treatment 
unit,  and  to  allow  for  turn  arounds  and  waiting  vehicles. 

•  Cover  and  Concealment,  forested  locations  are  desirable, 
placement  of  points  in  towns  or  in  the  vicinity  of  logical 
targets  would  be  avoided  (e.g.,  major  intersections,  bridge 
sites,  troop  concentrations). 

0  Good  Drainage  and  above  typical  flood  levels. 

•  Avoidance  of  Upstream  Sources  of  Contamination,  such  as 
industries,  major  towns,  sewer  outfalls,  or  stagnant  water 
(marshes,  swamps,  flooded  fields). 

The  above  doctrine  was  used  to  identify  typical  water  supply 
points  in  the  area  of  concern.  Three  points  were  identified  on  each  of 
10  map  sheets  (4  points  on  one  of  these  sheets),  to  permit 
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characterization  of  typical  points.  Other  points  could  have  been 
readily  identified,  due  to  the  amount  of  runoff,  number  of  streams  and 
excellent  road  nets  in  the  area.  Watersheds  which  would  provide  at 
least  5  liters  per  second  at  the  water  point  were  chosen  (approximate 
minimum  for  3000  gph  unit).  The  low  flow  values  of  1  to  2  liter  per 
second  per  square  kilometer  to  be  expected  in  the  region  indicated 
that  drainage  areas  of  at  least  5  square  kilometers  should  be  chosen. 
The  criteria  of  avoiding  upstream  sources  of  contamination,  where 
feasible,  and  the  distribution  of  towns  and  villages  indicated  that 
small  drainage  areas  be  chosen  when  available.  The  above  criteria  and 
considerations  of  vehicle  access  and  proximity  to  an  extensive  road 
network  were  used  in  identifying  the  water  supply  points  (WSP)  listed 
in  Table  A-l.  The  watershed  areas  were  delineated  on  copies  of  the 
identified  map  sheets. 


SECTION  A-3 


GENERAL  HYDROLOGIC  INFORMATION 

The  region  consists  of  rolling  terrain,  forested  or  under 
agriculture.  There  are  generally  sources  of  large  quantities  of 
surface  water  in  the  lower  part  of  the  area  of  concern  less  than  8 
kilometers  apart;  and  in  the  upper  part  of  the  area,  moderate 
quantities  of  surface  water  from  sources  less  than  16  kilometers 
apart.  There  are  generally  uniform  amounts  of  rainfall  throughout  the 
year  with  somewhat  more  rainfall  (by  30  to  50%)  occurring  in  the 
months  of  June,  July,  and  August  and  the  least  amounts  in  February, 
March,  and  April.  Due  to  the  differences  in  evaporation  losses, 
however,  the  largest  runoff  occurs  during  the  months  of  November 
through  April  (which  constitute  the  "Winter  Semester"  for  German 
hydrologic  studies).  This  as  shown  in  Figure  A-l  and  Table  A-2  for 
the  area  of  concern.  Summary  average  precipitation  data  for  selected 
stations  are  shown  in  Table  A-3. 

Typical  runoff  from  watershed  averages  from  5  to  15  liters  per 

second  per  square  kilometer,  with  the  lower  values  usually  associated 

with  smaller  watersheds.  Table  A-4  illustrates  mean  and  absolute  high 

and  low  runoff  rates  and  average  rates  for  several  watersheds  of  the 

Fulda  River  basin,  which  overlaps  the  area  of  concern.  Typically,  the 

flow  is  less  than  the  mean  flow  240  to  280  days  per  year,  due  to  the 

amount  of  runoff  directly  associated  with  storms.  The  flow  does  not 

exceed  50  to  60  percent  of  the  mean  flow  half  of  the  days.  Flow  less 

than  2  liters  per  second  per  square  kilometer  of  watershed  occur  very 

rarely,  and  for  streams  which  have  recorded  such  low  flows  the 

(A-l )* 

frequency  is  generally  less  than  2%  of  the  time. 


*TFe  number  in  the  parentheses  denotes  a  reference  that  is  identified 
in  Section  A-6. 
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Table  A -2.  Precipitation-runoff  audit,  Fulda  region. 


Jan 

Feb 

Entries 

in  inn 

of  Water  on 

the 

Drainage  Area 

Oct 

Nov 

Dec 

Mar 

Apr 

May 

Jun 

Jul 

Aug 

Sep 

N 

57 

50 

45 

55 

62 

72  . 

84 

76 

60 

63 

57 

59 

Av 

10 

12 

13 

13 

8 

6 

5 

5 

5 

6 

8 

9 

Ao 

22 

30 

23 

14 

5 

5 

6 

5 

5 

7 

17 

19 

V 

5 

5 

9 

28 

49 

61 

73 

66 

50 

24 

10 

5 

B 

+20 

+3 

-5 

-20 

-28 

-22 

-13 

-6 

-3 

+26 

+22 

+26 

A 

32 

42 

36 

27 

13 

11 

11 

10 

10 

13 

25 

28 

A/N 

.56 

84 

.80 

.49 

.21 

.15 

.13 

.13 

.17 

.21 

44 

47 

Ao/A 

.69 

.  71 

.64 

.52 

38 

,45 

.55 

.50 

.50 

.54 

.68 

.68 

Average  Runoff  in 

Liters 

Per 

Second 

Per 

Square 

Kilometer 

Av" 

3.  7 

5  0 

4.9 

5.0 

3.0 

2.3 

1.9 

1.9 

1.9 

2.2 

3.  1 

3.4 

Ao" 

8  2 

12.4 

8.6 

5.4 

1.9 

1.9 

2.  2 

1.9 

1.9 

2.6 

6.6 

7.1 

Legend 

N  -  Precipitation  (nun) 

Av  -  Ground  Water  Runoff  (mm) 

Aq  -  Surface  Runoff  (mm) 

V  -  Evaporation  (mm) 

B  -  Restoration  (+)  or  Depletion  (-) 

of  Ground  Water  Storage  (mm) 

A  -  Total  Runoff  (mm) 

A/N  -  Total  Runoff/Precipitation  Ratio 

Aq/A  -  Surface  Runoff/Total  Runoff  Ratio 

A  '  -  Ground  Water  Runoff  (f/s/km  ) 

A  "  -  Surface  Runoff  (;/s/km?) 
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able  A-3.  Selected  cl  i"i<i  to  logical  data  -  '.lest  Germany. 


As  illustrated  above,  the  most  runoff  is  surface  runoff;  the 
total  runoff  is  a  small  fraction  of  precipitation  most  of  the  year; 
and  precipitation  occurs  frequently  (on  almost  half  the  days) 
throughout  the  year,  but  on  only  an  average  of  one  or  two  days-  a  month 
is  the  precipitation  sufficient  to  cause  significant  runoff  occur. 

The  threshold  amount  of  precipitation  for  causing  runoff  varies 
during  the  year,  due  to  changes  in  temperature,  growing  vegetation, 
and  sunlight,  and  due  to  the  relative  permeability  of  the  surface  due 
to  prior  precipitation,  snow  cover,  or  frozen  ground. 

The  amount  of  area  covered  by  water  surface  is  a  small  fraction 
(1-3%)  for  the  region,  with  even  lower  values  typical  of  watersheds 
for  illustrative  water  supply  points. 

The  rivers  in  the  area  are  characterized  by  having  high  levels 
of  chemical  and  biological  contamination.  Lakes,  ponds,  and  small 
streams  contain  less  contamination,  however,  use  of  water  purification 
equipment  is  considered  essential.  Most  major  towns  and  industries 
are  situated  in  valleys  on  the  larger  streams  and  rivers.  There  are  no 
sizeable  reservoirs  or  lakes  in  the  area  of  concern.  Droughts  are 
rare,  but  even  in  the  most  recent  drought  period  (1959)  low  flow 
measurements  were  generally  above  0.5  liters  per  second  per  square 
kilometer  for  watershed  greater  than  5  square  kilometers. 


A- 14 


SECTION  A- 4 


WATERSHED  MODELING 

A-4.1  Occurrence  of  Surface  Runoff 

Most  stream  flow  is  due  to  surface  runoff,  however,  this  occurs 
only  intermittently.  The  principal  concern  is  the  likelihood  that 
radioactive  contamination  on  the  area  will  be  transmitted  to  the  water 
supply  point  (WSP).  Contamination  falling  on  all  of  the  stream 
surfaces  would  arrive  at  the  WSP  within  about  five  hours  for  most  of 
the  typical  WSP  watersheds.  The  average  time  for  the  contamination  to 
reach  the  WSP  would  be  approximately  half  this  value,  as  the  times 
developed  in  Table  A-l  are  based  on  the  arrival  at  the  WSP  from  the 
extreme  point  in  the  watershed.  The  area  of  running  water  in  the 
small,  typical  watersheds  would  generally  be  much  less  than  1%  (and 
possibly  less  than  0.1%). 

There  is  a  threshold  of  precipitation  which  will  result  in 
runoff.  This  threshold  varies  with  season  and  recent  rainfall.  New 
England  wa'  used  as  the  region  of  the  U.  S.  which  may  be  closest  in 
character  to  that  of  the  area  of  concern.  The  Appalachian  area  of 
Pennsylvania  may  also  approximate  the  region  and  could  similarly  be 
used  for  estimating  factors  not  immediately  available  for  the  German 
area.  Rainfall  and  runoff  records  for  New  England  have  shown  that  as 
little  as  0.15“  (approximately  4mm)  precipitation  in  the  optimum 
season  for  runoff  (February-March  in  New  England  and  Germany,  see 
Table  A-2)  may  cause  runoff,  as  noted  in  increased  stream  flow.  A 
typical  threshold  would  be  about  0.25"  (approximately  6mm),  while  0.5" 
(13mm)  may  be  necessary  under  very  dry  conditions  in  mid-summer  (July 
and  August,  see  Table  A-2).  Under  an  extreme  drought  condition  1.35“ 
of  precipitation  did  not  produce  an  increase  in  stream  flow. 

Another  factor  in  determining  the  occurrence  and  extent  of  storm 
runoff  is  recent  rainfall.  The  specific  factor  used  is  the  Antecedent 
Precipitation  Index,  which  relates  prior  precipitation  and  the  amount 
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of  time  since  it  fell.  In  summary,  the  saturation  effect  of  prior 

rain  is  assumed  to  be  reduced  by  a  factor  of  ten  percent  each  day. 

This  factor  and  the  seasonal  variations  are  used  with  rainfall 

(A-?  A-3  A-4) 

intensities  and  amounts  to  estimate  stream  flows.  ’  ’ 

The  general  occurrence  of  rain  in  the  region  of  concern  was 
summarized  in  Table  A-3.  Daily  precipitation  records  for  several  years 
for  the  area  are  given  in  Table  A-5  through  A-21.  These  data  permit 
the  determination  of  antecedent  precipitation  indices  for  individual 
days.  Typical  intensity  of  precipitation  in  the  storms  is  indicated 
in  Table  A-22.  Typical  occurrence  of  snow  as  a  component  of 
precipitation  is  illustrated  by  Table  A-23. 

A  simplified  relationship  for  use  in  the  runoff  analysis  for  the 
area  of  concern  is  necessary  due  to  the  limited  availability  of  stream 
hydrographs  in  which  precipi tation  and  subsequent  stream  flow  can  be 
directly  correlated.  Factors  used  for  the  simplified  relation  are 
those  typically  used,  but  developed  from  the  various  average  data  for 
the  area. 

An  antecedent  precipitation  index  (API)  factor  of  0.9  is  typical 

for  the  eastern  USA  and  has  been  recommended  as  an  appropriate  value 

( A-4 ) 

for  use  in  this  study.  A  threshold  precipitation  value  for  the 

occurrence  of  added  surface  runoff  (the  sum  of  the  API  and  the  day's 
precipitation)  of  from  4mm  to  I3mm  was  selected.  These  values  are 
based  on  the  experience  in  New  England.  The  threshold-time-of-year 
relation  selected  is  shown  in  Table  A-24  and  was  based  on  the  Fulda 
basin  precipitation  runoff  audit,  shown  in  Table  A-2. 

The  extent  of  precipitation  which  results  in  surface  runoff 
would  be  dependent  on  the  seasonal  surface  permeability  (partly 
reflected  in  the  above  threshold),  the  API,  the  short  term  intensity 
(e.g.,  Table  A-22),  and  in  winter  seasons,  the  extent  that  the 
precipitation  is  snow  (Table  A-23)  and  the  presence  of  snow  cover 
(approximately  mid-December  to  mid-March  in  the  Fulda  region.  Table 
A-3).  The  data  in  Table  A-23  indicates  that  most  of  the  precipitation 
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Table  A-5.  Daily  preci pi tation  data  (mm)  -  Bad  Hersfeld  (1972) 
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Table  A-6.  Daily  precipi tation  data  (mm) 


Bad  Kissinger)  (  1975) 


Day 

Jan 

rc5 

Oar 

Apr 

(lay 

Jill! 

Jul 

Au'i 

Sep 

Ccj. 

Nov 

pee 

1 

16  3 

9.0 

0.2 

0.4 

2.2 

■> 

11.5 

0.1 

2.6 

0.6 

5.7 

2.2 

4.6 

10.2 

3 

3.4 

5.3 

3.3 

1.6 

0.2 

0.3 

1.8 

0.6 

4 

2.5 

2.5 

0.5 

0.1 

0.8 

15.6 

1.5 

5 

7.9 

t 

2.5 

6 

6  3 

0  1 

.. 

0.5 

9.3 

7 

t 

0.9 

7.9 

7.3 

8 

t 

0.3 

5.8 

t 

0.3 

7.5 

9 

t 

13.6 

2.0 

10 

5.4 

2.0 

2.4 

0.5 

0.4 

13.7 

0.5 

11 

6  9 

1 . 5 

1.5 

5.1 

5.0 

12 

3  8 

5.9 

0.2 

1.4 

t 

0,3 

t 

C.6 

13 

7.  7 

6.6 

5.2 

0.  7 

t 

3.0 

6.5 

0.5 

14 

9.6 

0  2 

0.1 

t 

t 

0.9 

6.7 

0.8 

15 

2.2 

2.1 

0.  3 

t 

0.2 

16 

1,6 

0.6 

8.4 

t 

t 

17 

0.2 

4  7 

1,4 

t 

1.1 

0.3 

13 

t 

1.8 

t 

0.5 

t 

19 

2.3 

0.4 

1  -0 

t 

0.3 

1.0 

20 

10  0 

0,1 

4,8 

t 

21 

7.6 

0,7 

0.4 

t 

0.2 

;  2 

9.6 

0  3 

♦ 

0.1 

0.1 

2  3 

10  6 

1.6 

t 

1.5 

t 

24 

1,3 

2.4 

3  9 

12,1 

0.5 

16.3 

C  *J 

1.1 

t 

9  7 

3.2 

1.6 

0.4 

0-5 

26 

')  > 

t 

4.1 

0.2 

0.1 

27 

2 .  i 

0.1 

3.0 

0.2 

7.1 

1.2 

t 

26 

0,1 

8. 7 

t 

1.6 

0.5 

2  j 

0.6 

2.0 

0.1 

10.0 

t 

37 

t 

t 

3.9 

t 

0.7 

6.3 

10.7 

1  ) 

2.0 

0.8 

0.3 

0.3 

0.7 
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Table  A- 10.  Daily  precipitation  data  (mm)  -  Frankfurt  am  Main  (1976). 


f-'ar  Apr  p. 


10 

7  ] 

2  8 

1  9 

2  3 

0.  7 

12 

0  3 

17  6 

2.0 

13 

3  9 

0  5 

14 

0  7 

t 

15 

16 

1.0 

2.4 

17 

1 .  5 

1,2 

5,7 

H 

o  :> 

0  6 

15 

2  7 

0.2 

Jpn 

Jul 

An  c| 

Sep 
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fiov 

10  1 

0.9 

1.0 

0.2 

t 

0.2 

8.5 

2.7 

3.1 

4.4 

0.9 

0.6 

0.5 

1.8 

0.2 

0.9 

2.1 

t 

0.9 

0.5 

t 

9.2 

1.0 

t 

1,2 

t 

7.5 

3.5 

0.1 

t 

t 

1.3 

1.4 

0.4 

3.8 

2.3 

0.3 

0.9 

2  9 

3  6 

0.4 

0  4 

6.9 

2.0 

1.5 

5.2 

0.3 

3.7 

18  8 

t 

1.8 

t 

1.6 

0.2 

0.6 

1.9 

0.4 

19  9 

0.1 

3.6 

4  9 

0. 1 

o.: 

0  6 

1.2 

0.5 

0.2 

0.6 

3  7 

0.1 

0.4 

0.2 

1,6 

0.5 

0.3 

0.1 

0.0 

0.3 

1.8 

t 

1.4 

0.3 

2.8 

0.2 

2.1 

0.7 

t 

t 

823 
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Table  A-11.  Daily  precipitation  data  (min)  -  Frankfurt  am  Main  (1977). 


Vi: 

8.1'! 

fur 

‘ul  r 

Apr 

t)4y 

Jun 

Jul 

Auy 

Sep 

Oct 

llov 

1 

6  !' 

0.5 

2.3 

0.6 

3.0 

3.2 

1 

5.7 

0.4 

2.6 

1.3 

17.5 

3 

t 

L  -  * 

0.2 

7.0 

t 

0.2 

0.3 

12.5 

4 

7.6 

0.4 

1.4 

t 

1.5 

4.5 

2,4 

5 

0  1 

1.0 

1.6 

4.9 

16.0 

6 

0,  4 

4  9 

t 

0.6 

7 

3-  3 

1.1 

7.2 

3.6 

8.2 

0.3 

0.1 

8 

0.3 

0-3 

t 

3.1 

3.2 

0.3 

9.0 

0.1 

9 

5  1 

4.5 

2.0 

1.3 

1.4 

10 

6.;. 

9  6 

t 

0.1 

4.6 

0.1 

0.1 

11 

1.8 

1.7 

0.8 

4.9 

12 

11,3 

1.9 

5.3 

0.1 

3.4 

4.4 

13 

4.1 

0.  3 

1.4 

1.8 

1.3 

2,8 

40.2 

0,3 

14 

3  1 

0-1 

4.2 

1.5 

1.1 

0.2 

8,6 

15 

1.3 

58 

0.2 

0,5 

2.1 

5.2 

3.0 

16 

t 

0  - 

2.2 

6.2 

17 

- 
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22.2 

0.1 

0.1 

13 

.4 

2.7 

5.4 

6.7 

0.2 

0.4 

t 

19 

8.1 

:• 

1.9 

3.5 

31.6 

30 

?  ■  1 

0.6 

0.9 

4.1 

12.9 

0.8 

0.5 

31 

1.8 

1.1 

0.5 

t 

8.3 

5.0 

22 

t 

2.1 

0.) 

t 

24.0 
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1.3 

33 

0.  8 

11  6 
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5,0 

?: 

1.1 

*>  * 
i  ,  " 

t 

t 

0.7 

20,2 

2.1 

1.2 

1  > 

* 

t 
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6.6 

3.6 
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7.2 

1.7 

( 

8.1. 

t 
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4.3 

7.6 

1.3 

t 

2  7 

o. : 

11.7 

2.7 

0.9 

0.2 

f'i 

t 
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9.8 
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2.4 
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Table  A- 12.  Davy  precipitation  data  (mm)  -  Frankfurt  am  Main  (1  980) 


Dj  y 

Jan 

rob 

Mar 

Apr 

May 

Jt:il 

Jul 

Au  q 

Sep 

C:t 

r;-;v 

1 

t 

t 

t 

0.7 

t 

1.3 

No 

t 

2 

1,3 

t 

0.3 

0.4 

0.  3 

Du  t.  a 

3 

t 

2.0 

t 

1.2 

0.1 

for 

t 

4 

0,8 

1.4 

t 

t 

0.1 

MvM.h 

t 

v. 

5 

0  1 

0,6 

t 

0.5 

t 

6 

0,  1 

0.2 

0.9 

0. 1 

0  H 

0  7 

0.3 

t 

7 

0.2 

0.1 

0.  2 

t 

t 

0.5 

0.7 

0.9 

8 

t 

0.1 

0.  3 

0  7 

1.2 

0.3 

9 

t 

t 

0.1 

0  4 

0. 1 

0.3 

0.7 

10 

0  1 

0  4 

t 

0.7 

l;  !• 

t  ’ .  7 

11 

t 

0,  1 

t 

0  • 

t 

0.7 

12 

t 

0  5 

0.' 

0.  3 

13 

0  1 

G  1 

t 

0  2 

f  .1 

14 

1.0 

1.6 

0  : 

18 

0.3 

t 

n  ;> 

0.2 

0. 

16 

t 

t 

t 

0.  3 

1  . 

17 

0  1 

t 

0.7 

u  . 

18 

t 

t 

0  3 

19 

t 

t 

0.  1 

0.7 

0.2 

0. 

20 

t 

f 

0.1 

0  1 

1  . 

21 

0,  7 

0,7 

0.1 

t 

t 

t 

22 

0.  3 

0.1 

t 

0.4 

t 

23 

0.6 

t 

0.6 

24 

t 

0.  3 

1.3 

0.1 

O.o 

25 

0  1 

0.1 

1.2 

0.2 

0. 1 

t 

O.f 

76 

t 

0  9 

0  1 

0.4 

27 

t. 

U.  !■ 

0  1 

0.4 

0.1 

t 

t 

0.' 
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t 
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t 

t 
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2 .) 

t 

t 
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0. 1 

3.5 
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0.1 

O.i 

30 

o .  r> 

t 

t 
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0.  3 

t 

t 

i 

31 
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t 
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0.7 
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Table 


-13.  Daily  precipitation  data  (mm)  Fulda  (1975) 


t 


?i.  1 


i 


r,  ‘i 

0  5 


0.  I 
i  a 


i; 


33  7 

t  4.3  t 

t 

34.') 

;>  33  0  0 

3.3 

3,8  9  4 


0  ! 

14.3  0.0 

r  4  4 


0.0  0.7 

0.3 


0  ' 
1  3 
4  t 


0  7 
If-  7 


1  l) 
If.  • 
r 


1  ] 

9.6 

7.6 

1.7 


3.6 


0  . 


7.9 

16.0 

0.6 

1.3 

t 

0.1 

74.  . 


A- 2  5 


Table  A- 1 4.  Daily  precipitation  data  (mm)  -  Fulda  (1976). 


I'-jv 


3 

4 

6 

7 

8 

9 

10 

II 

1? 

1 3 

14 

15 

16 
17 
IB 

19 

20 
21 
22 

23 

24 

25 

26 
2  7 
7S 
?) 

If! 

?! 


,'.111 

F.-h 

Mil  r 

Apr 

h'Td 

13.2 

8.0 

t 

t 

5.9 

5.6 

1.3 

t 

0.2 

2.5 

t 

t 

1.5 

0.2 

1.8 

3  1 

2-3 

3. 1 

1.8 

0.  3 

4  3 

t 

2  6 

1.5 

0.1 

4  9 

0  2 

0.1 

£ 

4.4 

1.6 

0.9 

1.5 

0.2 

3  2 

t 

0.7 

2,2 

t 

1.4 

8  4 

3  1 

6  6 

3.  6 

0  7 

0.5 

t 

5.  2 

4.2 

2.  1 

i.n 

8.  7 

1.0 

0. 

0-7 

t 


!’:£Z 

Jun 

19  3 

Jul 

Auij 

6.1 

0.2 

2, 1 

1.3 

1.9 

4.9 

3.8 

0.5 

3.1 

6.6 

1.1 

1.8 

3.9 

25-3 

0,1 

0,8 

8,1 

t 

10.1 

t 

3.2 

0.2 

t 

4,9 

5.9 

4.4 

1.4 

0.8 

1.0 

0  7 

1.0 

t 

1.3 

fi,8 

0.2 

1 .8 

0.9 

t 

2.2 

5.6 

t 

Srp 

Oct 

flaw 

IV  r 

2. 4* 

oil 

6.5 

'  t 

5.9 

3.4 

1.1 

3.0 

0.3 

0.1 

1.6 

t 

0.2 

9.7 

1.6 

1.0 

0.2 

5,8 

7.6 

5.0 

t 

1.0 

3.1 

6.6 

1.1 

0.8 

1.0 

0,8 

5.5 

4.7 

t 

0  5 

1.1 

7.5 

0.  1 

0.4 

0.6 

0  1 

0  3 

2.1 

0.3 

0,3 

0.3 

t 

t 

2.3 

0.2 

t 

0.1 

t 

0.3 

2.0 

t 

0.4 

t 

t 

2.3 

t 

0,8 

0  2 

1.0 

t 

3,0 

0.2 

3.  1 

0.1 

1.0 

1.0 

0.2 

0.1 

8.7 

0.9 

15. 1 

t 


/\-26 
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Table  A- 

15.  Daily  precipitation  data 

(mm) 

-  Fulda 

(1977) 

Day 

j.ti 

Feb 

Mar 

Apr 

May 

Jun 

Jul 

Aug 

Sep 

Oct 

Nov 

Dec 

1 

1.3 

6.3 

1.9 

7.3 

2 

3.8 

0.7 

,  t 

0.5 

11.3 

3 

0  7 

1-1 

0.4 

3.4 

1.0 

0.2 

1.5 

31.2 

4 

t 

2.7 

1.6 

1.4 

t 

0.7 

3.1 

4.0 

5 

1.3 

t 

3.4 

8.7 

14.6 

t 

6 

1.6 

1.9 

3.0 

0.6 

1.2 

3.4 

7 

0.1 

7-0 

1.7 

0.3 

2.5 

0.1 

is.: 

2.8 

t 

8 

3.3 

t 

1.5 

0,6 

0.2 

0.3 

6.7 

0.6 

9 

2  0 

0.2 

t 

6.6 

0.3 

3.7 

0.4 

1.1 

10 

9.  0 

5  0 

t 

9.0 

t 

3.4 

t 

11 

0.3 

2.2 

1.1 

t 

t 

t 

0,8 

12 

3,4 

2 . 6 

2.1 

3.0 

5.2 

1.0 

13,5 

8.7 

13 

3. 1 

0.5 

4.0 

3.6 

2. 1 

6.3 

0.2 

0.7 

0.6 

14 

4.1 

7  0 

2.7 

16.2 

0.5 

14.3 

t 

15 

1.6 

0.3 

1.4 

2.3 

0.6 

5.5 

0.2 

0.6 

0.2 

5.5 

0.4 

16 

2,3 

O.c 

t 

0.2 

10.3 

0.1 

2.5 

17 

0.2 

0.2 

0.6 

23.7 

2.1 

18 

2.6 

10.2 

6.4 

0.3 

1.4 

19 

3  2 

13  8 

13.2 

8.4 

10.7 

0,7 

20 

L 

29.0 

1.0 

1.9 

3.3 

t 

21 

1.6 

0.  3 

t 

6.1 

0.2 

4.0 

22 

2.4 

t 

0  4 

5.0 

0.3 

0  6 

23 

0.2 

ua 

t 

0.1 

6.7 

1.1 

24 

0.2 

1,3 

30.4 

7.0 

9.9 

0.7 

25 

6.6 

1.3 

0  K 

2.3 

4.4 

3.2 

t 

3.7 

2.2 

26 

1 .?. 

2.7 

5.6 

0.9 

9.9 

1.2 

0.9 

t 

0.7 

2.1 

27 

2.6 

0.6 

11.6 

t 

t 

0.5 

f..> 

3.6 

1.1 

t 

t 

t 

IS.! 

0.7 

.  .» 

?  .  '« 

t 

0.1 

3.<> 

3.  ■: 

2.4 

1.5 

0.9 

3.4 

J ! 

t 

t 

6. 6 

1.2 

r\  ■>  7 
M-C  / 
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Table  A-17.  Daily  precipitation  data  (mm)  -  Giessen  (1975) 


D.) , 

Or: 

1-Vi 

r-'.ii- 

Apr 

Hay 

Jun 

Jul 

Aun 

Sep 

Oct 

Nov 

Dec 

i 

t 

0  4 

t 

_  t 

1.4 

672 

n 

1.3 

12.1 

5.6 

2.4 

3.9 

7.7 

3 

0.  1 

12.6 

1.5 

2.7 

t 

1.3 

1.0 

4 

t 

2.2 

0.4 

0.3 

4.2 

t 

5 

t 

0.5 

"3.7 

0.5 

1.0 

0.2 

6 

0,  2 

0.6 

5,1 

1.3 

0.9 

1.3 

7 

8,2 

4.0 

2.1 

4.0 

0.3 

0.3 

0.8 

8 

0,9 

4.0 

3.4 

t 

9 

t 

0.1 

10 

6,2 

0.1 

1.0 

3.7 

4.2 

0.2 

t 

11 

0.4 

t 

0-6 

2.3 

3.4 

0.3 

0.4 

12 

0  3 

0,8 

0.9 

5.0 

6.1 

2.3 

13 

t 

2.1 

1.0 

t 

7.7 

7.2 

1.7 

0.  7 

6 . 2 

14.6 

1.2 

0.5 

1.0 

1 5 

0  I 

0.8 

7.5 

0.8 

3.0 

t 

C,  5 

H. 

0,9 

1.3 

0.7 

t 

8.0 

t 

1  ’ 

1.3 

0.3 

l.C 

0.6 

15.4 

’2.7 

0.3 

6,6 

4.3 

H 

2  3 

16,2 

0  4 

2.6 

10.1 

t 

0.7 

4.0 

0.3 

0.3 

13 

0. 1 

1.4 

8.5 

2.9 

0  2 

t 

1.5 

0.5 

3.6 

20 

0  7 

0.1 

0  9 

26.4 

2.4 

0.2 

0,2 

21 

t 

4.6 

2.0 

6.5 

1.1 

0.1 

22 

5.4 

t 

35.4 

2  3 

0,6 

t 

24 

2.7 

2.8 

0  4 

2-0 

L  J 

7,6 

1.1 

0.6 

3.0 

0.5 

5.3 

0.1 

1  5 

crj 

0.3 

2.3 

0.1 

7.7 

3.1 

21 

11.9 

8  7 

0.9 

0.1 

1.5 

yn 

6.3 

0.6 

2.3 

2.3 

1.3 

t 

2.8 

0.  7 

4. 1 

v;: 

1.5 

t 

0.5 

0.6 

6.7 

1.0 

4.5 

0.5 

1.7 
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Table  A-19.  Daily  precipitation  data  (mm)  -  Giessen  (1977). 


v.1 .3  r  '• 

TeO 

Mar 

Apr 

(■.ay 

Jun 

Jul 

Auij 

Sep 

Oct 

fiov 

[v-C 

9  3 

t 

4.9 

o’.T 

T.o 

4.3 

9  ti 

0.5 

1.0 

t 

3.2 

17.2 

0.3 

1.9 

0,1 
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Table  A-23.  Percentage  of  frequencies  of  occurrence 
of  mixed  precipitation  at  Fluda. 


Hours 

Thunder¬ 

Rain  and/or 

Freezing  Rain 

Snow  and/ 

Percent  of  Observation 

(LST) 

storms 

Drizzle 

and/or  Drizzle 

or  Sleet 

with  Precipitation 

Jan 

00-02 

0.0 

7.7 

0,5 

12.1 

20.3 

03-05 

0.0 

13.8 

0.7 

10.8 

25.3 

06-08 

0.0 

17.0 

0.4 

10.7 

28.0 

09-11 

0,0 

15.3 

1.6 

9.1 

25.7 

12-14 

0,0 

13.3 

0.9 

8.6 

22.8 

15-17 

0,0 

12.1 

0.6 

8.5 

21.2 

18-20 

0.0 

11.0 

0.5 

14.3 

25.8 

21-23 

0.0 

5.4 

0.0 

19.4 

24,7 

Totals 

0.0 

13.7 

0.8 

10.2 

24.6 

Apr 

00-02 

0.0 

14.8 

0.0 

2.8 

17.6 

03-05 

0.0 

14.2 

0.0 

3.3 

17.2 

06-08 

0.0 

14.0 

0.0 

4.1 

17.7 

09-11 

0.0 

14.7 

0.0 

3.5 

17.6 

12-14 

0.5 

12.6 

0.0 

2.1 

14.7 

15-17 

0,4 

15.5 

0.0 

1.7 

17.1 

18-20 

0.6 

11.9 

0.0 

1.7 

14.1 

21-23 

0,0 

10.0 

0.0 

3.3 

13.3 

Totals 

0,2 

13.9 

0.0 

3.0 

16.7 

Jul 

00-02 

0.6 

3,  3 

0,0 

0.0 

3,9 

03-05 

0,3 

4.8 

0,0 

0.0 

5,0 

06-08 

0.3 

7.  7 

0.0 

0.0 

7.7 

09-11 

0.5 

6,8 

0.0 

0.0 

7.1 

12-14 

1.8 

6,3 

0.0 

0.0 

7.4 

15-17 

2.0 

6  6 

0,0 

0.0 

8.1 

18-20 

2.2 

3.4 

0.0 

0.0 

5,1 

21-23 

0.0 

3,2 

0.0 

0.0 

3,2 

Totals 

1.0 

6  1 

0.0 

0.0 

6.7 

Oct 

00-02 

0.0 

12.3 

0.0 

0.0 

12.3 

03-05 

0.0 

13.0 

0.0 

0.3 

13.2 

06-08 

0.0 

11.9 

0.0 

0.0 

11.9 

09-11 

0.0 

10ol 

0.0 

0.0 

10.1 

12-14 

0.1 

11.5 

0.0 

0.3 

11.6 

15-17 

0.4 

12.5 

0.0 

0.4 

12.6 

18-20 

0.0 

9.5 

0.0 

0.0 

9.5 

21-23 

0.0 

10.8 

0.0 

0.0 

10.8 

Totals 

0.1 

11.6 

0.0 

0.2 

11.7 
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Table  A-24.  Cummulative  precipitation  necessary 
for  added  surface  runoff. 


Month 

Thresh' 

January 

8 

February 

4 

March 

5 

Apri  i 

9 

May 

12 

June 

12 

July 

13 

August 

13 

September 

12 

October 

12 

November 

10 

December 

9 

in  the  area  occurs  as  rain,  even  in  mid-winter.  The 
precipitation-runoff  audit.  Figure  A-l,  shows  that  snow  cover  has 
little  inhibiting  effect  on  runoff,  the  fraction  of  precipitation 
which  becomes  surface  runoff  being  greatest  in  the  winter  months. 
Snow  cover  may  have  a  delaying  effect  upon  when  the  actual 
precipitation  enters  the  streams,  however,  as  the  surface  runoff  would 
consist  of  both  melted  snow  from  prior  precipitation  and  recent 
rai nfal 1 . 


The  single  seasonal  precipitation-surface  runoff  ratio  selected 
to  be  used  for  this  study  is  shown  in  Table  A-25.  The  relationship  is 
an  approximation  as  the  proportional ity  would  actually  be  a  function 
of  intensity  and  total  precipitation  in  the  individual  storm.  The 
entries  in  Table  A-25  were  estimated  by  considering  the  general 
frequency  of  precipitation  sufficient  to  directly  increase  surface 
runoff  (about  1  per  month  in  winter  to  2  per  month  in  summer.  Table 
A-13),  the  average  rainfall  in  the  month  and  the  average  surface 
runoff  (Table  A-2),  and  subtracting  the  thresholds  (Table  A-24)  from 
the  precipitation. 


The  factors  developed  above  can  be  used  to  provide  a  rough 
approximation  of  surface  runoff  which  might  occur  following  a  given 
level  of  precipitation  on  a  day  of  a  given  month,  with  the  additional 
data  of  recent  prior  daily  precipitation.  This  relationship  is  shown 
below: 

SRO  -  [SR0/Pg]  x  Pe 


wi  th 

and 

where 


P  =  P  +  API  -  CPT 
e  o 

API  =  IP0_n  x  (0. 9)n 

SRO  -  surface  runoff  (mm) 

[SRO/P  ]  =  surface  run-off  excess  precipitation  ratio, 

'  given  in  Table  A-25 

Pg  =  excess  precipitation  (mm) 

PQ  =  24-hour  precipitation  (mm);  given  in  Table  A-13 
API  =  antecedent  precipitation  index 
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Table  A-25.  Estimated  proportion  of  excess  rainfall 
which  becomes  surface  runoff. 


Month  Surface  Runoff/Excess  Rainfall 

January  0.45 
February  0,65 
March  0  58 
April  0  34 
May  0  1 1 
June  0  10 
July  0.13 
August  0.10 
September  0.12 
October  0.16 
November  0.40 
December  0  42 
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CPT  =  cumulative  precipitation  threshold  (mm),  given 
Table  A-24 

PQ  n  =  24-hour  precipitation  n  days  before  the  day  of  PQ 

(but  if  Po-n>  CPT  use  Po-n  =  CPT,  to  avoid  double 
counting  of  prior  precipitation  that  resulted 
surface  runoff),  as  given  in  Table  A-13. 

The  above  expression  provides  the  possibility  of  double  counting 
24-hour  precipitation  below  the  level  of  the  CPT  which  occurred  prior 
to  both  the  current  storm  and  an  earlier  storm  which  resulted  in 
excess  precipitation.  Further  refinement  of  the  expression,  however, 
does  not  appear  justifed  in  view  of  the  assumptions  and  approximations 
associated  with  the  factors  it  contains. 

Application  of  the  expression  is  illustrated  by  the  following 
example. 

Problem:  Determine  surface  runoff  (SRO)  from  a  watershed  in 

the  Fulda  river  basin  region  due  to  the  storm  of 
14  July  1975. 

Solution:  [SR0/Pe]  =  0.13  (Table  A-25) 

CPT  =  13  mm  (Table  A-24) 

PQ  =  33.7  mm  (Table  A-13,  14  July) 

API  =  [.5  x  .92]  +  (Table  A-13,  12  July) 

[6.9  x  . 93]  *  (Table  A-13,  11  July) 

(prior  precipitation  ignored) 

API  =  5.4  nin 

P  =  33. 7+5.4-13=26  mm 
e 

SRO  =  0.13  x  26 
SRO  =  3.4  mm 

2 

If  the  watershed  area  is  10  Km  ,  the 
surface  runoff  passing  the  water  point 
resulting  from  the  precipitation  on 
July  1975  would  be: 

(3.4)mm  x  (10)km2  x/  1  \  m  x  (1000)2  jn2  =  34000m2. 

'1000/  mm  kr12 
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A  4.2  Time  Factors  in  Surface  Runoff 

The  expression  developed  above  does  not  indicate  the  time, 
following  the  occurrence  of  precipitation,  when  the  surface  run-off 
would  pass  the  point  of  ex' t  of  the  watershed  (that  is  the  water 
supply  point).  This  time  is  a  function  of  the  size,  shape,  and  average 
slope  of  the  watershed;  surface  conditions  (time  of  year);  and  the 
parameters  of  the  storm  iti.’f. 

This  study  is  primarily  concerned  with  small  watersheds  (Table 

A-l)  with  negligible  surface  storage.  Surface  runoff  resulting  from 

rain  on  the  watershed  which  exceeds  the  cumulative  precipitation 

threshold  (Table  A-24)  would  start  to  arrive  at  the  water  supply  point 

(WSP)  immediately  from  the  area  draining  immediately  upstream.  The 

amount  arriving  from  a  continuing  rainfall  would  increase  rapidly 

until  the  "time  of  peak  flow"  was  reached,  which  corresponds  to  the 

longest  time  that  it  takes  runoff  originating  anywhere  on  the 

watershed  to  arrive  at  the  WSP.  Snyder's  Synthethic  Procedure  can  be 

used  to  develop  this  longest  time  (also  referred  to  as  the  "time  of 
( A-5 ) 

concentration" ) .  In  this  procedure  the  time  is  related  to  the 

watershed's  relief,  shape,  and  length  by  the  following  expression: 

tp  =  CT  •  (L  •  Lc)  03 

where 

tp  =  time  of  concentration  (hr), 

Cy  -  constant, 

L  =  length  of  watershed  (miles),  and 

Lc  =  distance  from  the  center  of  gravity 

watershed  to  the  watershed  outlet  (water 
point)  (miles). 

Values  of  Cy  (see  Table  A-l)  are  based  on  those  for  the  Appalachian 

Mountain  regions,  modified  for  the  relative  average  steepness  of  the 

watershed  slope.  The  calculated  values  of  t  are  shown  in  Table  A-l. 

P 
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The  rate  of  runoff,  from  rain  falling  at  a  uniform  rate  and 
considering  such  factors  as  velocity  of  sheet  and  stream  flow,  would 
increase  with  time  as  the  ground  absorbed  less  of  the  precipitation 
and  flows  became  greater.  As  a  result,  the  amount  of  flow  at  the  WSP 
would  continue  to  increase  beyond  the  time  of  concentration  ("time  to 
peak  flow"),  at  least  until  the  precipitation  rate  reduced. 

The  flow  of  surface  runoff  following  a  rainfall  decreases 
approximately  exponentially.  This  is  illustrated  for  typical  streams 
in  Figure  A-2,  and  some  specific  rivers  in  Figure  A-3.  The  gradual 
reduction  in  flow  following  heavy  rain  for  points  in  the  Fulda  River 
basin  is  shown  for  three  cases  in  Figures  A-4,  A-5,  and  A-6.  Note 
that  the  stream  records  are  in  terms  of  gauge  height.  The 
precipitation  records  are  cumulative,  thus  the  periods  of  most  intense 
rainfall  are  where  the  precipitation  curves  have  the  greatest  slope. 
Figure  A-4  shows  the  results  of  one  winter  storm.  Figures  A-5  and  A-6 
illustrate  stream  heights  during  and  following  two  successive  storms. 


In  consideration  of  the  above  and  the  relatively  small  areas  of 
the  WSP  watersheds  (Table  A- 1 ) ,  it  may  be  assumed  that  surface  runoff 
occurs  essentially  within  four  days  of  the  rainfall.  An  assumed 
distribution  of  this  flow  is:  first  day,  38%;  second  day,  44%;  third 
day,  11%;  and  fourth  day,  7%.  To  approximate  flows  from  successive 
days  of  excessive  precipitation,  the  flows  can  be  accumulated.  For 
convenience  of  analysis  in  this  effort,  the  flow  may  be  assumed  to  be 
distributed  over  the  entire  24  hours.  Thus  surface  runoff  of  1  mm  in 
one  day  would  be  equal  to  a  flow  of: 


(1) 


mm 

day 


_ ni 

mm 


x  (1000)  y- 

ITl 


x  (1000) 2 

Km 


11.6 


■'  /  sec 
2 

Km1" 


The  approximate  total  flow  would  be  the  combination  of  ground 
water  runoff  for  the  corresponding  month  (Table  A-2)  and  the  surface 
runoff  caused  by  the  precipitation. 
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a.  Red  River  Basin  above  Denison,  Texas 

b.  Delaware  River  Basin  above  Port  Jervis,  New  York 


c.  Susquehanna  River  at  Tawanda,  Pennsylvania 


££  >y  rc<r>f£H  m  inct>€S 
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Figure  A-3.  Actual  U.  S.  storm  hydrographs 
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In  the  case  where  radioactive  contamination  occurs  within  the 
fourth  day  of  excess  precipitation  (which  produced  surface  runoff), 
the  residua)  flow  would  produce  some  immediately  contaminated  stream 
flow.  The  area  from  which  such  flow  would  occur  could  be  expected  to 
be  less  than  even  the  fractional  amount  of  the  runoff,  as  the  entire 
surface  area  generally  only  contributes  runoff  while  it  is  raining. 
Subsequent  areas  from  which  drainage  is  occurring  are  temporary  ponds 
and  lower  areas  of  sheet  flow.  In  consideration  of  the  Fulda  River 
basin  topography  and  to  provide  a  basis  for  the  study,  use  of  the 
following  fractions  of  the  surface  area  as  contributors  of  surface 
runoff  are  recommended:  day  of  precipitation,  100%;  second  day,  40%; 
third  day,  8%;  and  fourth  day,  2%. 

A-4.3  Volume  of  Water  on  Watershed 

An  approximate  volume  of  surface  water  on  the  WSP  watershed  is 
required  for  contamination  and  dilution  assessment.  The  approximation 
developed  below  is  based  on  prior  estimates  of  normal  and  post-rain 
area  of  surface  water,  runoff  of  storm  flow  over  time,  ground  water 
flow,  and  assumption  of  average  stream  velocities. 

The  volume  of  water  on  the  watershed  directly  associated  with 
rainfall  depends  on  the  volume  of  rain,  the  proportion  of  rain  which 
becomes  surface  runoff,  and  the  duration  and  shape  of  the  storm 
hydrograph.  An  approximate  flow  distribution  and  time  were  assumed 
above.  This  is  repeated  in  Table  A-26  with  the  presumed  area  of 
watershed  covered  by  surface  water.  It  is  assumed  that  the  rain 
occurs  on  only  the  first  day.  Rain  on  the  successive  days  or  prior  to 
complete  recovery  from  a  prior  precipitation  should  be  treated  as 
discussed  in  conjunction  with  the  antecedent  precipitation  index. 

For  the  purpose  of  this  study  and  with  recognition  that 
approximations  are  sought  which  would  characterize  the  region  and  not 
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Table  A-26.  Storm  runoff  and  area  coverage  factors. 


I  tem 


Distribution  of  Storm  Runoff 
(day  of  passing  the  WSP) 


Area  Covered  by  Water 

(area  contributing  contamination 

to  stream  via  surface  water) 


Day  of 

2nd 

3rd 

4th 

Rai  n 

Day 

Day 

0..33 

0.44 

0.11 

0  07 

1 .00 

0.40 

0..80 

0  02 
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necessarily  be  exact  for  one  specific  watershed,  it  is  assumed  that 
the  storm  surface  runoff  and  the  runoff  due  to  ground  water  may  be 
treated  independently  and  are  additive.  The  principal  considerations 
associated  with  this  assumption  relate  to  definitions  of  the  types  of 
runoff  and  stream  flow. 

There  are  several  forms  of  surface  runoff  from  a  storm. 

Precipitation  on  flowing  water  has  an  immediate  impact,  however,  the 
surface  area  is  a  very  small  percentage  of  the  total.  Some  of  the 
storm  runoff  flows  directly  through  the  stages  of  sheet  flow  into 
intermittent  (ephemeral)  streams  (essentially  carrying  only  surface 
runoff  and  therefore  only  existing  while  there  is  post-storm  drainage) 
into  perennial  streams  (which  carry  ground  water  and  surface  runoff 
and  thus  flow  essentially  all  of  the  time).  Other  storm  runoff  is 

"interflow"  which  is  surface  runoff  which  is  delayed  due  to  ground 
surface  conditions  such  as  heavy  vegetation,  forest  floor  cover,  or 

very  localized  ponding. 

Interflow  may  also  occur  as  lateral  flow  just  below  the  surface 
which  joins  the  stream  flow  in  time  to  form  part  of  the  storm 
hydrograph  (the  runoff  above  that  which  would  have  occurred  had  the 
storm  not  taken  place).  The  ground  water  flow  is  generally  somewhat 
elevated  after  the  storm  discharge  (4  days  for  the  watersheds  of  this 
study),  which  can  be  attributed  in  part  to  slower  interflow.  The 

ground  water  recession  is  a  much  longer  and  continual  process,  as 
ground  water  runoff  is  dependent  on  periodic  recharging  the  ground 
water  level  by  precipi tation.  Another  form  of  direct  ground  water 
influence  on  the  storm  hydrograph  results  from  the  temporary  sharp 
rise  in  the  water  table  adjacent  streams  carrying  increased  flow.  The 
higher  water  level  can  raise  the  water  table  level  in  the  banks 
sufficiently  to  reverse  the  hydraulic  gradient  to  away  from  the 
stream,  surcharging  the  adjacent  ground  and  deferring  the  normal 
ground  water.  The  impact  of  this  temporary  in-ground  storm  storage  is 
to  prolong  the  storm  hydrograph  and  to  raise  the  post-storm  ground 
water  flow.  That  flow  of  water  that  entered  the  channel  as  stream 
runoff  and  exited  the  ground  adjacent  to  the  stream  banks  during  the 
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period  of  the  storm  hydrograph  can  be  considered  to  be  surface  runoff 
for  the  purposes  of  this  study. 

The  volume  of  storm  water  on  the  surface  of  a  watershed  may  be 
determined  from  the  definition  of  surface  runoff  (which  is  inexact), 
the  volume  of  precipitation  and  the  rate  of  runoff.  The  volume  of 
surface  runoff  is  the  fraction  of  the  excess  precipitation  which 
results  in  runoff.  Table  A-24  and  Table  A-25  provide  average  values 
by  months  of  the  year  of  the  cummulative  precipitation  necessary  for 
added  surface  runoff  and  the  proportion  of  that  excess  rainfall  which 
runs  off.  These  provide  the  expression  for  surface  runoff  (SRO): 

SRO  =  (SR0/Pe)  x  (P  +  API  -  CPT) 

where  SR0/Pe  is  from  Table  A-25,  CPT  is  from  Table  A-24  and  the 
antecedent  precipitation  index  (API)  is  determined  by  the  time  and 
amount  of  the  preceding  precipitation. 

The  volume  of  water  on  the  watershed  due  to  precipitation 
occurring  on  day  one  may  be  calculated  from  the  value  of  SRO  (in  mm) 
derived  using  the  expression  above  and  the  factors  in  Table  A-26.  The 
result  is  shown  in  Figure  A-7.  There  is  surface  water  on  the 
watershed  during  and  immediately  following  a  storm  which  does  not  run 
off  as  surface  runoff,  but  enters  the  ground  or  is  evaporated  (a  very 
small  proportion).  This  surface  water  is  not  included  in  Figure  A-7. 
The  amount  of  precipitation  which  is  at  the  surface  and  then 
immediately  or  later  enters  the  ground  is  essentially  the  total 
precipitation  (P  )  less  the  surface  runoff  (SRO).  This  may  be  derived 
from  the  equation  shown  in  the  preceding  paragraph. 

The  volume  of  ground  water  that  is  on  the  watershed  surface  at 
any  time  is  dependent  on  the  ground  water  flow  (Table  A-2),  and  shape, 
topogaphy,  and  size  of  the  watershed.  (summarized  in  Table  A-l  for 
the  representative  WSP).  This  volume  is  akin  to  "Channel  Storage", 
except  that  the  term  is  associated  with  the  excess  volume  of  water  in 
the  streams  due  to  a  storm  (part  of  the  volume  shown  in  Figure  A-7). 
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Cubic  Meters  Surface  Water  on  Watershed 
per  Square  Kilometer  per  1mm  SRO 


P 

o 


Total  Precipitation 

(Storm  assumed  to  occur  in  6-hour  period) 


Time  After  Storm  (days) 


Figure  A-7.  Volume  of  storm  runoff  on  watershed. 
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The  volume  in  the  streams  when  all  flow  is  due  to  ground  water  (i.e., 
after  the  4th  day  following  a  storm  for  the  WSP  watersheds  being 
studied)  can  be  equated  to  the  ground  water  flow  and  average  time  that 
it  takes  emerging  ground  water  to  exit  the  watershed.  The  time  of 
concent  rat  ion  (or  lag  time)  discussed  above  and  shown  in  Table  A-l  is 
an  indicator  of  the  time  that  emerging  ground  water  would  take  in 
transit  to  exit  the  watershed.  (Note:  it  is  assumed  that  the 
watershed  boundaries  apply  to  both  surface  and  ground  water  runoff. 
This  is  not  necessarily  the  case  but  these  boundaries  generally 
coincide  in  the  Fulda  River  basin).  The  time  of  concentration  may  not 
be  used  directly  as  it  relates  to  the  time  for  precipitation  falling 
at  the  extreme  point  on  the  watershed  to  contribute  to  the  surface 
runoff  at  its  exit.  By  definition,  ground  water  would  not  enter  the 
stream  at  that  distance.  Ground  water  becomes  surface  water  at 
springs,  which  then  become,  essentially,  the  head  of  perennial  flowing 
streams;  or  by  seepage  at  the  banks  of  a  stream  where  a  water  table 
with  hydraulic  gradient  toward  the  stream  is  intercepted  by  the  stream 
banks.  There  can  be  some  sheet  flow  associated  with  ground  water 
emergence,  however,  in  a  cultivated  area  such  sheet  flow  would 
generally  have  been  intercepted  by  drainage  ditches. 

Values  of  time  of  concentration  are  partially  determined  by  the 
difference  in  elevation  and  are  based  on  storm  flows,  which  are  larger 
and  at  much  increased  velocities  than  when  the  streams  are  fed  only  by 
the  ground  water.  Based  on  the  average  ground  water  flow  (Table  A-2), 
an  average  slope  for  the  lower  basin  of  about  .03  (Table  A-l),  an 
assumed  mean  water  depth  of  1/4  foot  taken  is  equal  to  the  hydraulic 
radius,  a  value  of  n  of  0.0B0  for  a  typical  natural  stream  channel, 
and  the  Manning  formula,  an  approximate  average  stream  flow  velocity 
at  the  watershed  exit  of  0.6  meters  per  second  can  be  derived. 

Use  of  this  value  and  distance  to  the  centroid  of  the  watershed 
(Column  (6)  Table  A-l)  provides  a  basis  for  a  rough  approxiinat i on  of 


A-52 


the  average  time  for  ground  water  to  be  on  the  surface  of  the 
watershed  prior  to  exit  (t  ).  The  resulting  expression  for  volume  of 
ground  water  on  the  surface  at  any  instant  (Vg,  in  m  )  is  as  follows: 

V  { m3 )  =  A  (;/sec/km^)  x  Area  (km^)  x  t  (sec)  /  1000 
g  u  ga 

where  A^  is  from  Table  A-2  for  the  corresponding  month 

Area  is  from  Table  A-l,  column  (4) 

t  =  distance  to  center  of  area  (km)  (Table  A-l,  Col.  (6)) 
ga  x  1000/. 6 

The  volume  for  ground  water  on  the  surface  at  any  time  may  be 
added  to  the  volume  of  storm  surface  runoff  determined  from  Figure  A-7 
to  provide  an  estimate  of  the  total  volume  of  surface  water,  flowing 
towards  the  exit,  on  a  watershed  at  a  given  time. 
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SECTION  A- 5 


CONCLUSIONS 

The  above  discussion  and  tabulated  data  provide  bases  for 

determining:  (1)  the  time  following  ground  contamination  that  is 

likely  to  elapse  before  soluble  components  of  the  contaminant  are  apt 
to  arrive  at  a  water  supply  point,  (2)  the  volume  of  water  that  may  be 
expected  to  arrive;  and  (3)  the  sources  of  that  water.  Soluble 

material  that  enters  the  ground  water  would  take  much  longer  to 
arrive,  and  its  arrival  would  be  distributed  over  years.  The  estimated 
distribution  of  its  arrival  has  not  been  estimated  and  application  of 
such  empirical  and  theoretical  bases  as  exist  for  such  calculations 
would  be  highly  speculative. 

In  the  absence  of  more  detailed  and  coordinated  German 

precipitation  stream  flow  records,  and  with  recognition  of  their 
intended  application,  the  developed  expressions  appear  appropriate. 
The  times  of  concentration  could  generally  be  ignored  due  to  their 

relative  brevity,  however,  the  estimates  are  provided  (Table  A-l)  for 
use  in  any  cases  of  surface  runoff  concurrent  with  or  immediately 
following  creation  of  neutron  induced  radioactivity,  fallout,  or 
rai nout. 
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SECTION  B-l 


INTRODUCTION 


This  appendix  presents  a  technical  description  of  the  computer 
code  WSWCM  (Watershed  Water  Contamination  Model).  WSWCM  has  been 
developed  by  Science  Applications,  Inc.  ( S A I )  under  contract  to  the 
Defense  Nuclear  Agency  (DNA)  for  use  on  the  project  entitled  "Nuclear 
Warfare  Water  Contamination  Threat  Assessment". 

The  purpose  of  the  computer  code  WSWCM  is  to  determine  the 
radiological  water  contamination  that  would  occur  due  to  the 
radioactive  fallout  from  a  nuclear  weapon  detonation.  The  focus  is  on 
the  water  contamination  threat  to  U.  S.  Army  field  forces  that  would 
arise  in  the  event  of  nuclear  warfare  in  Europe. 

A  review  of  the  literature  has  shown  that  many  water 
contamination  models  have  been  developed,  and  are  being  developed,  to 
address  specific  situations  such  as  the  potential  water  contamination 
associated  with  chemical  and  nuclear  waste  storage  facilities;  land 
use  changes  like  industrial ization  and  urbanization;  and  the 
application  of  herbicides,  pesticides,  and  fertilizers  in  agricultural 
areas. Typically,  these  models  are  very  complex; 
require  a  considerable  amount  of  detailed  information  on  the 
character i st ics  of  the  contamination  source,  the  relevant  chemical  and 
physical  processes,  and  the  water  source;  and  are  often  "calibrated" 
so  as  to  reproduce  the  results  of  actual  field  data.  The  application 
of  such  sophisticated  models  for  this  water  contamination  threat 
assessment  is  not  appropriate  nor  necessary  because  of  the  scoping 
nature  of  the  assessment  and  the  lack  of  detailed  input  information. 


*The  number  in  the  parentheses  denotes  a  reference  that  is  identified 
in  Section  B-5. 
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For  the  purposes  of  this  assessment,  a  simple  model  which 
incorporated  the  major  factors  that  affect  water  contamination  and 
required  a  minimum  of  input  data  was  deemed  sufficient.  Accordingly, 
it  should  be  understood  that  WSWCM  is  a  very  simple  water 
contamination  model  intended  only  for  scoping-type  calculations. 

The  technical  approach  used  in  WSWCM  is  discussed  in  Section  B-2. 
A  discussion  of  the  computer  programming  for  WSWCM,  including  input 
and  output  descriptions,  is  given  in  Section  B-3.  Section  B-4 
provides  a  sample  problem  that  illustrates  the  use  of  WSWCM. 
Referenced  material  is  identified  in  Section  B-5. 
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SECTION  B-2 


TECHNICAL  APPROACH 


B-2.1  Overview 


The  Watershed  Water  Contamination  Model  (WSWCM)  calculates  the 
time-dependent  activity  concentration  of  fission  product  radionuclides 
dissolved  in  water  that  could  result  from  the  deposition  of  nuclear 
weapons'  fallout  on  a  watershed.  WSWCM  considers  both  the  prompt 
water  contami  nat  i  on  that,  would  result  from  the  fallout  material 
deposited  directly  in  the  water  and  the  delayed  water  contamination 
that  would  result  from  the  fallout  material  initially  deposited  on  the 
land  surface  and  subsequently  transported  to  the  water  by 
precipitation  runoff.  All  activity  is  assumed  initially  to  be 
associated  with  solid  particulate  fallout.  The  activity  may  leave  the 
watershed  only  by  radioactive  decay  or  by  being  dissolved  in  water 
which  flows  past  the  water  supply  point. 

The  principal  characteri s t i cs  of  WSWCM  are:  (1)  tne  watershed  is 
modeled  as  a  mix''  .  tank,  (2)  radionucl  i de-spec i f  ic  distribution 
coefficients  are  used  to  address  fallout  solubility,  and  (3)  the 
model  treats  radioactive  decay  including  daughter  in-growth.  It  is 
important  to  note  that  WSWCM  addresses  radioactive  material  in 
solution  but  does  not  incorporate  any  modeling  of  particulate  or 
sediment  transport. 

The  technical  approach  adopted  for  WSWCM  is  described  in  the 
following  material.  Section  B-2. 2  discusses  the  stream  water 
con' ami nat i on  model  that  treats  the  prompt  water  contamination.  The 
runoff  water  contamination  model  that  treats  the  delayed  water 
contamination  is  discussed  in  Section  B-2. 3.  A  discussion  of  the 
i r 1 1  j ’ .  data  required  for  the  models  is  given  in  Section  B-2. 4. 
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B -2 . 2  St  redin  Water  Contami  net  i  on  Model 


The  stream  water  contamination  model  addresses  the  radiological 
contamination  that  results  from  the  dir-wi  de;  is  1 1 1  on  of  fallout 
material  into  the  stream.  The  model  Jet  ermines  the  time-dependent 
concentrations  of  rad  i  onuc  I  i  des  m  •  hi  s' red  w.Prr  that  passes  the 

water  supply  point. 

The  stream  water  com  arinatior  node  1  uses  four  separate 

compartments  to  represent  i  ne  arount  of  parent  and  daughter 

radionuclides  in  solid  and  1  iquid  ,  liases.  The  mechanisms  for  transfer 
from  one  compartment  to  another  compartment  include:  radioactive 
decay,  dissolution  from  the  solid  phase  to  the  liquid  phase,  and 
transport  of  liqu’d  phase  contaminants  out  of  the  system  by  water 
out-flow.  The  concentration  of  a  radionuclide  in  the  stream  water  is 
determined  by  dividing  the  activity  resent  in  the  liquid  phase  by  the 
volume  of  water  in  the  stream. 

The  basic  equations  of  the  stream  water  contamination  model  are 
described  in  the  following  material.  Section  B-2.2.1  discusses  the 

f our-compartment  model  and  the  water  concentrat i on  calculation.  The 
modeling  of  the  dissolution  from  the  solid  phase  to  the  liquid  phase 
is  discussed  in  Section  B -2.2.2.  Section  B-2.2.3  addresses  the  mixing 
tank  model  that  is  used  to  treat  the  liquid  phase  transport  of 
contaminants  otd  of  the  system. 

B-2.2.1  Compartment  Model  1 ng  and  Concentration  Calculation 

The  f our-compartment  model  used  to  keep  track  of  the  amount  of 
parent  and  daughter  radionuclides  in  the  solid  and  liquid  phases  is 
shown  in  Figure  B-l.  For  each  compartment,  a  differential  equation 
can  be  written  based  on  the  rates  at  which  material  enters  and  leaves 
the  compartment;  the  differential  equation  can  then  be  solved  to 
obtain  the  amount  of  material  in  t.he  compartment  as  a  function  of 
time.  These  equations  and  their  solutions  are  shown  in  Figure  B-2. 
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Figure  B-l.  Four-compartment  model. 
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Fiqure  B-2.  Compartment  equations  and  solutions. 
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The  model  and  equations  given  in  Figures  B-l  and  B-2  are  used  to 
determine  the  amount  of  a  radionuclide  present  in  the  liquid  phase; 
for  example,  Qw-(t)  for  radionuclide  i.  The  activity  concentration  of 
the  radionuclide,  C  ,(t),  is  determined  by  its  radioactive  decay 
constant,  .  .,  and  the  volume  of  water  in  the  stream,  V,  by 

M  (t) 

Cwi  ^  V 

B-2.2.2  Dissolution  of  Radionuclides* 


At  equilibrium,  the  distribution  of  a  radionuclide  between  the 
solid  phase  and  the  liquid  phase  is  expressed  by  a 
radionuclide-specific  distribution  coefficient,  Kd.  By  definition, 


Kd  =  amount  of  radionuclide  sorbed  onsolid  phase, 
amount  of  radionuclide  left  in  solution 


so 


where 


Kd 


Q|q  /m 
Qeq  TV 


0^  -  amount  of  the  radionuclide  present  in  the  solid  phase 

at  equi 1 i bri urn, 

Q  -  Amount  of  the  radionuclide  present  in  the  liquid  phase 
at  equilibrium, 

m  =  mass  of  solid  phase  material,  and 


V  =  volume  of  liquid  phase  material. 


*In  this  discussion,  the  subscript  i,  denoting  a  specific  radionuclide, 
has  been  intentionally  omitted  to  avoid  unduly  complicating  the 
equations. 
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The  dissolution  of  a  radionuclide  from  the  solid  phase  to  the 
liquid  phase  is  modeled  as  a  first  order  reaction  with  a  constant 
coefficient , 


so 


''■•elide  ; s  initially  present  only  in  the  solid 


amount  of  the  radionuclide  present  in  both 


-1  t 


‘  W 


I  <\ 

'  '  is  assum'd  that  the  equilibrium  between  the  solid 

.  achieved  wthin  a  time  corresponding  to  the  time 
;  *  '  ;  *  a n k  model  ,  .  So 
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thus 


1  in  [1+  V 


mKd 


B-2.2.3  Mixing  Tank  Model 

The  portion  of  the  stream  that  is  up-stream  from  the  water  supply 
point  is  considered  to  be  a  mixing  tank  in  which  i nstantaneous  and 
uniform  mixing  of  the  fallout  material  and  the  stream  water  occurs.  A 
time  character! st ic  for  this  mixing  tank  model  is  obtained  by  dividing 
the  volume  of  water  in  the  stream  by  the  stream  flow  rate. 

With  a  mixing  tank  model,  the  amount  of  a  contaminant  in  the 
stream  is  given  by 


X(t) 


Xoe 


-t/T 


where 


Xq  is  the  initial  amount  of  the  contaminant,  and 
T  is  the  time  characteristic  of  the  model. 


B - 2 . 3  Runoff  Water  Contaminat ion  Model 

The  runoff  water  contaminat ion  model  addresses  the  radiological 
contamination  that  results  from  the  fallout  material  that  is  initially 
deposited  on  the  land  surface  and  is  subsequently  transported  to  the 
water  by  precipitation  runoff.  The  model  determines  the 
time-dependent  concentrat i ons  of  radi onucl i des  in  the  runoff  water  and 
couples  the  watershed  runoff  with  the  stream  water  flow  to  give  the 
contamination  of  the  water  that  passes  the  water  supply  point. 


The  runoff  water  contami nat i on  model  follows  the  activity  of 
parent  and  daughter  redionucl ides  during  the  time  after  initial 
surface  deposition  when  they  are  subjected  to  successive  rains.  The 
radioactive  material  in  the  solid  phase  is  affected  by  radioactive 
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decay  and  dissolution  into  the  liquid  phase.  The  radioactive  material 
in  the  liquid  phase  is  also  affected  by  radioactive  decay  and 
transport  off  the  land  into  the  stream  by  precipitation  runoff. 

The  basic  equations  of  the  runoff  water  contamination  model  are 
described  in  the  following  material.  Section  B-2.3.1  discusses  the 
method  used  to  determine  the  amount  of  parent  and  daughter 
radionuclides  in  each  phase  as  time  progresses.  The  runoff  water  flow 
model  is  addressed  in  Section  B-2.3.2.  The  equations  used  to  determine 
the  stream  water  contamination  that  results  from  the  contaminated 
runoff  water  are  discussed  in  Section  B-2.3.3. 

B-2.3.1  Contamination  Time  Sequence  Modeling 

The  time  sequence  modeling  used  for  the  runoff  water 
contamination  model  assumes  that  the  fallout  material  is  deposited  on 
the  land  surfae  at  time  tQ  and  rains  subsequently  occur  at  times  t^, 
t^,  . ..t  .  Each  rain  causes  a  partioning  of  the  radioactive  material 
between  the  solid  phase  and  the  liquid  phase.  The  material  in  the 
solid  phase  remains  on  the  land  surface  and  is  subjected  to  further 
phase  partioning  by  successive  rains.  The  material  in  the  liquid  phase 
is  transported  off  the  land  surface  and  into  the  stream  by  the 
precipitation  runoff. 

Figure  B-3  shows  the  sequence  of  events  for  the  radioactive 
material,  both  parent  and  daughter  radionuclides,  present  in  the  solid 
phase.  Each  of  these  events  is  described  below: 

(a)  t  =  t„ .  The  initial  amounts  of  parent  and  daughter  radionuclides 
present  in  the  solid  phase  are  given  by  Q°^^A  and  Q0^ ^ ^ ,  respectively. 

( b )  to  <  t  <  1 1 .  During  the  time  period  prior  to  the  first  rain,  the 
process  of  radioactive  decay  affects  the  amount  of  the  parent  and 
daughter  radionuclides  present  in  the  solid  phase.  So, 

QS1AU)  =  Qs1Ae“V  and 

-  °  I1]  *  °V'b\ 
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where  and  'p  are  the  radioactive  decay  constants  for  the  parent  and 
daughter  radionuclides,  respectively. 

( c )  t  =  tj.  the  first  rain  occurs  at  time  tj.  It  is  assumed  that  the 
radionuclides  instantaneously  reach  an  equilibrium  distribution 
between  the  solid  phase  and  the  liquid  phase. 


As  pointed  out  in  Section  B-2.2.2,  the  distribution  coefficient  Kd  is 
given  by  the  expression* 


so 


Kd  - 


Qs/m 

cT/V 

w 


but  the  total  amount  of  the  radionuclide,  '■} .  ,  i  s 


Ur  - 

I 

so  UT  - 

or  U$  ■ 

with,  of  course. 


Q  +  '} 

‘ S  'W 


This  dissolution  of  ,1  radionuclide  from  the  solid  phase  to  the  liquid 
phase  results  in  the  following  equations: 


♦Note  that  different  values  of  m  and  V  are  used  for  the  stream  water 
contain! nut i on  model  and  the  runoff  water  contami nat 1  on  model. 
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QS2A(V  r  °S1A( fcl5 '  t1  +  mKd,  ] 

M 


and 

QS2B(  ll )  =  QSlB^tl)  '  l 1  +  -Fl—  ] 

mKdn 

b 

where  specific  distribution  coefficients  are  used  for  the  parent  and 
daughter  radionuclides  (Kd^  and  Kdg,  respectively),  and  indicates 
the  volume  of  potential  runoff  wafer  associated  with  the  first  rain. 


(d)  t  ]_  s  t  s  to  •  during  the  fire  period  between  the  first  and  the 
second  ram,  the  process  of  r.il inactive  decay  affects  the  amount  of 
the  , a  rent  and  daughter  radun  ic lidos  present  in  the  solid  phase.  So, 


bch  bi.t\  1 


jnd 


Qs?i,(t)  -  'WV* 


*  0,„.(1. e  B 


A 

B  •;A 
(t-t, ) 


'■A(t_tr  -  e"  l:.^  t_tl }] 


( e )  t  -  t  p  .  The  'ream:  rain  occurs  at  time  t„.  ’his  rair  causes 
another  part :  »r>  i  »'g  of  he  radioactive  material  mm  the  solid  and 
liquid  phases. 


So, 


a  no 


l)S3B^t?)  '  QS2B-  4' 


l> 

[> 


i' 


-1 


-1 


This  model in-i  of  decay  and  dissolution  can  be  continued  as  long 
as  the  radioac*.  l  v«*  Mat  erial  ,  resent  in  the  solid  phase  serves  as  a 
source  for  the  radiological  .  -nt  .i.ainj!  ion  e*  the  water. 
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As  was  seen  in  Figure  B-3,  radioactive  material  enters  the  liquid 
phase  at  the  time  of  each  rain,  (i.e.,  tj,  t^,  etc.)-  Each  rain  is 
treated  as  in  individual  event  and  the  amount  of  radioactive  material 
present  in  the  liquid  phase  is  affected  only  by  radioactive  decay. 


For  the  first  rain,  the  amount  of  radioactive  material  initially 
present  in  the  liquid  phase  is  given  by 


and 


QW1A( tl)  =  QS1A(V  ‘  j1  "  [l  +  niK JA  ] 
QWlB(tl]  =  °SlB(tl)  '  j1  -  [l  +  jj^-]  ]  • 


At  any  subsequent  time,  the  amount  of  parent  and  daughter 
radi  onucl  ides  present  in  the  liquid  phase,  as  a  result  of  the  first 
rain,  is  given  by 


and 


^WIB^  =  ^WlA^ 


(t-t,)  _ 


B'  A. 


""  QW2A^t?^  ’  0  ^  2 ^ 


n  for 

the  second 

rain 

>- 

i,  *  .Iq 

I" 

1 

1  j 

{.- 

h+  !i.] 

|- 

1 

*•  mKdgJ 

1  1 

QW2B(t)  r-  Qw2A^2^ 


■fe)f 


and 

e"‘-Alt‘L2l  -  e_'& 


->n(t-t9)  „  0-'R(t-t2)l 


+  ^B^  V  ’  0  B 


_ ‘ r( E' tp) • 
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B-2.3.2  Runoff  Water  flow  Model 


Figure  B-4  shows  the  time  history  of  runoff  water  on  the 

watershed.  The  figure  is  based  on  an  engineering  estimate  of  the 

hydrological  character i st ics  of  a  typical  watershed.  The  key  features 

of  the  time  history  are:  (1)  the  volume  of  water  on  the  surface 

reaches  a  maximum  at  8  hours,  (2)  the  surface  runoff  is  completed  in 

4  days,  and  (3)  the  integral  of  the  time  history  curve  is  1000 
3  2 

m  /Km  per  mm  of  SRO  (surface  runoff). 

An  approximation  to  the  time  history  curve  of  Figure  B-4  can  be 
developed  by  considering  a  mathematical  model  for  the  volume  of 

_  1  f 

surface  water,  V(t),  that  has  an  input  rate  of  RQe  and  an  output 
rate  of  kV. 

Mathematical ly, 

4V-  =  +  R  e"' t  -  kV 

dt  o 

with  V ( o )  =  o 

thus  V(t)  =  Rq.  Ri_  •  [e"  lt  -  e~kt  j 

The  values  of  the  parameters  of  the  equation  ( i .  e . ,  R  ,  k,  and  \)  can 
be  determined  by  trial  and  error  with  the  objective  of  matching  the 
three  key  features  of  the  time-history  curve  shown  in  Figure  B-4.  The 
parameter  values  of  Rq  -  313./  m  /hr  per  mm  of  SRO,  k  =  0.0348  1/hr, 
and  >  -  0.3011  1/hr  provide  a  suitable  fit,  as  shown  in  Figure  B-4. 

Using  this  model  approximation,  the  rate  at  which  the  runoff 
water  flows  off  the  watershed  given  by 

r  =  kv 

f(t)  -  Rq  •  k-k,  •[e"'t-o"kt] 
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Engineering  Estimate 


Time  After  Storm  (Hour) 

1  nine  of  storm  runoff  on  watershed. 


where 


3 

F  (t)  -  flow  rate  of  runoff  water,  m  /hr  per  mu  of  SRI), 

Ro  =  313.7  m'Vhr  per  nun  of  SRO  (a  fitted  constant), 

k  -  0.0348  1/hr  (a  fitted  constant),  end 

\  1  0.3011  1/hr  (a  fitted  constant). 


A  plot  of  the  runoff  water  flow  rate  as  a  function  of  time  is 
shown  in  Figure  B-5. 

B-2.3.3  Water  Contamination  Calculation 


The  equations  presented  in  Section  B-2.3.1  determine  the  amount 
of  parent  and  daughter  radionuclides  present  in  the  liquid  phase  as  a 
function  of  time  for  a  specific  rain.  The  activity  concentration  of 
the  radionuclides  in  water  is  determined  by  the  radionuclide 
radioactive  decay  constant  and  the  volume  of  potential  runoff  water 
associated  with  the  rain.  Tor  example,  the  i^'  rain, 


and 


wm  •  :a 


c  (t)  =  i  a  aril0 

LwiBlti  V. 

i 


The  material  presented  in  Section  B-2.3.2  described  a  runoff 
water  flow  model  and  presented  an  equation  for  the  time-dependent  flow 
rate  for  the  i'  rain, fa  ( t - 1 ) . 


When  the  contami nat.ed  runoff  water  enters  the  stream,  it  mixes 
with  the  normal  stream  flow,  T,  and  the  existing  stream  water 
contamination  C-w^(t)  and  Cw^(t),  as  determined  by  the  model  described 
in  Section  B-2.2 


The  resulting  water  contaminat ion  is  given  by 

n 


r  •  cwA(t)  *  £  rpt-tp  ■  cwjfl(t) 


cat) 


F,X, 
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whore  the  L  indicates  the  inclusion  of  up  to  n  rains. 
i  =  l 

B - 2 . 4  Input  Data  tor  Mode  1 s 

The  streuiii  water  contamination  model  and  the  runoff  water 
contain! nat i on  model  require  three  basic  types  of  input  data  or 
information:  (1)  fission  product  radionuclide  data,  (2)  information 
on  watershed  and  precipitation  characteristics,  and  (3)  data  and 
information  for  solubility  modeling.  The  following  material  discusses 
the  approach  used  to  provide  the  necessary  data  and  information  for 
each  of  the  three  categories. 

B-2.4.1  f_i  s_s  i  on  :'roducf  Radi  onucl  ides 

The  radioactive  fallout  from  a  nuclear  weapon  explosion  contains 
several  hundred  fission  product  radionuclides.  For  a  variety  of 
reasons  (e.y.,  relative  yield,  radioactive  decay  characteri st i cs , 
water  solubility,  radiation  dosimetry,  etc.),  not  all  of  the  fission 
product  radionuclides  are  of  significance  to  the  water  contamination 
threat  analysis.  Table  B-l  lists  those  fission  product  radionuclides 
that  have  been  identified  as  of  importance  to  the  threat  analysis.* 

In  Table  B-l,  the  radionuclides  are  given  as  pa  rent -da ughter 
pairs  (e.g.,  Sr-90,  and  Y-90),  where  relevant.  This  is  an  important 
aspect  of  rad,  'active  contami nat i on  modeling  since  there  are  cases 

*This  identification  was  bast'd  on  an  assessment,  of  the  relative 
importance  of  specific  radionuclides  present  in  unf ract i onated 
fission  products  in  terms  of  ingestion  dose  commitments  using 
data  from  Reference  B-S. 


Table  B-l.  Fission  product  radionucl ides . 


Parent 
ladionucl ide 

Normalized 

Ground  Concentration* 
(Ci/km2  for  1  R/hr  at  H+lJ^ 

Daughter 

Radionuclide 

Normalized 

Ground  Concentration* 
(Ci/km2  for  1  R/hr  at  H+l) 

Sr- 89 

4.1 

Sr-90 

.027 

Y-90 

.0060 

Sr-91 

660. 

Y-91 

.21 

Zr-95 

4.8 

Nb-95 

.0030 

Zr-97 

410. 

Nb-97 

220. 

Mo -99 

100. 

Tc-99m 

0.0 

Ru- 103 

4.2 

Rh-103 

0.0 

Ru-105 

270. 

Rh' 105 

5.7 

Ru- 106 

.079 

Rh-106 

0.0 

Sb- 127 

2.2 

Te- 12  7 

0.0 

Te-  129m 

.070 

Te- 129 

0.0 

Te-131m 

13. 

1-131 

15. 

Te- 132 

71. 

1-132 

30. 

1-133 

300. 

1-134 

4400. 

1-135 

990. 

Ba-140 

24. 

La-140 

.63 

Ce- 141 

.93 

Ce- 143 

200. 

Pr-i43 

.30 

Ce-144 

.94 

Pr-144 

0.0 

Cs-137 

.28 

ild- 147 

11. 

Pm- 147 

0.0 

*The  normalized  ground  concentrations  were  obtained  using  the  SAI  computer 
code  F 1 1  DOS  described  in  Reference  B-6.  The  calculated  ratic  of  exposure 
rate  to  the  surface  contamination  was  13.4  R/hr  per  Ci/m2  at  H+l  hour. 
This  ratio  is  equivalent  to  2400  R/hr  per  KT/mi2,  which  compares  quite 
well  with  the  value  of  the  K-factor  reported  in  Reference  B-7. 


! 


| 
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where  the  initial  .amount  of  a  daughter  radionuclide  is  small  but  large 
amounts  are  subsequently  added  by  the  radioactive  decay  of  the  parent 
radionucl ide. 

2 

Table  B-l  also  shows  the  normalized  ground  concentration  (Ci/Km 
for  1  R/Hr  at  H  +  l  Hour)  for  each  radionuclide.  These  concentrations 
are  based  on  the  calculated  inventory  of  unfractionated,  U-?35  fission 
products  present  at  one  hour  after  the  fission  event  occurs.  The 
concentrations  are  referenced  to  a  fallout  deposition  contour  that  has 
an  above  ground  external  radiation  exposure  rate  of  1  R/Hr  at  H+l  hour. 

B - 2 . 4 . 2  Watershed  and  Precip itation  Characteri sties 

The  watershed  and  precipitation  characteristics  used  in  WSWCM 
pertain  to  possible  field  water  supply  points  lucated  in  the 
preselected  scenario  area.  This  area,  shown  in  figure  B-h,  is  bounded 
by  Marburg,  Giessen,  Frankfurt  am  Main  on  the  wos* ,  and  the  Fulda 
River  valley  on  the  east.  The  descriptive  information  on  the 
watershed  and  precipitation  character! st ics  for  the  scenario  area  is 
contained  in  Appendix  A,  “Water  Source  Information." 

Tables  B-2  and  B-3  provide  information  on  the  watersheds  that 
support  specific  potential  water  supply  points  sited  in  the  scenario 
area.  This  information  is  used  by  WSWCM  to  determine  the  mixing  tank 
time  characteri st ic ,  :  ,  and  the  volume  of  water  in  the  stream,  V,  for 
the  stream  water  contamination  model. 

The  time  character i st ic  (hr),  is  determined  by  the  distance  from 
the  center  of  the  watershed  to  the  water  supply  point,  D  (Km),  and  the 
estimated  average  stream  flow  velocity,  S  (m/s),  by 

P  _1_ 

S  '  3.6 
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Scenario  area 


Table  B-2.  Selected  data  on  water  supply  points. 


Water  Supply 

Point  No. 

Area 

(Km-) 

Distance  from  Center 
of  Area  to  WSP  (Km) 

1 

63  6 

5  8 

2 

10  0 

1  4 

3 

12  6 

3.2 

4 

27  1 

2  7 

5 

12  8 

4.0 

6 

28  1 

3  5 

7 

12  2 

3  6 

8 

30.9 

2  0 

9 

13.6 

2  5 

10 

7  9 

3  0 

11 

12.4 

4  5 

12 

16.6 

5.3 

13 

9  3 

3  2 

14 

21  4 

3  9 

15 

33  0 

7  2 

16 

17  4 

3.2 

17 

10  8 

2.3 

18 

8.4 

2  6 

19 

20,5 

3.9 

20 

22  1 

3  3 

21 

8  7 

3.1 

22 

14  3 

3  1 

23 

19  9 

4  0 

24 

31  7 

4  6 

25 

45  3 

6  4 

26 

15  4 

2  6 

27 

9  5 

2  0 

28 

12  5 

2  7 

29 

25  8 

4  2 

30 

24  7 

4  0 

31 

40  2 

4.2 
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Table  B-3.  Ground  water  areal  flow  rate. 


Ground  Water  Runoff 


Month 

( ' /s/Km- ) 

January 

3  7 

February 

5  0 

March 

4  9 

Apri  1 

5  0 

May 

3.0 

June 

2  3 

July 

1  9 

August 

1  9 

September 

1  9 

October 

2.2 

November 

3.1 

December 

3.4 
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where  the  factor  1/3.6  converts  Km/in/s  to  hr.  Values  for  D  can  be 
found  in  Table  B-2  for  each  water  supply  point  watershed.  The  average 
stream  velocity  at  the  watershed  exit  has  been  estimated  to  be  0.6 
meters  per  second. 

The  volume  of  water  in  the  stream  V  (r).  in  the  absence  of 
prec i pi t at  i  on  runoff,  is  determined  by  the  ground  wafer  areal  flow 
rate,  i  (•  /shm"),  t  he  area  of  the  watershed,  A  (K;:;‘"),  and  the 
watershed  fie  character i st i c,  (hr),  by 

V  =  f  •  A  .  -  .  (3600) 
whore  the  factor  3600  converts  s  to  hr. 

Tables  B-4,  B-f>,  and  B-o  provuie  information  on  the  precipitation 
characteristics  of  the  scenario  area.  This,  information  is  used  by 
WSWCM  to  determine  the  volume  of  surface  runoff  water  caused  by 
precipitation  for  the  runoff  wafer  cent  ami  nation  model . 

The  volume  of  s  ;>'face  runoff  water  resulting  from  given  rain.  V 
(  ■),  is  determined  by  the  area  of  the  watershed,  A  (Km'' ) ,  and  the 
linear  amount,  of  surface  runoff  wafer,  SRU  (mm),  by 

V  -  A  •  SR0  .  in6 

where  the  factor  in  converts  km'-mm  *  o  .  . 

The  linear  amount  of  surface  runoff  water  resulting  fro::  a  given 
rain  is  determined  by  the  amount  of  water  depus  ted  by  the  rain,  the 
prior  precipitation  history  of  the  watershed,  and  the  hydrological 
characteristics  of  watershed.  The  equations  used  for  calculating  the 
linear  amount  of  surface  runoff  are: 

SRI!  [ SR  )/P  J  x  I' 

e J  e 

wi  th 
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Table  B- 


4.  Estimated  proportion  of  excess  rainfall 
which  becomes  surface  runoff. 


Month 

Surface  Runoff/E) 

January 

0.45 

February 

0.65 

March 

0.58 

April 

0.34 

May 

0.11 

June 

0  10 

July 

0.13 

August 

0.10 

September 

0.12 

October 

0.16 

November 

0.40 

December 

0.42 
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Table  B-5.  Daily  precipitation  data  (tnm)  -  Fulda  (  1975) 


n'.  *' 

[■'. . 

Jl.M 

Jiil 

A  .  : 

Si.-:' 

Oc  *. 

t 

t 

t 

G.i 

*5  ^ 

L.  m  / 

O.’l 

2.3 

o.: 

o.a 

0.1 

1.3 

4.-'. 

U  - 

1.4 

0.3 

1.0 

t 

Of 

3.3 

t 

9.7 

0.9 

t 

t 

6.1 

K.S 

t 

1.7 

4.1 

4 . 

iu.: 

t 

t 

t 

3.9 

t 

2.5 

0.9 

1.2 

4.2 

6.9 

t 

1.1 

2.3 

0.5 

9.6 

6.4 

3.  ; 

7.5 

11.2 

3 

35  7 

1.2 

1.2 

: . : 

t 

4.3 

t 

0,1 

O.o 

1.7 

t 

1.3 

t 

2 , 7 

n  • 

r 

4.1 

0.  1 

4 . 

0.2 

r  ■ 

3.6 

0.' 

3 .  C 

3.3 

0.2 

0.3 

c  : 

3.r» 

9.4 

If ..." 

0.4 

0.6 

0.2 

07.5 

If  .•■ 

t 

t 

0.1 

0." 

t 

14.3 

o.r 

0.3 

0.2 

4.3 

1  • 
c  .  : 

t 

15.  ' 

t 

O.f 

1  } 

0.1 

0.  ' 

0.7 

t 

2 .  ‘ 

0.  • 

n.  ’ 

5. 

t 
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Table  B-6.  Cummulative  preci pi ta tion  necessary 
for  added  surface  runoff. 


Month 


Threshold  (nun) 


January 

8 

February 

4 

March 

5 

Apri  1 

9 

May 

12 

June 

12 

July 

13 

Augus  t 

13 

September 

12 

October 

12 

November 

10 

December 

9 
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and 


API  Y. 


p  x  m.9) 
o- n  ■  ' 


where 


SRO  -  surface  runoff  (mm) 


[SRU/i'J  surface  runoff-excess  precipitation  ratio, 
as  given  in  Table  B-4. 


>\  -  excess  precipitation  (mm) 

24-hour  precipitation  (mm);  given  in  Table  B-5 

AIM  -  antecedent  precipitation  index  (mm) 

CRT  =  cumulative  precipitation  threshold  (mm),  given 
in  Table  B-(> 


-  24-hour  precipitation  n  days  before  the  day  of  P 

(but  if  Po-n  >  CRT  use  Po-n  -  CPT,  to  avoid  double 
counting  of  prior  precipitation  that  resulted  in 
surface  runoff),  as  given  in  Table  B-5. 


(3-2.4. 3  Sol  ubi  1  i  ty  Model  i ng 


Very  little  information  is  available  on  the  solubility  of 
specific  radionuclides  present  in  nuclear  weapons  fallout.  Information 
is  also  loci  ini)  on  the  rate  at  which  the  fallout,  dissolves  and  the 
radionuclides  enter  the  liquid  phase.  Most  of  the  statements  found  in 
the  available  literature  regarding  fallout  solubility  refer  to  fallout 
as  basically  insoluble  but  cite  some  radionuclides  present  in  the 
fallout  as  soluble. 


To  handle  solubility  modeling  in  WSWCM  use  is  made  of 

radionuclide-specific  distribution  coefficients.  Information  on  these 

distribution  coefficients  is  found  in  those  portions  of  the  nuclear 

power  literature  dealing  with  environmental  contami nat i on  of  water 

bodies  by  routine  releases  from  nuclear  power  plants  and  with  the 

impacts  of  potential  releases  of  radioactive  material  from  long-term 

(B-p,  B-Q  B-10) 

nuclear  waste  storage  facilities.  '  ’  ’  ■  The  use  of  these 

distribution  coefficients  for  fallout  solubility  modeling  provides  a 
reasonable  approach  in  the  absence  of  adequate  data  on  actual  fallout. 
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A  distribution  coefficient  is  defined  as: 


amount  of  radionuclide  sorbed  on  solid  phase 
amount  of  radionuclide  left  in  solution 

Since  the  solid  phase  activity  is  usually  expressed  in  units  of  Ci/g 
and  the  liquid  phase  activity  in  units  of  Ci/m.  ,  Kd  typically  has 
units  of  m  /g.  Table  R-7  shows  selected  values  of  distribution 
coefficients  lor  those  elements  whose  radioisotopes  are  considered  in 
WSWCM.  for  a  specific  element,  the  value  of  Kd  is  dependent  upon  the 
chemical  state  of  the  element,  the  type  of  solid  matrix  in  which  it 
exists,  the  physical  charact er i sf ics  of  the  solid  and  liquid  phases, 
and  the  nature  of  the  dissolution  process;  however,  the  actual 
relationship  of  the  value  of  Kd  to  these  is  generally  not  known. 
Values  of  Kd  are  normally  determined  by  laboratory  or  field 
experiments  and,  as  shown  in  Table  B-7,  these  values  exhibit  a  wide 
range. 

It  is  important  to  note  that  the  distribution  coefficients  refer 
*  o  the  phase  dis'nbufion  of  the  radionuclide  at  equilibrium.  In 
KSWC.M,  assumptions  made  alw,*  reaching  equilibrium  conditions  are,  in 
effec*  ,  assumptions  regarding  the  rate  at  which  the  radionuclides 
dissolve.  Wh ;  1  •  t  he;  e  .issu:  get  ons  do  not  appear  unreasonable,  it  is 
not  really  possible  t>  validate  them  because  of  the  absence  of  actual 
fallout  data. 

To  use  distribution  coefficients  in  WSWCM  it  is  necessary  to 

siecify  the  amount  of  Solid  phase  material  with  which  the  radionuclide 

is  associated.  If  nas  been  assumed  that  this  solid  mass  is  the  amount 

of  soil,  or  fallout  material,  represented  by  a  uniform  deposition  over 

the  witershed  ure.i  of  material  with  a  thickness  of  100  microns  and  a 

density  of  l.T  g/cc.  These  assumptions  give  the  solid  phase  material 

5  ? 

a  uniform  areal  density  of  l.S  x  10  Kg/Km".  This  density  is 
multiplied  by  the  area  of  concern  to  determine  the  mass  of  solid  phase 
material,  for  the  stream  water  contamination  model,  the  area  is  1“  of 
the  watershed  area,  for  the  runoff  water  contain i nat ion  model  the  area 
is  ho*  0f  *he  watershed  area. 
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Table  B-7.  Selected  distribution  coefficients*. 


Element 

Ba 

Ce 

Cs 

I 

La 

Mo 

Nb 

Nd 

Pm 

Pr 

Rh 

Ru 

Sb 

Sr 

Tc 

To 

Y 

Zr 


Kd 

fallal 

500 

10,000 

1,000 

10 

500 

25 

10,000 

10,000 

10,000 

10,000 

5,000 

5,000 

100 

1,000 

1 

100 

1,000 

1,000 


*The$e  values  were  selected  from  References  B-8  and  B-9 
values  exhibit  a  wide  range;  for  example  the  value  for 
to  10000  and  the  value  for  Sr  ranges  from  8  to  4000. 


Reported  Kd 
Zr  ranges  from  1000 
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SECTION  B-3 

COMPUTER  PROGRAM 


B-3. I  Introduction 


The  computer  program  WSWCM  is  written  in  “°RTRAN- 1 V  for  operation 
on  the  PDP-11  mini-computer  at  SAI -Schaumburg.  WSWCM  requires  16K 
memory  storage.  A  typical  WSWCM  problem  has  a  running  time  of  5 
minutes.  User  problem  input  is  provided  from  a  keyboard  terminal. 
Code  output  is  directed  to  a  printer  and  a  disk  file  for  subsequent 
data  plotting. 

WSWCM  consists  of  a  main  program  with  seven  subroutines.  Stored 
data  arrays  in  the  code  contain  probl em- i ndependent  data. 
Problem-dependent  data  is  input  by  the  user.  Program  output  is 
printed  and  processed  for  input  to  a  separate  data  plotting  package. 
Figure  B-7  illustrates  the  basic  organization  of  WSWCM. 

The  following  sections  provide  a  more  detailed  description  of 
WSWCM.  The  user  input  is  discussed  in  Section  B-3. 2.  Section  B-3. 3 
describes  the  main  program,  subroutines,  and  stored  data  arrays.  The 
program  output  is  discussed  in  Section  B-3. 4.  A  Fortran  listing  of 
WSWCM  is  provided  in  Section  B-3. 5. 

B-3. 2  User  Input 


The  user-supplied,  problem-dependent  input  to 
program  WSWCM  consists  of  three  types  of  data:  Basic 
precipitation  data  and  radionuclide  data.  Table  B-8 
input  data  required. 


the  computer 
problem  data, 
indicates  the 
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Table  B-8.  User  input  data  required  by  WSWCM. 


________  Basic  Problem  Data  -------- 

1.  Fallout  Exposure  Rate  of  Watershed 
2  Problem  Simulation  Time 

3.  Identification  Number  of  Water  Supply  Point 

4  Month  in  which  Problem  Starts 

5.  Number  of  Rains  during  Problem  Simulation  Time 

6  Number  of  Parent-Daughter  Radionuclide  Pairs 

_______  Precipitation  Data  ------- 

7  Time  at  which  i—  Rain  Occurs 

L_ 

8.  Amount  of  Precipitation  in  i—  Rain 

9.  Antecedent  Precipitation  Index  for  i—  Rain 

hi 

10  Month  in  which  i  —  Rain  Occurs 

_______  Radionuclide  Data  ------- 

11.  Identification  Name  of  j  —  Parent  Radionuclide 

12  Radioactive  Half-Life  of  j—  Parent  Radionuclide 

13  Distribution  Coefficient  of  j—  Darent  Radionuclide 

14.  Normalized  Ground  Concentration  of  j —  Parent  Radionuclide 

15.  Identification  Name  of  j—  Daughter  Radionuclide 

f-  hi 

16.  Radioactive  Half-Life  of  j—  Daughter  Radionuclide 

h  hi 

17.  Distribution  Coefficient  of  j—  Daughter  Radionuclide 

t  h 

18.  Normalized  Ground  Concentration  of  j  —  Daughter  Radionuclide 
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B-3.2.1  Basic  Problem  Data 


The  basic  problem  data  is  9 1 ven  on  a  single  card  image  containing 
6  data  entries.  These  data  entries  correspond  to  the  first  six  items 
1 i sted  i n  Table  B-8. 


(1) 


(2) 


(3) 


(4) 


(3) 


(6) 


Fallout  Lxposure  Rate  of  Watershed 
Name:  R 

Format:  [10.3 

Remarks:  The  fallout  exposure  rate  of  the  watershed 
is  the  normalized  external  exposure  of  the 
fallout  expressed  as  R  per  Hr  at  H+l  Hour. 

This  rate  is  assumed  to  be  uniform  over  the 
watershed  area. 

Problem  Simulation  Time 

Name:  LTMAX 

format:  15 

Remarks:  In  the  model,  time  starts  at  0  and  continues  to 
LTMAX.  LTMAX  is  expressed  in  hours. 

Identification  Number  of  Water  Supply  Point 

Name:  NWSP 


Format:  1? 

Remarks:  The  model  contains  data  (see  Section  B-3.3 

regarding  stored  data  arrays)  for  31  specific 
water  supply  poi nts/watersheds. 

Month  in  winch  Problem  Starts 

Name:  MOS 

format:  12 

Remarks:  The  months  of  the  year  are  numbered  from 
1  (January)  to  12  (December). 

Number  of  Rains  during  Problem  Simulation  Time 

Name:  N 

Format:  12 

Remarks:  None 

Number  of  Parent-Daughter  Radionuclide  Pairs 
Name:  NNUCS 

format:  12 

Remarks:  None 
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B-3.2. 2  Precipitation  Data 


The  prec i pi  tat  i  on  data  consists  of  4  data  entries  (Items  7 
through  10  of  Table  B-8)  for  each  rain;  the  data  for  each  rain  is 
given  on  a  single  card  image.  The  user  must  input  N  precipitation 
data  card  images  (see  Section  B-3.2. 1,  Item  5). 

(7)  Time  at  Which  it*1  Rain  Occurs 

Name:  LT1(I) 

Format:  15 

Remarks:  I.T1(I)  is  expressed  in  hours. 

(8)  Amount  of  Prec  i  pi  tat  i  on  in  i^  Rain 

Name:  XPO 

Format:  Cl  0.3 

Remarks:  XPO  is  expressed  in  mm. 

(9)  Antecedent  Precipitation  Index  for  i^  Rain 

Name:  API 

Format:  F10.3 

Remarks:  API  is  expressed  in  mm. 

(10)  Month  in  which  i^1  Rain  Occurs 
Name:  MOR 

Format:  12 

Remarks:  The  months  of  the  year  are  numbered  from 
1  (January)  to  12  (December). 


B-3.2. 3  Radionuclide  Data 


The  radionuclide  data  is  provided  for  parent -daughter 
radionuclide  pairs  on  two  card  images;  the  first  card  image  is  for  the 
parent  radionuclide,  the  second  card  image  is  for  the  daughter 
radionucl  nie.  Each  card  image  contains  4  data  entries;  these  entries 
correspond  to  items  11  through  14,  or  items  15  through  18,  of  Table 
B-8.  The  user  must  input  NNUCS  pairs  of  radionuclide  data  card  images 
(see  Section  B-3.2.  1,  Item  (>).  If  a  radionuclide  is  not  a  member  of  a 
parent-daughter  pair,  the  user  enters  null  data  for  the  parent  and 
inputs  the  specific  radionuclide  data  for  the  daughter  radionuclide. 
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(11)  Identification  Names  of  jth  Parent  Radionuclide 

Name:  NAMLA 

Format:  A4 

Remarks:  If  the  radionuclide  pair  duesn' t  nave  a 
enter  NON! . 

(12)  Radioactive  Half-Life  of  h  Parent  Rudi'u  .<  I  i  de 

Name:  THAL  f A 

format:  L1U.3 

Remarks:  THALFA.  is  expressed  in  hours.  If  the 

radionuclide  pair  doesn't  hu.e  a  ,  oivnt  enter 

1.00011*00. 

(13)  Distribution  Coefficient  of  .i11'  Parent  Radionuclide 

Name:  XKPA 

F  o  rma  t :  L 1 0 .  3 

Remarks:  XKDA  is  expressed  in  units  of  /kg.  If  the 

radionuclide  pair  doesn't  have  a  parent,  enter 
1 . OOOL  +00. 

(14)  Normalized  Ground  Concentrat i on  or  j'1'  Parent  Radionuclide 

Name:  GA 

Format:  '.10.3 

Remarks:  GA  is  expressed  in  .amts  of  Ci  per  R/Hr  at  d+1 

hour.  If  the  radionuclide  pair  doesn't  have  a 
parent  enter  U. 0001  • 00. 

(15)  1  dent  i  f  i  cat  i  on  Name  of  Diughfer  Radionuclide 

Name:  NAMFB 

format:  A4 

Remarks:  None 

(lh)  Radioactive  Half-Life  of  j*' 1  Daughter  Radionuclide 
Name:  THALFB 

format:  1.10.3 

Remarks:  THALFB  is  expressed  in  hours. 

(17)  Distribution  Coefficient  of  jl  Daughter  Radionuclide 
Name:  XKDB 

Format:  LI  0.3 

Remarks:  XKDB  is  expressed  in  units  of  ’'/Kg* 
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(18)  Normalized  Ground  Concentration  of  Daughter  RAD  I ONUC  L IDE 
Name:  GB 

Format:  CIO. 3 

Remarks:  GB  is  expressed  in  units  of  Ci  per  R/Hr  at 
H+l  hou>'. 

B - 3 . 3  Computer  Program 


The  computer  program  WSWCV  consists  of  a  main  program  and  seven 
subroutines.  WSWCM  also  has  stored  data  arrays  that  contain 
problem- i ndependent  data . 

B-3.3.1  Main  Program 

The  main  program  controls  ! he  operation  of  the  model  and  performs 
some  of  tv a  problem  calculations  and  problem  i nput/output  operations. 
The  basic  flow  of  the  '•■.tin  program  is  illustrated  in  Figure  B-8. 

As  seen  in  Ci  cure  R-  ' ,  'he  mam  program  initially  assembles  the 
basic  problem  data  and  precipitat ion  data  provided  as  user  input,  and 
some  of  the  watershed  and  |  recipi  t  at  ion  data  provided  by  the  stored 
data  array;  the  data  are  used  'n  calculate  some  of  the  watershed  and 
preci( itaticn  parameters.  The  main  program  then  begins  to  calculate 
water  contamination  values  by  processing  each  parent-daughter 
radionuclide  pair  in  succession.  Tor  each  parent-daughter  pair,  the 
main  program  reads  the  user-provided  radionuclide  data,  calculates  the 
initial  water  contamination  and  delayed  water  contamination,  and 
outputs  results.  The  water  contamination  calculations  are  performed 
using  a  time  step  of  one  hour  for  the  tirst  hours  of  the  problem 
and  a  time  step  of  3  hours  lor  the  remainder  of  the  problem  simulation 
time.  The  calculated  results  include  the  radionuclide  concentration 
and  the  time  integral  of  the  radionuclide  concentration,  for  both  the 
parent  and  th<-  daughter  radionuclide,  and  a  composite  radionuclide 
concentration  which  aggregates  all  of  the  radionuclides. 
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Each  parent-daughter  radionuclide  pair  is  processed  through  the 
model  for  the  problem  simulation  time.  After  all  radionuclide  pairs 
have  been  processed,  the  problem  is  completed  and  the  computer  program 
terminates. 

B-3.3.2  Subrouti nes 

The  main  program  flow  diagram  presented  in  Figure  B-8  showed  the 
calls  made  to  the  subroutines  in  WSWCM.  Table  B-9  lists  the  seven 
subroutines  and  identifies  their  purposes. 

B-3.3.3  Stored  Data  Arrays 

WSWCM  has  five  stored  data  arrays  that  contain 
probl em- i ndependent  data.  Table  B- 10  lists  the  arrays  and  identifies 
the  data  contained  in  them. 


Program  Output 


The  output  provided  by  the  computer  program  WSWCM  includes: 

o  the  radionuclide  concentrations  in  water  as  a 
function  of  time  for  each  parent-daughter  radio- 
m.cl ide  pair, 

o  the  radionuclide  concentration  in  water  as  a 

function  of  time  for  all  radionuclides  present,  and 

o  the  t ime- i ntegrated  radionuclide  concentrat i ons  in 
water  as  a  function  of  time  for  each  radionuclide. 

This  output  is  directed  to  a  disk  file  (Unit-2)  to  provide  results  for 
subsequent  processing  to  a  graphical  form.* 


♦The  output  could  also  be  directed  to  a  printer  (Unit-3)  to  provide 
results  in  a  tabular  form.  FORTRAN  Write  statements  for  this  purpose 
are  contained  in  the  WSWCM  program,  but  have  been  Commented  out  since 
they  are  not  needed  if  the  disk  file  is  used. 


1 , 
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Table  8-9.  VJSWCM  subroutines. 


HEAD  -  This  subroutine  sets  up  the  headings  and  labels  on  the 
data  plots 

RESET  -  This  subroutine  initializes  data  values  at  the  start  of 
the  program  and  at  the  beginning  of  each  new  month. 

STEADY  -  This  subroutine  calculates  the  radionuclide  concentrations 

in  the  steady  flowing  stream  in  the  absence  of  precipitation 
runoff 

RAIN  -  This  subroutine  calculates  the  radionuclide  concentrations 
in  the  runoff  water. 

FLOW  -  This  subroutine  calculates  the  time-dependent  flow  rate  of 
the  runoff  water  into  the  stream 

SUM  -  This  subroutine  performs  trapezoidal  integration  of  the  time- 
dependent  radionuclide  concentrations  to  determine  the 
integrated  concentration 

TAIL  -  This  subroutine  writes  the  data  arrays  into  the  data  plot 
fi  1  e 
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Table  B- 10. 


W SW CM  stored  data  arrays. 


XB (31)  -  This  array  contains  the  area  of  each  watershed, 
see  Table  2-2 

XD( 31 )  -  This  array  contains  the  distance  from  the  center 
of  the  watershed  area  to  the  water  supply  point, 
see  Table  2-2. 

CPT( 12)  -  This  array  contains  the  monthly  data  on  the 
precipitation  threshold  for  surface  runoff, 
see  Table  2-6 

SR0PE(12)-  This  array  contains  the  monthly  data  on  the  ratio 
of  the  surface  runoff  to  the  excess  rainfall, 
see  Table  2-4 


FC ( 1 2 )  -  This  array  contains  the  monthly  data  on  the 
ground  water  areal  flow  rate,  see  Table  2-3 
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The  processing  of  the  disk  file  to  produce  data  plots  is 
performed  using  a  generalized  plotting  program  written  for  the 
Versaplot-07  System  ^  ^  )  A  FORTRAN  listing  of  this  program  is 
provided  in  Section  B-3.5.  It  should  be  noted  that  the  program  is  not 
part  of  the  WSWCM  program  and  was  not  developed  as  part  of  this  work. 

B-3.5  Computer  Program  Listings 

A  FORTRAN  listing  of  WSWCM  is  provided  in  Figure  B-9.  A  FORTRAN 
listing  of  the  figure  plotting  program  used  with  WSWCM  is  provided  in 
Figure  B-10. 
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PROGRAM  WSP 


A  PROGRAM  TO  ESTIMATE  ACTIVITY  CONCENTRATIONS  IN  A  STREAM 
FOLLOWIN'  A  CONTAMINATING  FALLOUT  EVENT.  CONTAMINATION 
DUE  TO  R<|"OFF  FOLLOWING  RAINS  IS  INCLUDED. 
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C  . 
C 
C 

C  . 

c 

c . 
c 

c . 
c 

c . 
c . 
c . 
c 

c . 
c 

c . 
c . 
c . 
c 

c . 
c . 
c 

c . 
c 

c . 
c . 
c . 
c 

c . 
c . 
c . 
c 
c 

c . 
c . 
c 

c . 

c . 
c 

c . 
c 

c . 
c 
c 
c 

c. 
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..THIS  PROGRAM  CREATED  FEBRUARY,  1982  BY  JIM  A.  ROBERTS  AT  SCIENCE 
..APPLICATIONS  INC.,  SCHAUMBURG.  ILLINOIS.  ITS  INTENDED  USE  IS  FOR 
. .SCALING  CALCULATIONS  CONCERNING  SMALL  WATERSHEDS  IN  EUROPE. 
..HYDROLOGY  DATA  SPECIFIC  TO  THOSE  WATERSHEDS  ARE  INCORPORATED  IN 
..THE  MODEL  UPON  WHICH  THIS  PROGRAM  IS  BASED.  USE  OF  THIS  PROGRAM 
..FOR  OTHER  CALCULATIONS  IS  DISCOURAGED. 


•FC  IS  THE  MONTHLY  AVERAGE  GROUNDWATER  FLOW  APPEARING  ON  THE 
SURFACE  OF  THE  WATERSHED  (L/S*KM**2). 

•SROPE  IS  THI  FRACTION  OF  EXCESS  PRECIPITATION  WHICH  RUNS  OFF 
ON  THE  SURFACE,  AVERAGED  BY  MONTH. 

.  C  P  T  IS  THE  MONTHLY  AVERAGE  PRECIPITATION  THRESHOLD  WHICH 
M J ST  BE  EXCEEDED  BEFORE  SURFACE  RUNOFF  BEGINS  (MM). 

.HLO.DAF  IS  THE  OATA  INPUT  FILE. 

.  WAT  E  R . DAT  IS  AN  OUTPUT  FILE  FOR  INTEGRATED  ACTIVITY  CONCENTRATIONS. 

.  P ..  .DAT  IS  AN  OUTPUT  FILE  FOR  PLOTTING  ACTIVITY  CONCENTRATIONS 
VS.  TIME  USING  PLOT  PROGRAM  CY  EGBERT. 

.R  IS  R/HR  AT  1M  ABOVE  GROUND  AT  HM  HR. 

H-l  HR  IS  THE  START  TIME  FOR  THIS  PROGRAM. 

-B  IS  AREA  0-'  THE  WATERSHED  <KM**?>. 

.0  IS  THE  DISTANCE  TO  THE  CENTER  OF  THE  WATERSHED  (KM). 

.TAU  IS  THE  AVERAGE  RESIDENCE  TIME  (HR)  FOR  GROUNDWATER  TO  FLOW  ON 
SURFACE  OF  STREAM  BEFORE  IT  PASSES  THE  WATER  SUPPLY  POINT  (WSP). 

•S3  IS  THE  AREA  OF  THE  STREAM  SURFACE  <KM**2>. 

,X4  IS  THE  MASS  (KG)  OF  100  MICRON  SEDIMENTS  WITH  DENSITY 
1.4  G  /CM*  *  3  IN  THE  STREAM  BED. 

.LTMAX  IS  THE  NUMBER  OF  HOURS  FROM  START  TO  END  OF  PERIOD 
0-  INTEREST  FOR  THIS  PROGRAM. 

.MOS  IS  THE  NUMBER.  OF  THE  MONTH  IN  WHICH  THE  EVENT  OCCURS. 

.N  IS  THE  NUMBER  OF  RAINS  FOR  WHICH  DATA  WILL  BE  ENTERED. 

•NNUCS  IS  THE  NUMBER  OF  P AR E NT- DAUGHT C R  NUCLIDE  PAIRS  PlUS 
THE  NUM3ER  OF  NUCLIDES  WHICH  WILL  BE  ENTERED  WITH  PARENT  'NONE1. 

.  L  T 1  IS  THE  TIME  (HR)  AT  WHICH  THE  RAIN  OCCURS. 

.  Xp0  IS  THE  AMOUNT  OF  RAIN  IN  A  DAY  (MM). 

.  Ap I  IS  THE  ANTECEDENT  PRECIPITATION  INDEX  (MM),  A  MEASURE  OF  THE 
DEGREE  TO  WHICH  the  GROUND  IS  SATURATED  FROM  PREVIOUS  RAINS. 

IT  IS  CALCULATED  USING  THE  SUMMATION  PRESENTED  BY  SilVERS. 

.MDR  IS  THE  NUMBER  OF  THE  MONTH  IN  WHICH  THE  RAIN  OCCURS. 

.SPO  IS  THE  AMOUNT  OF  PRECIPITATION  FROM  A  DAY'S  RAIN 
WHICH  APPEARS  AS  SURFACE  RUNOFF  « MM ) . 

.  L  H  R  IS  THE  NUMBER  OF  HOURS  IN  A  MONTH. 

.NAMEA  AND  NAME  B  -RE  ABBREVIATED  NAMES  FOR  THE 
PARENT  AND  DAUGHTER  NUCLIDES  RESPECTIVELY. 

■THALFA  AND  THALFS  ARE  THE  RADIOACTIVE  HALF-LIVES 
FOR  THE  PARENT  AND  DAUGHTER  IN  HOURS. 

•XKDA  AND  XKDS  ARE  THE  DISTRIBUTION  COEFFICIENTS 
FOR  PARENT  AND  DAUGHTER  WHICH  RELATE  THE  QUANTITY 
OF  THE  NUCLIDES  SORBED  ON  SEDIM'NTS  TO  THE  AMOUNT 
IN  WATER  SOLUTION  AT  EQUILIBRIUM  (CI/KG  /  CI/L). 

•GA  AND  GB  RELATE  THE  INITIAL  ACTIVITIES  OP  PARENT  AND 
DAUGHTER  TO  THE  R/HR  AT  H*1  HR.  UNITS  ARE  Cl/  R/HR. 

.  X  L  A  AND  XLB  ARE  THE  RADIOACTIVE  DECAY  CONSTANTS 
FDP.  PARENT  AND  DAUGHTER  (HR-1). 

.AA0  AND  AB0  ARE  INITIAL  ACTIVITIES  OF  PARENT  AND  DAUGHTER 


C  FALLING  AS  SEDIMENTS  INTO  THE  STREAM. 


Figure  B-9.  USWCM  program  listing. 
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C...CA  AND  CB  ARE  THE  T I  HE - DE P E ND E N T  ACTIVITY  CONCENTRATIONS 
C  Or  PARENT  AMD  DAUGHTER  IN  THE  STREAM  IN  THE  ABSENCE  OF 
C  RAIN.  UNITS  ARE  PCI/L. 

C...AWA  AND  AWB  ARE  ACTIVITY  CONCENTRATIONS  OF  PARENT  AND  DAUGHTER  IN 
C  STREAM  WATER  AT  THE  END  OF  A  MONTH. 

C...MON  IS  THE  NUMBER  OF  THE  CURRENT  MONTH. 

C...L  IS  A  COUNTER  FOR  THE  NUMBER  OF  TIME  STEPS. 

C...SUMFCA  AND  SUMFC3  ARE  THE  SUMS  OF  THE  PRODUCTS  OF  CONCENTRATIONS 
C  TIMES  FLOW  RATES  AT  A  GIVEN  TIME  FOR  THE  STEAD'.'  FLOWING  STREAM 

C  AID  ANY  RUNOFF  RAIN  WATER  FOR  PARENT  AND  DAUGHTER  RESPECTIVELY. 

C...SUHE  IS  THE  SUM  uF  THE  FLOW  RATES  FOR  THE  STEADY  FLOWING  STREAM 
C  AND  RUNOFE  RAIN  WATER  AT  A  GIVEN  TIME. 

C...NWSP  IS  THE  NUMBER  OF  THE  WATER  SUPPLY  POINT  Or  INTEREST  FOR  THIS  RUN. 
C 

c 

DIMENSION  F  C (  12). SROPE (  lOl.CPTi  I  2  > . LT 1 (  100  > . SSA<  1 300  )  , YB<  3 1  )  . 
lXDt31).SSB<  100U  )  . SRO<  100  )  . X!  1000  )  . VA(  1 000 ) . YB <  1000)  . YC <  1 000  ) 

COMMON /C 1/3. R.GA.GB.J. XKDA, XKDB 

COMMON /C  2 / V .XKA.XkB.XL1 ,  XL2.XI.  3.XL4.EL42.EL21  .  E  L  4  1  . 
IEL31.EL43.CA.CB.XLA.XLB  ,  AA0  ,  Alia  ,  AWA  .  AWE 
OAT A  XB/63. 6, 14.  .12. 6. 27. 1.12. 8, ’8. 1.12. 2, 30. 9. 13. 6. 7. 9, 

112.4.  16. 6. 9. 3. 2  1.  4, 33.  .17. 4. 10. 8,8.4.20. 5. 22. 1.8.7. 

214.  3,  19. 9. 31. 7. 45. 3,  15. 4. 9. 5.  '2. 5. 25. 8. 24. 7. 40. 2 / 

DATA  XD/5. 3, 1.4, 3. 2. 2. 7, 4.  ,3.2,3. 6. 2.  .2. 5, 3. .4. 5. 5. 3. 3. 2, 3. 9, 7. 2, 

3.  2.  2. 3. 2. 6. 3. 9. 3. 3. 3.1.3.  1.4.  .4. 6. 6. 4, 2. 6. 2.  ,  2  .  7  ,  4.2,4.  .4.2/ 

DATA  FC/3. 7. 5. 0.4. 9. 5. 0.3. 0.2. 3. 1.9. 

11.9.1.9.2.2.3.1.3.4/ 

DATA  SROPE/. 45.  .66.  .SB. .34, .11.  .10.  .13. 

! .10,  .12,. 16,  .40,  .42/ 

DATA  CPT/8..4..5..9..12..12..13.. 

13.  .12.  ,12.  ,10.  .9./ 

I N T E  G E  R  *  4  NAME  A . NAME  B 

OPEN! UNIT  =  I  . NAME  * • H20 . DAT '  , TYPE  * ' OLD  1  .READONLY  ) 

OPENt  UN  I  T  =  3  .  NAME  =  'WATER  .DAT  1  .  TYPE  =' NEW  ) 

OPE  N(  UNI  T  =  2  .  NAME  =  '  P  i  .  DAT  ,  TYPE --'NEW  ) 

WRITE!  3.333  ) 

333  FORMAT! 8X .'  INTEGRATED  CONCENTRATIONS  PC  l  *0AY /l '  , 7 M 
C...RFA0  IN  NUCLIDE  INDEPENDENT  DATA. 

READ!  I. 1  l'1. L  T  MAX  .  NWSP  .MOT  .  N  .  N'-UCS 
1  FORMAT! E10.3.IS.4!2> 

3  =  X  B ! NWSP  ) 

0 “ X 0 ! NWSP  ) 

T  AU  =  ( D  *  I  000 .  )/( 0.6*3600.  > 

S  B  =  0 . 0 1  *  8 
X  M  =  3  B  *  1 . 4 1  ♦  05 

C . . .READ  IT  RAIN  DATA  AND  CALCULATE  3RO  FOR  EACH  RAIN. 

DO  20  1=1 . N 

READ!  1  .7  > T  1  (  I  1.XP0.API  .MOR 
7  F ORMAT!  I  5 . 2E 10. 3 .  1  2  > 

APC “AP I -CPT ( MOR  ) 

I F  <  A P C . G T . 0.0  )APC  .  0 
SPO!  I  ) - SROP E ( MOR  )*<  XP0+APC  ) 

20  CONTINUE 

C... CLEAR  GROSS  ACTIVITY  ARRAY. 

DO  33  I  -  I  . I  000 
YC( 1  >'0.0 
30  CONTINUE 
L  H  P,  *  7  2  9 

C  .  .  . BEGIN  NUCL IDE  LOOP 
DO  300  1=1, NNUCS 

C. . . INITIALIZE  MON  FOR  EACH  PASS  THPOUGH  LOOP. 

MON  =  MOS 

C .  .  .READ  IN  NUCLIDE-DEPENDENT  INFORMATION  FOR  A  PARENT  DAUGHTER  PAIR. 

READ!  1  . 2  )  NAME  A , THAL  FA . XKDA . GA 


Figure  B-9.  WSWCM  program  listing  (cont). 
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2  FORMAT1A4.3E10.3) 

READ  i  1.2 >  NAME  B . THALEB.XKDB.GB 
XL  A  =  AL  OG ( 2.0  > /THAI  FA 
XL  B  =AL  OG  < 2 . 0  >  /  T HAL  F  B 
AA0--6B*R*GA 
ABU- BB*R*G8 

C...CALL  HEAD  TO  WRITE  HEADINGS  FOR  PLOTS 
CALL  HEAD'NAMEA.  NAME  B  .  L  T.'1AX  .  NWSP  ) 

L  =  3 

KN  TR - -  1 
AWA -0.0 
AWB -3 . 3 
C  A  ~  0 . 0 
C  B  ■  0 . 0 

C... BEGIN  OUTER  TIME  LOOP.  UHICH  CAUSES  THE  INNER  LOOP  TO  FIPST  GO  FROM 
C  1  TO  24  IN  STEPS  OF  1  HR.  THEN  FROM  27  TO  LTMAX  IN  STEPS  OF  3  HP. 

DO  230  I  P  E  P  -  1  .3.2 
ISTP  =  IREP 

I F ( I STP . EO . 1  >GO  TO  50 
I  T  F  -  2  7 
I Tl -LTMAn 
GO  TO  51 

50  ITF-I 

I TL -  24 

51  CONTINUE 

c..  .begin  inner  t:ve  loop. 

DO  1  TJ  ■  !lf  .  ITL  ,  I  S  'P 
C... CLEAR  SUMS  -OR  EaCH  TIME  STEP 
sumf'.:a  =  0.  o 
SUi>lFCB=0.0 
sum?  --P .  0 

LNOW=J/lHP 

i r  lnow.  l:  .fnt3  go  to  ill 

>  N  1  R  -•  L  NOW 
-••'0-FC>  KjN  ) 

C.  . *  CAu  L  RESET  ro  INITIALIZE  values  IF  IT'S  the  FIRST  TIME 
C  T -  ROUGH  AND  TO  RE-INITJALIZE  WHEN  IT'S  A  NFW  MONTH. 

CALL  RESET! F .NM.LNOW.LHR.F MO .B.TAU.XKDA.VF,  DB ' 

C... AT/ANCE  MON  FOR  NEXT  CALL  T0  RESET. 

MON-MON  *  1 

I  F  <  MJM  .  C.T  .  1  2  >M  TN^  I 
111  CONTINUE 

JT-J-; NOW'LHR 

C...J'  IS  NI.MBER  CF  HOURS  SINCE  START  OF  PROGRAM  0»  NEW  MONTH. 

c...o.L  s r e a u v  tq  calculate  conce n tpat ions  in  steady  stream  plow 


call  STEA'tv'JT; 

S 11 M  F  C  A  -  0  A  *  f 
SUMFCB-CB-T 
SUMF - F 

C . . .CALL  RAIN  TM  CALCULATE  THE  CONCENTRATIONS  AND  ASSOCIATED 
C  FLOW  RATES  INTO  ThF  STREAM  FOR  EmCH  RAIN  AT  TIME  J. 

DO  6  0  F  I  .  N 
;  MT1I  i.  > 

I  F  <  F  .  GT  .  1  1LT-Lt1(K)-LT1(f. -1  ) 

I F ( 0  .  GE  .  LT  1  ( F  > . AND . - J-LT  1  t k  )  )  .  LE . 96  ) 

ICAIl  RAIN!F,LT.lT1(‘  )  ,SR0(  F  )  . ,'i!  MFC  a. SUMF  0  6. SUMF  .XlA.XLB  ) 

60  CONTI NUE 

CONA-SUMF CA/SUMF 
CONB-SUMFOB/SUMF 

C...CDNA  AND  CONS  ARE  THE  T I  ME - D t P E N D E NT  ACTIVl’v  CONCENTRATIONS  Or 
C  PvRENT  AND  DAUGHTER  IN  THE  STREAM  INCLUDING  THE  CONTRIBUTIONS 
C  FROM  RAINS. 

C... SAMPLE  WRITES  FOR  THOSE  NOT  PLOTTING 


figure  !!-9.  WSWf.H  progr.im  listing  (cont). 
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C  WRITE! 3.3 (NAME A ,  J , CONA 

C  WRITE! 3.3  > NAME  8 , J , CONB 

3  FORMAT) IX, 'CONCENTRATION  OF',A4,'  AT  TIME', 15. 

I  1  HOURS  15  '  . E 12. 5.  '  PCl/L  '  ) 
l  =  l*\ 

X  (  L  > '  J 

IFICONA.lt. 1 .0E-21  >CONA= 1 .0E-21 
YA< L  '  =CONA 

IF ! CONB.LT. 1 ,0E -21  1C0NB*  t  .  0E-2 1 
V 3 ( L  >  =  CONB 
CONT=CONA*CONB 
YC <  L  >  =  YC ! L  )  ♦  CONT 
100  CONTINUE 

200  CONTINUE 

C ...  SET  UP  CALLS  TO  SUM  FOR  INTEGRATED  CONCENTRATIONS. 

SUMA=0.0 
SUME'0.0 
PSA!  l  1  =  0.0 

I  F  <  GA  .  EQ  .  <3 . 0  >G0  TO  201 
CALL  SUM! YA< 1 > . 1 .0. 24 ,SA1 > 

$  UMA  =  SA 1/24 . 

SSA<  1  )  *  SUMA 

C  WRITE! 3.298 ) NAME A , SUMA 

201  CONTINUE 

CALL  SUM( V8< 1  )  . 1 .0.24  .  SB  1  > 

SUMB=SB 1/24 . 

SS8  <  1  )  =  SUMB 

C  WR ITE < 3 .298 INAMEB , SUMB 

298  FORMAT! IX. ' INTEGRATED  CONCENTRATION  OF  .A4,'  AT  DAY  1  IS', 
1E12.5.  '  PC1*DAY/L  '  ) 

DO  91  L  DA  V  =  2 , L  TMAX / 2 4 
L  OC  =  24  +  !  ( LDAY-2  >*8  ) 

SSA! L  DAY >=0.0 
!F(GA.£Q.0.0) GO  TO  202 
CALL  SUM! YAILOC >.3.0.9. SA ) 

SA-SA/24 . 

SUMA  =  S  UMA»  SA 
SSA!  L  DAY  )  ■=  SUMA 

C  WRITE! 3,239 ) NAME A , LDAY . SUMA 

202  CONTINUE 

CALL  SUM! YB( LOC  )  , 3 .0. 9. SB ) 

SB  =  SB/24  . 

SUMB=SUMB*SB 
SSB! LDAY )• SUMB 

C  WRITE!  3.299  INAMEB. LDAY. SUMB 

299  FORMAT! IX. ' INTEGRATED  CONCENTRATION  OF  '.A4.'  AT  DAY'. 15, 

1 '  IS'  .E 12.5.  ‘  PCI ’DAY /L '  ) 

91  CONTINUE 

C . . .CALL  TAiL  TO  DO  FINAL  WRITES  TO  PLOT  FILES 
II -I 

CALL  TAIL! I .X  ,  YA , YB . YC , L  TMAX . 1  I  .NNUCS  ,  NAME A . NAME B . LDAY , SSA . SSB . 
1  NWS  P  I 

300  CONTINUE 

CLOSE i UN  I T* 2. D ISP*' SAVE  '  ) 

CLOSE ( U  N I T  = 1  ) 

CLOSE! UNIT=3.DISP= 'SAVE '  > 

STOP 
£  ND 
C 
C 

SUBROUTINE  RESET! F. XM . L NOW . L HR , F MO . B . TAU . XKD A . XKD B ) 

COMMON /C2 / V. XKA .XKB.XU.XL2.XL7.XL4.EL42.EL2I.EL41. 

I  EL  3  l  .EL43.CA.CB.XLA.XLB.AA0.A90.AWA.AWB 


Figure  8-9.  WSWCM  program  listing  (cont). 
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C. . .SUBROUTINE  TO  INITIALIZE  VALUES  AT  START  OF  PROGRAM  AND 
C  BEGINNING  OF  EACH  NEW  MONTH. 

C 

IF( LNOW.LT. 1 >G0  TO  10 
AWA=CA*V/1 . 0E ♦  l  2 
AW8  =  CB  * V /  1 . 0E ♦  1  ? 

C...AWA  AND  AWB  ARE  ACTIVITIES  OF  PARENT  AND  DAUGHTER  IN  STREAM 
C  WATER  AS  COMPUTED  AT  THE  LAST  TIME  STEP  IN  PREVIOUS  MONTH. 
E1HR=-XL1*LHR 
IF(E1HR.LT.'50.  >E 1  HR  =  -50 . 

EIHP=EXP(EIHR) 

E2HR  =  -'/L2*LHR 
:f(e:hr.lt.-50.  >e2hr=-50. 

E2HR--EXP<  E2HR  > 

A  B  3  =  XL B * AA0 / < XL2-XL 1  >  >*<  E 1HR-E2HR  >+AB0*E2HR 
AA0 - AA0*  E l HR 

C...AA3  AMD  AB0  SET  TO  ACTIVITIES  OF  PARENT  AND  DAUGHTER  IN  SEDIMENTS 
C  AFTER  ONE  MONTH  OF  DECAY  AND  INGROWTH. 

10  CON  F I NUE 
F  =FMC*B 
V  =  F  *TAU*3600 . 

XKA=ALOG(  1  ,0+V/( XKDA«XM>  l/TAU 
XKB  =  ALOG(  1  . 0  +  V/ ( XKDB*XM>  l/TAU 
XL  1  * XLA  +  XKA 
XL2=XLB+XKB 
XL3-XLA+1 .0/TAU 
XL4=XLB+1 .0/TAU 
EL42=XL4-XL2 
EL21 -XL2-XL 1 
EL  4 1 =XL4-XL 1 
E  L  3 1 =XL  3-XL 1 
EL43=X14-XL3 
RETURN 
END 
C 
C 

SUBROUTINE  STEADY! JT  ) 

COMMON /C2 / V, XKA. XX B .XL 1.XL2.XL3.XL4.EL42.EL21.EL41. 
IEi.31.EL43.CA.CB.XLA.XLB.  AA0 .  AB0 .  AWA ,  AWB 
C 

C. . .SUBROUTINE  TO  CALCULATE  CONCENTRATIONS  IN  STEADY  FLOWING  STREAM. 

C 

E 1 -< -XL  1 *  JT  ) 

I T ( E  t .LT.-50.  ) E 1  - - 50 . 

E  1  -EVP ( E  1  1 
EC*!  -  XL  2* JT  ) 

I r  < E2.LT. -50.  ) E 2  =  - 50 . 

E  2  "  E  X  P  (  E  2  > 

E  3  =  < -  X  L  3  * J  T  ) 

IFIE3.LT. -50.  IE  3  =  -50 . 

E  3  =  E  ;P(  E  3  ) 

E  4  - ( -  XL  4  * JT  ) 

IFIE4.LT. -50  .  > E 4 - -59 . 

E  4  =  E  X  P  (  E  4  ) 

TERM! =XLB*XKB*AA0/EL21 
TERM2-E 1 /EL4 1 -E2/EL42 
TERMl'-yKB*AB0*E2/EL42 
TERM4  =  XLB*XXA*AA0/El  31 
TERM5-E1/EL41-E3/EL43 
TE  PM6  =  XL  B *AWA*E  3 / E  L  43 
TERM ' - ( 1 . 0/EL41-1 .0/EL42) 

I ERMB  =  XXB  AB0/EL  42 
T E P M 3 - (  1  .0/CL  41 -1 .0/EL43  ) 

TERM10=XLB*AWA/EL  4  3 


Figure  B-9.  WSWCM  program  listing  (cont). 


B  -  5 1 


o  o 


E  LC8  =AWB- < TERM 1 *TERM7* TERMS *TERM4*TERM9 +  TERM10  > 
CB-(1.0E+12/V>*<  TERM1  *TE  RM2-*  TE  RM3  ♦  TE  RM4  *TE  RM5  *TE  RM6 +  E  L  C  B  *E  4  ) 
CAM  1  .0E  +  1 2/V >*< < XKA ’ AA0 /  E  L  3  1  )*( E 1 -E3  >*AWA*E3 ) 

RETURN 

END 


SUBROUTINE  RA I N <  K , L T , L T 1 . SRO , SUMF CA , SUMF CB . SUMF . X L A . XL 8  > 

C 

C. . .SUBROUTINE  TO  CALCULATE  T I  ME - DE P E ND E NT  CONCENTRATIONS  OF  PARENT  AND 
C  DAUGHTER  IN  RUNOFF  RAIN  WATER. 

C 

COMMON  /C 1 / B , R , GA .GB.J.XKDA.XXDB 
DIMENSION  AA( 100 >.A8< 100 > 

V l =  SRO*  B* 1 .flE  +  06 

C...V1  IS  The  VOLUME  OF  WATER  <  L  >  ASSOCIATED  WITH  A  GIVEN  PAIN 
C  =  MM* ( M/ 1 0*  *  3  MM  )  * ( AR  E A  IN  XM* * 2  ) * <  1 0* * 6  M* * 2 / KM* * 2 > *  1 0* *  3  L/M**3 
XM1~. 99*8*1. 4  E  *05 

C...XM1  IS  the  MASS  OF  100  MICRON  PARTICLES  WITH  A  DENSITY  OF  1.4  G/CM**3 
C  COVERING  THE  SURFACE  OF  THE  WATERSHED  OUTSIDE  THE  STREAM  BED. 

IF ( K  .GT  .  1  >GO  TO  7  1 
AA< X  >-R*0 . 99*B*GA 
A6< K >=R*0.99*B*GB 

C .  .  . AA (  1  )  AND  A3(  1  )  ARE  THE  ACTIVITIES  OF  PARENT  AND  DAUGH.ER  ASSOCIATED 
C  WITH  PARTICLES  AT  THE  START  OF  THE  PROGRAM. 

71  CONT.NUE 

C...SET  A I  A  AND  A I B  TO  ACTIVITIES  LEFT  IN  SEDIMENT  AFTER  LAST  RAIN. 

A  1  A  =  A  A  (  K  ) 

A 1 B  =  AB<  X  . 

ELA--XLA*'_T 

I F  (  E .A . L  T . -50.  )  E  L  A- -50 . 

E IA-C  XP ( E  l A  ) 

ELB  -  X  L  B  *  l  T 

1FIElS.lt. -50.  >  E  L 6= -  50 . 

ELB  =  EaP( ElB  > 

A1B--1XLB/(XLB-XIA)1*A1A*'ELA-EIB)*A1B*ELB 
A  I A  =  UA'iLA 

C...A1A  AND  A 1 B  ARE  NOW  DECAYED  TO  time  OF  PRESENT  RAIN. 

C...E1A  AND  E1B  ARE  AfllVITIES  LEACHED  INTO  RUNOFF  WATER. 

E 1  A  -  A  1 A / <  1 +  X  X  DA  *  XM 1  -VI) 

E  1 B -  A  1 B / t  1 +XF  DB*XM1  'Vi  ) 

C...JT!  is  time  FROM  start  OF  THIS  RAIN  TO  PRESENT  TIME  STEP. 

0  T  1  '  D  -  L  T  1 

E 1  LA- - < L A*( UT 1  I 

I  F  (  E  I  LA  .  l  T  .  -50  .  IE  1L  “---50  . 

E  1LA-EXP(  E  1 L  A  ) 

E  1  L  £  '  -  X  L  B  *  (  0  T  I  - 
I  F  (  E  I  LB  .  lT  .  -50.  TUB  -  -50  . 

EILB'EXPIEILB) 

C...CIA  AND  C.B  APE  ACTIVITY  CONCENTRATIONS  IN  RUNOFF  WATER  AT  TIME  0. 

C  I  B  =  <  1  .01*1  2 /  VI  1*1  i  <  l  B  *  E  !A/l  X  L  B  -  X  L  A  '  >  *  <  E  UA-EUB  >  *  £  i  B  *  E  1  L  B  > 

C 1 A  = {  1 . 0E  +  1 2/V l  )  *  E 1  A  *  E  U  A 

C...SET  AA  AND  AB  TO  ACTIVITIES  LEFT  IN  SEDIMENT  AFTER  THIS  RAIN. 

AA<  X  * !  )  -  A  I  A- E 1  A 
A  E  <  L  •  1  )  A  I  3  - 1  1  B 

C . . .CALL  FLOW  TJ  CALCULATE  FLOW  RATF  INTO  STREAM  FOR  RUNOFF 
C  FROM  THIS  PTIN  AT  TIME  J. 

CALL  FLOW.  SRO.F  1  ..IT  l  ) 

C  WR IT:<  3. 72 )JT1 .C 1A.C IB  ,F  1 

72  FORMAT!  IX.  •  J-LT1  -  .13.'  CIA  =-,E12.5.'  C1B  ■=  '  .  E  1  2 . 5  .  ‘  F  1  =',E12.5) 
C...ADD  CONTRIBUTIONS  FROM  THIS  RAIN  TO  PREVIOUS  SUMS 

S  UMF  CB-SUMFCB*..  IB*C  I 
SUMFCA-SUMFCA*C I A*F i 
SUMF -SUMF ♦ F  1 


f  i  (jure  B-9,  WSWCM  pro'jtuiP  listing  fcont). 


oooooo  oooo 


RETURN 

END 

C 

c 

SUBROUTINE  FL0WCSR0.F1 ,0T1) 

...SUBROUTINE  TO  CALCULATE  THE  T I ME - D E P E NDE NT  FLOW  RATE  INTO 
THE  STREAM  FOR  CONTAMINATED  RUNOFF  RAIN  WATER. 

EP3=-.301 1  * JT  1 
IFCEP3.LT. -50.  >  EP  3  =  -50 . 

E  P  3  =  E  X  P  <  E  P  3  ) 

EP03-- ,0343*OT1 
I F ( EF03 . LT . -50 . >EP03=-50. 

EP03=EXP< EP03  ) 

F  l=-5RO* 1  1  . 3 9 * < EP3-EP03  > 

RETURN 

ENO 


SUBROUTINE  SUM< V . DX , NO . S  ) 

..TRAPEZOIDAL  INTEGRATION  ROUTINE  FOR  COMPUTING  INTEGRATED 
..CONCENTRATIONS  IN  PCI*DAY/L 

DIMENSION  Y<  1  ) 

C...  V  ARRAY  HAS  Fill  VALUES 

C...  UX  =  EQUAL  DISTANCE  DELTA  X  VALUES  (TIMESTEPS) 

S  -0 . 0 

DO  1  1*2. NO 

1  S  =  S*(  C  Y (  I  )  +  Y(  1-1  >  )  / 2 . 0  )  *  DX 
RETURN 
END 

C 

c 

SUBROUTINE  HEAD  < NAME A . NAME B . L TMAX , NWSP  I 
C 

C.. .ROUTINE  TO  SET  UP  HEADINGS  ON  PLOTS 
C 

I NTE  GE  R*4  NAME  A , NAME  B 
WR ITE  <  2 , 8 i NAME  B . NWSP 

8  FORMATC 3X . • 70 • . UX.A4 , IX . 'WATER  CONTAMINATION  WSP  #'.I2> 
WR ITE ( 2 , 9  )  NAME  A 

9  FORMATC 3X.  '70'  ,23X,  PARENT  IS  .A4> 

WR ITE l 2 . 1 0 (NAME  B 


10 

F  0 P M  \ T  t  3 X .  ' 35 '  . A4 , 

IX.' CONCENTRATION 

IN  WATER  ( PCI/L  ) 

WP  I  TI  <2,11  > 

1 1 

F  ORMATC 3X .  '  31 T I  ME 

SINCE 

INITIAL 

FALLOUT  ( HR ) ' ) 

WR I TE ( 2  .  1 2  > 

12 

"ORMATC 9  X ,  ' 0 '  , T  y .  * 
WP I TE ! 2  .  1  3  > 

1  .  9  X 

. ' 2 ' . 9X. 

'  3  '  ) 

1  3 

" ORM-vT C  4 X  .  '  2  I 

WR 1 T" <  2 . 1 7  > 

1  7 

FORMA  T (  1 X .  ■ 3 ■  .  1 X  .  • 
WR I TE ( 2 .  1 4  >L  TMAX 

10.  '  . 

1  X  .  '  1  3  .  ' 

.IX.' 

.650'  ) 

1  4 

form-,  rc  1  X  .  '  0 . 0  '  .  IX 

WR I TE ( 2  .  1 5 > 

.  I  5 

15 

FORMATC  IX.  ' 9 . E  *06 ‘ 

.IX.' 

1 . 0E  <-00  ' 

) 

RETu^N 

END 

C 

C 

SUBROUTINE  TA I L ( L . X . YA . YB . YC . LTMAX . I  I .NNUCS. 
1  NAME  A .NAME B , L DAY . SSA . SSB , NWSP  ) 


figure  B-‘l.  WBWCM  proqr.im  listing  (cont). 


B-r»B 


..ROUTINE  TO  WRITE  DATA  ARRAYS  TO  PLOT  FILE 

?MTr 'f iCN  x<  1  1  ’  YAl  1  ' . VB<  I  ),YC(  i  i.SSAI 1  )  , SSB (  1  > 
!NTEoER*4  NAME A , NAME  B 
WRITE; 3,30)  NAMEA . NAME  B 
30  FORMAT! /, 10X,  ' DAY '.  7X.A4 . 10X  A4  /> 

00  4.1  K=1  ,  L  DAY 

WRITl(  3.032  ) k , SSA! X ),SS8(K> 

23  format; 8X, 15. •  . . . ' . 1 P . E 1 0 . 3 . '  ...'.£10  3) 

40  FONT  I HUE  J 

WR I TE ( 2 . 4  >L 


'*  R  1  T  E  '  2  ,  *  1 
WR I TE  (  2 , * ) 
W  R  I  T  E  i  2  ’.  5  > 
WR I TE ( 2 . *  ) 
WR ITC ( 2 . •  ) 
5  FORMAT; 15. 
•i  "CRM  \T(  15. 


*  >  <  X  (  L  L  > ,  L  L  -  I  ,  L  ) 

*  '<  YA(  LL  )  ,LL  “I  ,  L  ) 

5  >L 

*  X  X( LL  > ,LL*  1  ,  L  > 
"MVB(LL).LL'l.L) 

5 .  '  1  4  5  ‘  ) 

5  .  '  1  4  1  ■  ) 


WRITE; 2.6> 

6  F  ORMATi  IX  .  •  (*■  ) 

I  *  I  I 

IF(  I  .  LT  .  NNUCS  (PETURN 

■F!WRITE?2U!9”nws0  NUClIDES-  SET  Up  PL0T  0F  GROSS  ACTIVITY 
19  ,F°R^(T!,3>;  j  ^  ^  •  1  4  X  .  'TOTAL  UATER  CONTAMINATION  WSP  #'.12 

21  FORMAT; 3X ,  70 ',  15X .  SUMMED  OVER  ALL  NUCLIDES') 

WR  I  i  t  (  2  .  .T  J  ) 

22  fORfUTi^y .  ' 31GP0SS  ACTIVITY  IN  WATER  <  PC  I /L  >  '  > 

WR  I  i  l  • r ,  1  i  ) 

11  WRITE(R3X;  3iTIME  S!NCE  :nIT1al  FALLOUT  (HR)') 

12  FORM  at;  9X  .  '  0  '  ,  9V  .  ■  1  '  Qy  •  •>  •  g»  .3  , 

WRITE  t  2.23 ) 

.’3  F  0RM»T;  i  x  .  '  1  '  > 

write; 2, 1 7  ) 

i?  FORMAT;  IX.  ' 3'  . IX.  '  10.  '  .IX.  '  13.  '  . IX,  '  .650'  ) 

WRITE; 2.14  1LTMAX 

4  r 0 D m r <  ix.  'p.er  .  1  x ,  1 5  > 

WPiT£.  2,25  t 

15  I  OPM'.T;  IX  .  9. E*  3 •  » 

WR  IT;;  3.5  >  L 

W  R  I  *  ;  .  3  ,  *  ■  <  V  :  .  .  '  .  L  L  -  1  .  L 
wpI'F-  2.*  VC 


Figure  B-9.  WSWf.M  pro.-r.r  liMm- 


A  PLOT  OF  V  VS.  X  ON  L I N ( LOG  1 0 ) - L I N ( L OG 1 Z )  FOR 

=  vr  ?  >;n  ,s  /'.no  on  a  rot/. i  l y  coo. 

■  :!  n. 

0  "  J .  L  PO’-'(  •*«  > ,  VHARK  (  l S  > , T V ?  C  <  5  > . ! P TV? E  <  7  > 

.  - 1  ■  i  > .  v  1 1  .  v.  l;  if!  !i ).  v ;  ::r  (/.  > 

.  .  ....  {■:).<.  .  i  / '  .  !/> 

v  ....  ’ !  .  •  ?  / 

.  *  i  . ; .  / .  \  " 

■  - .  .  :  !  i  i  ' 

: . .  .  .  •  ■')  t  •  t : 

-,l  •'*,  ;r  ‘  \v:~’ '  rAC h  vay: 


,i  i 

r  .)  I  ; 

.  •.  »  «  .  1  .  •  .  i  .  S,  \  .  7S  .  ?  .  5  ,  ?.  .  T» ,  f-  .  r' ,  S.  5 


.  .  v  i.  *'  r  ?: v  *  t  r  giv.d 


;•  v  ~  \  p l,\ ■  ?;  on  rn.v 

I.,  T  Y  V  ■  J  i  ' :  :J 7  Tr:‘  > 

v-  ■’  i  »..v  \  > ■  :  •  . .  v* ;  ;  V;*:::  cls  too 


t'l’'-1  *: \  •••  ’ . y.  ^ 


t  i  v:.- :/  -.?T  ViA'i v.  ( ^  a  * :  ?-*> 

;  *  •-  r  •  .  •; i  t • 


7rt<  f  '•  S  ) 


lng  ;  in';:  <•  .  *:  i.  .  1  c 


(  '  o  i 


X  .  A  ' 


■c^ro  (  :  t 


■  f  *)  f  ‘TT 


' 'JOT*:  •  I"  T 1 I  S  T  >  !  >  TvrE  *cunvr  ♦W?’' 3  V  > 
:*  rs  ;  i,  ,nj  V  n  a  \  \ 


i  r  r  r  t  n  •*;  ' 


Figure  B-10.  Data  plot  program  listing. 
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TME  V  DATA  LIST  DIRECTED  "  "  "  . 

'"■P  i:  l  '■■  !  I  E  LINER  NT.'  T  TIN'S) 

'  *  *s.\  •  •  \  •.*.  .  r>.;;  .  '  7  .M.i.oV#  V  t  yX-  AN'^TH  .  T  C':'"P’_cTT  IN.'JT  HE 


;  >r 


.  or  d/.t >  rot  ,t>. 

o  "i  \  It  :  ’  St  :  I  L  .  I  '  ,  TV  ?E  =  '  OLD  ’  ) 

•  /  .  •  . '  >  i '  ::e  i 

■ .  •  i  . 


]  - 

1  •  V  ,  ) 

1  .  >  . 

ir'.  .  • . 

i  t  .  .  IL'.; 

'  i  v  ;  I'  l  .t 

>  *.  <  :  i  m  ;•  -  1 

)  .*  1  i. 

.  in'.iDv 

■  ■)  .  .  F.  0  .  "  '  "0  TO  TOl 

n  v 

If. 

■.  n 

(  7  ,  i 

?  •  •  •  •  *  1  "  1 

M..  ••  -.7  ) 

r .  i  j 

'  7  ,  ; 

•  <  v  v/’-.y 

factr 

.  1  /  ..H/.' 

.  r  •  . '  -  C> f. 6  )  N ?  T  .  I  T  Y P  C  ,  !  H  I  s T  . 

nr  n 

M  1 

\  7  s  , 

\ 

•  M  w  T 
M  V  '  j 

.  i  •  1 4  * .  p  t  ) 

.  [  - : .  i  -  .  > 

-  ; 

>■  1 

!  >  f  i  _  '1  7 

’I  '  \  *  **  ;j  ' 

[r  ij-„-  i 

A D I '  f  •  C  J " V  E  Mywnra  '.IT) 

r* 

{  :  \  ■ 

•  <  * 

.  IT.-) 

*  .  ■  > 

yi.  . 

..  > 

;  T, 

v'  ::  i.ogi~  •  r  iv  0\  iv  :  ■  \ot  =  3 

'■■(!)  -,‘,!.Qr 
>  .  v  ?  :  7  • 

{'!(])) 

■  '  mm 

.  1  ;  T'  10  ■ 

1 

,  -  \  .  C  ■  V' 

t  . ./  it  . 


j\-rc  ,  7  •  V  r:;.,"  -  r-  p  •  r-  C.-vL  E  TO  INT-A-? 

'  •  -i  •'  \  v  •  I r  I  ’  <  :  ;  '-0T  *  A . 


(-  i  i .  i  -<!>.•-  ''-p-:  ,  i ,  y  •  •  .  "i 

F  :  I  !  .  ■  .  ■  ■ 


Figure  B- 10.  Data  plot  program  listing  (cont). 
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YXTV-0. T 

v  .  Y  VrV  5  ) 

.  ,  T  , 

.  .,<  I- n  U  TiX.<  VtlEV- .?!?  >.  .2  ,  16,21.0.-1  ) 

.  .  .  ■ .  sv  u- 1 )  u . :  v.- .  1 ,  ( v:  ■./- .  <■■2  > , .  l?  . t at -  l  > 

u;;  i  i 

■  •  •  .  i  ,  i. :  o'  i )  > 

i.  •  .  .  ;i  (  . I V  I  ) .  (  YK'W-  .3  ) .  .  '~9.  .  P?WR  .fl.0,-1) 

'(.i  'i  o  <;•  \ 

.  ■  I  -  1  >  OM  I TX 

:  v.  v;:: .jCTT  > . .  is . )  6  .a.fi.  - 1  > 

v: i  ;/?.r 

:  l  l.  '  I  1  S  1  :5  ,  NAMEY  .IT.. 9. LMGX  > 


2  ,  17  .  c~.  ■  >CaU  AH  I  ■(  0.  .0.  .  f.’.VlEY  .  LiiGY  .  VAX  .9.0.  ,  Y<  N">T*1  >.  YINPT  +  2  >  > 

T - . *  7.1)  CO  TO  1 ' 

1  ; j  1-2,1 3 

-.5  j-u  oio-'ui  f l o.f.T <  i ) 

..11  i-'.Ii  IS  ./.RPEK  THAI  VAX. 

.  ’<  1  '-7  ;• 

l  \ - 1  - 1 

l  •  ‘i  ;  >-  L  '■;•■(  l-i  )  + 1 
>  u:iiT 

t  ,7  2".-  .  '  YA  ;  1  GO  TO  <1 

•  n  •  •  -  •  •••'.* 

.  co . 

■'  i  c  "'  rr.; : 

;  1  z  r-  * ; .  2 

■  >"Xax 
. ■  :  A  2  J  0  J  fSD  Y  AXIS. 

■  •  i  nor.  -ivj.yax.S) 

-  '  l  572  y.  t  :<r>.'.0. 0.3. 1 . 16.0. jf. -2  > 

'  )  -  3  :  =  1  .  t.  I  mV 

7  SY.m  .*•_  !  (  X  \ZU-  .  12*<  I  REP-  1  )  ) ,  <  I  -1  >«0^1T.  .2. 15.0.0.-1  > 

. .  :t  <  x:i ,g*(  1  re p - 1  >  > . '  1- 1  > -uhit.  . 12 . ten . j.a, - 1  > 

.  7-.  12*  !  -  1  >  UNIT 
_  •  -'LCATi  L  ?  0  ’  <  I  >  > 

•:<  1  3*.  6*<  irep-i  >  i.yhyt.  .aq, pour, 0. 0,-1  > 

- L  .  IV  )  10  TO  43 

13  '3  : .  -  1  , i 

1  .Vi-'. m >*<  i - 1  >*unit 

v.  L  S',  7  '.i  (  '.J-.r;r:,(,-.C75*t  IT.EP-1  n.TVNARK.  .  15,  '5. 3. 0,-1  > 


'  7  .  V  I  0  I  2  5  .  NAME  v  ,  9X? .  rr ,  l  NG  Y  ) 

. ■  "■ .  van  ♦ .  5 . .  1 5 ,  n;.  t  :  .  r  .r ,  l  ijg  n 
I.'.,,  i  .  VAX  '  .  IT  .  .  12S  .  .C  .13  3r2> 


TH 


I  !  !  !  1  !  !  !  !  !  7???????????************ 


!• -  iro.tOY.FO.  ?  >C0  TO  pft.o 

>  '  •  1  -  i  i  i  1  1 


1  . LMASX2 . LKAGX3 


Figure  B-10. 


Data  plot  program  listing  (cont). 
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IF (  ICRIDX.EQ.0)GO  TO  501 
KX  X  =  <T 

■:  •  i 

v'  (  i )  =  u n  I T x 

■ ;  <  i  x  xx< 1 >  =  i . 

to  5  I  i 
m  x-i 

V  -  f 

,  <  !  >-'vN;TX*(  LIMX-1  ) 

:  ■  <  i  :< .  to.  j  >'::<<  \  >m  t  imx-  l .  ) 
•  10  '■  ~  <  I  i R  1 .1 V  .  L 0  .  O  > GO  TO  511 


•.  .'I  1  1  - X  I  T 
i ,  i  < .  5 .1 .  ; :  y  y  :  i  >  - 1 . 

Y  T  r 

11  .  .  .  =  1 

■  ■  .  i  :  - <  l  r;y- 1  i’um  it 

If ;  IV . to.  ) VV<  t  )  =  <L I MV- 1 .  > 

•  20  tV  LL  5  U0;.J.f>.J'..0.KX*(LIMX-l  >  *K  X  X  .  X  X  .  K  Y  *  (  L  I  MY-  1  >  +  KYY , YY  ,  LMASK1  1 

!  T  (  1  -;x  i  0  X  . E  .  2  100  TO  52  1 

■  ■  ;  ' '  y 

-■  3  .1  if’.  .11  MX 
. 0  V  :  2  0-1.9 

31  ( <  1-1  l-UXITXM ALOG10< J*1 .  ) -ALOG 1 0( 0 +0 .  )  ) 

.  >  «...  r»M-t 

■',>  ro  3  jot 

-l  i  :  :  ,=  iox.lt. 2  >  no  to  52.; 


5;  5  .  LIMX 

1  !  ■  1-3  '-0  liT'.-*(AlOG10<  2.  )-ALOG10(  1  .  >  ) 

•  1  NI  TM  <  ALOGlfM  3  .  )  -  AL  OG  1 0<  2  .  )> 

<  1  •  '-1  )-  J  I !  T  !*<  ALOGlfM  6.  >-A'.OG10<  3.  >  ) 
i  :  -  •  -  1  !  I  TX  ■*(  ALOG  10(  10.  )-AlOG10<6.  )  > 

■;  ;■  .7 


,  .  ' - . T 1 ~  X ” ( L I  MX  -  1  ) 

.<;•>.  i)-i.  / 1  n  r  0  x 
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Figure  8-10.  Data  plot  program  listing  (cont). 
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Figure  B- 10.  Data  plot  program  listing  (cont) 
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Figure  8-10.  Data  plot  program  listing  (cont). 
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SECTION  B-4 


SAMPLE  WSWCM  PROBLEM 


B-4. 1  Problem  Statement 


The  watershed  of  Water  Supply  Point  Number  8  is  contaminated  on 
July  4  with  a  fallout  intensity  of  1  R/hr  at  H+l  hour.  Determine  the 
water  contamination  that  occurs  out  to  August  14.  Use  the 
precipitation  data  given  in  Table  B - 1 1 . 

B-4. 2  WSWCM  Input 

The  user  input  requirements  for  the  computer  program  WSWCM  were 
described  earlier  in  Section  B-3.2.  Figure  B-ll  illustrates  the 
sample  problem  input  by  giving  examples  of  the  basic  problem  data,  the 
precipitation  data  for  one  of  the  rain  events,  and  the  radionuclide 
data  for  one  of  the  parent-daughter  radionuclide  pairs.  The  sources 
of  the  data  values  used  for  the  sample  problem  are  discussed  below. 
Figure  B - 1 2  shows  the  actual  data  input  for  the  sample  problem. 

B-4. 2.1  Basic  Problem  Data 

i 

i  The  problem  statement  specifies  that  the  Fallout  Exposure  Rate  of 

|  the  Watershed  is  1.0  R/hr  at  H+l  hour,  R  -  1.0. 

The  Problem  Simulation  Time  must  cover  the  period  of  July  4  to 
August  14;  this  period  is  about  999  hours,  LTMAX  =  999. 

The  problem  statement  gave  the  Identification  Number  of  the  Water 
Supply  Point,  NWSP  =  8. 

The  month  in  which  the  Problem  Starts  is  July,  MOS  -  7. 
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pie  problem  input. 
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The  Number  of  Rains  during  the  Problem  Simulation  Time,  from  July 
4  to  August  14,  is  determined  from  Table  B-ll,  N  =  21. 

The  Number  of  Parent-Daughter  Radionuclide  Pairs  included  in 
Table  B-l  is  22,  NNUCS  =  22. 

B-4.2.2  Precipitation  Data 

The  precipitation  data  to  be  used  in  the  sample  problem  is  given 
in  Table  B-ll.  Between  July  4  and  August  14  there  are  21  rain  events. 
The  procedure  for  developing  the  precipitation  data  for  WSWCM  can 
be  shown  by  considering  the  first  rain  event. 

The  first  rain  event  occurs  on  July  6,  this  is  2  days  after  the 
problem  starts;  the  Time  at  which  the  1st  Rain  Occurs  is  thus  48 
hours,  LT1 ( 1 )  =  48. 

As  seen  in  Table  B-ll,  the  Amount  of  Precipitation  in  the  1st 
Rain  is  1.1  mm,  XPO  =  1.1. 

The  Antecedent  Precipitation  Index  for  the  1st  Rain  is  calculated 
using  the  method  given  in  Section  B-2.4.2.  An  example  of  this 
calculation  is  given  in  Table  B-12,  API  =  21. 

The  Month  in  which  the  1st  Rain  Occurs  is  July,  MOR  =  7. 

B-4.2.3  Radionuclide  Data 


The  set  of  data  presented  in  Taule  B-l  includes  22 
parent-daughter  radionuclide  pairs.  The  procedure  for  developing  the 
radionuclide  data  for  WSWCM  can  be  shown  by  considering  the 
radionuclide  Te - 1 3 1  in  (from  the  Te-131m,  1-131  pair)  as  the  first 
parent  radionucl ide. 

The  Identification  Name  of  the  1st  Parent  Radionuclide  is  an 
abbreviation  of  Te-131m,  NAMBA  =  T31M. 


Table  B-12.  Antecedent  precipitation  index  calculation. 
(Example  for  July  6) 
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The  Radioactive  Half-Life  of  the  1st  Parent  Radionuclide  is  the 
half-life  of  Te- 1 3 1  in ;  the  value  can  be  found  in  numerous  nuclear  data 
reference  books,  THALFA  =  30. 

The  Distribution  Coefficient  ot  the  1st  Parent  Radionuclide  is 
the  Kd  value  for  Te  given  in  Table  B-7,  XKDA  -  100. 

The  Normalized  Ground  Concentration  of  the  1st  Parent 
Radionuclide  is  the  value  given  in  Table  B-l  for  Te-131m,  GA  =  13. 

13-4,3  WSWCM  Output 

The  output  provided  by  the  computer  program  WSWCM  for  the  sample 
problem  includes: 

o  33  plots,  one  for  each  parent-daughter  radionuclide  pair, 
of  the  radionuclide  concentration  in  water  as  a  function  of 
time  (Figures  B-l 3  to  ii 

o  1  i lot  of  the  composite,  or  total,  radionuclide  concentration 
in  water  as  a  function  of  time  (Figure  B - 3  b ) ,  and 

o  1  tabular  listing  of  the  1 tme- i ntegrated  radionuclide 

concentration  in  water  as  i  function  of  true  for  the  parent  - 
daughter  radionuclide  pairs  (Table  B-13). 
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Figure  B-24.  Zr-95,  Nb-95  viater  contamination 
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Sr-90,  V-90  water  c on tani nation 
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Figure  B - 2 0 .  Ru-103,  Rh- 103  water  contamination 
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Fiqure  B-30.  Cs-137  water  contamination 


RH06  WATER  CONIWKIMU  ! ON  KSf  »  8 


Figure  B -32.  Ru-106,  Rh-106  wa ter  contamination 
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Figure  B-34.  1-134  water  contamination 
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